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Gamma-Ray Spectroscopy 


By S. C. Curran, F.RB.S. 
Department of Natural Philosophy, The University, Glasgow 


§1. INTRODUCTORY 


Many major advances have taken place in the field of gamma-ray 
spectroscopy in the last seven years. During this period, on the 
experimental side, we have seen the introduction of the proportional 
counter and the scintillation spectrometer, in a variety of forms, together 
_ with what might be termed the re-introduction of the crystal spectrometer 
in the form of the large bent crystal. At the same time, and even more 
significantly perhaps from the point of view of ultimate progress on a 
broad front, radical changes have been made in the more orthodox 
techniques of magnetic beta and gamma spectroscopy. Among the 
notable achievements here we have the successful solution of the problem 
of double focusing which has now proved very capable of application to 
work at high resolution and recent developments along the lines of lens 
spectroscopy are no less notable. In parallel with these advances the 
introduction of new methods of measuring magnetic field strength, for 
example magnetic resonance and nuclear induction, constitutes a really 
powerful aid to the achievement of improved standards of accuracy in 
energy measurement. Both in magnetic spectroscopy and in the bent 
crystal method new standards of great accuracy have been set. In brief 
the tools now available for the analysis of gamma radiations at low and 
at high resolutions are very considerably increased in number and in 
power. Preliminary studies with the new low resolution instruments 
can be followed up when necessary with high resolution equipments. 
On the theoretical side there has likewise been a correspondingly 
rapid progress. ‘The shell model of nuclei has proved of great value in 
the elucidation of the spins and parities of the lower-lying excited states 
of nuclei, the states which feature mainly in the transitions investigated 
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by the methods of y-spectroscopy, at least as applied to radioactive 
substances. The accurate evaluation of theoretical conversion coefficients, 
including those for energies approaching the K-shell binding values, has 
proved very important, as has the deduction of radiation probabilities 
distinguishing between magnetic and electric multipole transitions. This 
work has been supported and extended empirically by the systematic 
study of K conversion coefficients (Rose 1952) and K/L. ratios 
(Goldhaber and Sunyar 1951), especially in the examination of isomeric 
transitions. Internal pair production has been put to use experimentally 
and the theory checked (Rae 1949, Slatis and Siegbahn 1951). 

Closely associated with these advances in theory and experiment is 
the advent of the atomic pile as a high intensity neutron source capable 
of yielding a great variety of strongly radioactive elements. A great 
abundance of material for study has thus been created largely by the 
process of neutron capture. The relatively limited number of substances 
which must be prepared by bombardment with high energy deuterons 
or other particles do not seriously affect this happy situation. It is 
true that the process is not too well suited to the problem of preparation 
of carrier-free sources but even here the development of isotope 
separators, in particular the electromagnetic process, gives a method of 
overcoming this limitation in some cases. This same separation 
technique has proved useful in preparing sources of well-defined thickness. 

Many rather more mundane advances have contributed to the 
improvement and some of these will be considered in the following 
pages. Here we may note the improvements in methods of source 
preparation and support, the total thickness being sometimes several 
peg/cm* at most. 

As regards coincidence techniques, the spectrometers have been 
made sufficiently fast and efficient to allow coincidences between the 
radiations analysed at the same time in two spectrometers to be studied 
satisfactorily even when the resolving power of both analysers is fairly 
high. For the same reason angular correlation studies have proved 
much more successful. In this same connection we should note that 
scintillation spectrometers are capable of directly measuring the life-times 
of excited states down to values less than 10-® sec. This extension of 
the range of measurable life-time of a y-ray state has proved valuable 
in many investigations. 

It is hoped that these remarks indicate sufficiently well the nature 
and scope of the progress that has been and is being made in this field 
of nuclear physics. In the following pages we shall concern ourselves 
mainly with the critical examination of the available experimental 
techniques, attempting to indicate some of the advantages and limitations 
of the various methods. An effort will be made to assess the relative 
performance when this seems desirable and practicable. 

We propose to discuss the nature of the different methods and to 
follow this with consideration of the low resolution instruments. 
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A description of the performance of the higher resolution instruments 
concludes the review. 


§2. PRINCIPLES OF THE METHODS 


Many methods of analysing gamma-radiation are practicable and 
here we include in our list those which seem to have been used fairly 
widely. There are the following : 


2.1. Determination of Absorption Coefficient 
This technique does not lend itself to the examination of complex 
radiations and it is now practically obsolete except in a few special 
cases where the discontinuity in coefficient occurring at the K-absorption 
edge can be exploited successfully to estimate the energy of a radiation 
within rather narrow limits. It applies therefore mainly to x-rays and 
y-rays lying in the energy range 0 to 100 kev. 


2.2. Absorption or Range-analysis of Secondary Electrons 

The most practical form of instrument in this method consists of a pair 
of counters placed side by side with a space between in which absorbing 
foils can be placed. Secondary electrons (photoelectric or Compton) 
released from a sheet of metal or the face of one of the counters nearest 
to the y-source pass through both and produce coincidences. The energy 
of such secondaries can be estimated by examining the variation of the 
coincidence counting rate with the amount of absorber put between the 
detectors. This technique has proved valuable in the analysis of low 
intensity radiations and it can yield rather accurately the energy of the 
hardest quanta present. It fails to separate the components of a 
heterogeneous beam. 


2.3. Photodisintegration 

Nuclear reactions of the (y, p) type can be used since the energy or 
range of the released protons is related to the quantum energy of the 
radiations. The substance employed is usually deuterium and the 
forward-going protons have a kinetic energy close to the quantum 
energy. The cross section, particularly near the threshold, is rather low 
and the evaluation of the proton range not very precise. It has been 
employed fairly widely by loading nuclear emulsion plates with heavy 
hydrogen. Background difficulties tend to be acute but it has proved 
extremely sensitive in certain instances. 

At higher energies the method is of value in determining approximately 
the intensity lying within defined energy bands. A variety of materials 
of different thresholds can be exposed to the radiation and their induced 
radioactivity used to estimate the spectrum of the y-ray source 
(e.g. synchrotron). 
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2.4. Momentum of Secondary Electrons 


The measurement of the curvature of secondary electrons in a magnetic 
field of known strength forms the basis of most of the methods of energy 
analysis. The principle is most clearly illustrated in the cloud chamber 
studies where the radius of curvature p of the tracks in field H can be 
measured directly to give Hp(=mvje). The efficiency in this case is 
rather low since thin radiators are placed at the side or within the 
chamber and the general background of tracks involves the experimenter 
in a process of rather stringent selection to ensure that true photo- 
electrons, Compton recoils or pairs are selected. In addition the 
magnetic field cannot readily be made homogeneous and the curvature 
itself is modified by eddy processes, diffusion, and scattering within the 
gas. The technique is not a high resolution one and its chief advantage 
lies in the directness of the evidence. The analogous method employing 
electron-sensitive nuclear emulsions offers little promise as the estimation 
of curvature must be replaced with the more hazardous estimation of 
range. Moreover available emulsion thicknesses would appear to limit 
the application to radiations of relatively low energy. 

2.5. Total Ionization 

The energy required to produce an ion-pair has been shown to be 
remarkably constant for at least two gases, argon and nitrogen, over a 
large range of energy of the primary electron. It is probably safe to say 
that the same is true for all the rare gases and it is indeed unlikely that 
marked variation will occur with any of the well-known gases and 
vapours. If then the primary electron is released in the photoelectric 
process by a y-quantum and if it is brought to rest within the gas the 
total number of ion-pairs formed will be a good measure of its energy 
and hence of that of the quantum. This is the principle on which the 
proportional tube spectrometer operates and it has proved remarkably 
successful in the analysis of soft y-radiations (and x-rays). With 
photoemission in the gas itself it operates well up to ~100 kev and with the 
use of thin radiators of heavy elements it can be used satisfactorily up 
to ~500 kev. It is possible that further studies will extend the range 
as high as several Mey. 

Another ionization method is possible where a liquid ionization 
chamber takes the place of the proportional counter. Here the absence 
of multiplication of the electrons of the ion-pairs initially created makes 
the position difficult at low energies (<100 kev) but it should operate 
relatively well at higher energies when suitable thin radiators are 
employed. ‘The liquid that appears most suitable, argon, seems to have 
most of the necessary qualities but the writer is not aware of anv 
experiments performed along these lines. ; ° 


2.6. Total Light Emission 


There is clear evidence that in many fluorescent substances the 
remarks made in §2.5 above regarding ion-pairs apply equally well to 
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the yield of fluorescent quanta. If this total yield is measured successfully 
we have a measure of the energy of the quantum giving rise to the 
primary photoelectron. The measuring device used in the scintillation 
spectrometer is the photomultiplier and, since the fluorescent substances 
commonly employed are solids or liquids within which the primary 
electron (or electrons) are released, the efficiency can be maintained at 
a very high value over a big range of quantum energy, and the range 
of the primary (or primaries) confined within the volume of the 
scintillator. Partly for this reason and partly because fairly heavy 
elements are present in good scintillators (e.g. I in NalI(Tl)) the pair 
process as well as the photoeffect can be utilized. 


2.7. Wavelength 


Selective reflection of y-radiation from the atomic lattice of a crystal 
can be used in y-spectroscopy in the same way as in x-ray spectroscopy. 
This very fundamental method of determining quantum energy 
(hv=hc/A) has been employed from time to time since its early application 
by Rutherford and Andrade (1914 a, b). Two serious difficulties are 
encountered—the need for strong sources of high specific intensity and 
the accurate measurement of the small angles involved. The first of 
these has been alleviated to some extent by exploiting the bent crystal 
and by the introduction of the scintillation spectrometer as the detector. 

Clearly y-spectroscopy is a subject with wide ramifications and in the 
following pages it is proposed to discuss two of the relatively new 
approximate (low resolution) methods and two high definition instruments. 


§3. THE PROPORTIONAL TuBE SPECTROMETER 


Early work with the proportional tube as a spectrometer was concerned 
as much with soft B-emitters as with y-radiation but here we are concerned 
with the latter. Curran ef al. (1948, 1949) showed that well-defined 
peaks appeared on the pulse distribution observed with the counter 
when it was exposed to homogeneous x-rays or soft y-rays, e.g. the 
fluorescence x-radiation of copper and the y-radiation of RaD. Since 
then the technique has proved a powerful and flexible one and the 
energy range has been extended to beyond 500 kev. Some of the factors 
affecting closely its application as a spectrometer will now be considered 
in more detail. 


3.1. Linearity 


It was soon shown, for the gases, argon and nitrogen, that the energy 
dissipated by an electron in forming an ion-pair was constant to a close 
approximation over a wide range of energy of the primary. Initially 
this was established for the range 2 to 50 kev but recent work shows 
that it holds from about 0-2 to 150 kev in nitrogen, and, to a somewhat 
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less close limit to about 600 kev in argon. The constancy is good to about 


2° over most of the ranges given here but at about 0-2 kev the 
measurements are less certain than this. This linearity is at least a most 
convenient feature as the calibrating radiation need not be close in 
energy to those investigated. Useful calibrators may take the form 
of radioactive sources emitting K-capture radiation or sources in which 
strong internal conversion gives rise to K and L rays or sources emitting 
soft y-radiation, e.g. RaD y-rays of energy 46-5kev. Many suitable 
sources have been proposed but lacking these the fluorescence x-rays 
of different elements can be excited by any suitable means and 
transmitted into the spectrometer. 


3.2. The Pulse Profile 

Franzen (1951) has given a lucid description of the nature of the 
signals from both pulse ionization chambers and proportional counters. 
In the latter the voltage pulse appearing at the collecting electrode (the 
positive wire) is caused largely by the motion to the cathode of the 
positive ions formed in the avalanche, and to a much less important 
extent by the collection of the electrons. The collection time 7' of the 
positive ions is given by 


T'~pb? log (b/a)|(2uV,), . . . . . ~ « (I) 


where p=gas pressure, b=radius of outer cylinder, a=radius of inner, 
=mobility of the ions and V),=applied voltage. For typical counters 
T lies in the range 0-1 to 1010-8 sec. This suggests that the pulse 
is of long duration but if the motion of the positive ions is considered as a 
function of the time t from the instant of their formation it is found that 
the signal S(t) at the collector is 
é 
S(t)=— 2C log (b]a) log{1+[(b/a)?—1]t/T}, = rar Ghote vie (2) 
where e=the electrcnic charge and C=the capacity of the chamber. 
The expression for S(t) shows that the signal rises very rapidly as 
¢ increases from zero and at t=aT'/b it has the value —e/(2C) to a close 
approximation, that is, one half of the final value —e/C. However for 
practical counters a/b is ~10-3 so that in times between 0-1 and 
1-0 10° sec the signal reaches half value. In recent forms of the 
spectrometer relatively fast amplifiers have been used in conjunction 
with the counters and the late slower development of the signal to full 
amplitude is rejected. Such rejection does not appreciably affect the 
proportionality of the pulses and very considerable advantages accrue 
from the speedy response of the instrument. The presence of the 
localized ion cloud within the tube after the signal has decayed must be 
regarded as a potential source of error at high counting rates 
(> 10% per sec say). 
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3.3. Statistical Fluctuations 

The fluctuations in amplitude of pulses arise in two ways : 

(1) The number of ion-pairs released along the total track of the 
primary particle fluctuates about an average value. 

(2) The multiplication of the electrons in the vicinity of the wire varies. 

Fano advanced arguments to support the view that when an ionizing 
agent expends the whole of its energy within a volume of gas the 
fluctuations about the mean value of the number of ion-pairs formed, 
v say, should be less than that given by Poisson distribution 


ME eer VeCR Te ved, oe ae Ssh Srl acs 01343) 
where P(n) is the probability that n pairs are formed. Following Fano 
(1947, 1946), Franzen expresses this argument in the form of the equation 


(J’—Jo)?=(n—«/W)?J,/n, Payee neal Leer (4) 
where J’ is the actual number of ion-pairs formed, J, is the average 
number, W the energy dissipated for each pair produced, « the energy 
lost per collision and n is the number of ion-pairs liberated per collision. 
For a Poisson case (J’—J,)?=J, and the factor 
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by which the expressions differ is generally less than unity. By his 
definition of « Fano takes into account excitation and ionization in 
which the liberated electron may or may not be sufficiently energetic 
to produce ionization in the surrounding gas. As an example he deduces 


that for atomic hydrogen bombarded with electrons 

ge Ace hale ae) 4) oe ce Pe At at aby oe 2k; (0) 
As yet it appears impossible to check these conclusions except with the 
proportional tube as an analyser of rather soft homogeneous radiations. 
In the pulse ionization chamber detecting more energetic particles the 
statistical fluctuations are considerably less than the instrumental 
contributions to line width, e.g. amplifier noise. On the other hand the 
lines observed in the proportional counter, assuming otherwise near 
perfect conditions, have widths determined by the variations of the 
gas multiplication factor as well as by the fluctuations in the number of 
primary ion-pairs. Curran, Cockroft and Angus (1949) succeeded in 
separating these two contributions to the width. They examined the 
spread in the pulse amplitudes produced by single slow photoelectrons 
released by light from the cathode (see fig. 1) and from this spread. 
deduced the spread due to the gas multiplication process alone when a 
finite number of primary electrons were multiplied near the wire. By 
comparing the width of a semi-empirical line of this kind with the 
observed larger width found in practice when primary electrons absorbed 
in the gas released on the average the same finite number of electrons 
(as ion-pairs) they separated out the spread introduced by fluctuations 
in the ion-pair formation. The results are shown for a typical case in 


Pulses per energy interval. 
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fig. 1. It appears certain from the results that (J’—J,)? is less than Jo 
and the value of y is between 2/3 and 1 for a 50: 50 mixture of argon and 
methane. The two contributions appeared approximately equal and 
the total statistical fluctuation could be expressed” in the form | 
(2n/3-+2n/3)¥/2/n. The analysis is rather difficult experimentally and 
further work with more gases would be useful. The results of Hanna et al. 
(1949) are in reasonable agreement with the above. 

Franzen has pointed out that the distribution of amplitude of fig. 1(a) 
has a maximum while the distribution which he deduces on a simple 
duplication theory should decrease monotonically. For the probability 
about a mean number WN he finds, for V>1 
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(a) Pulse size distribution for single electrons in proportional counter. 


(6) Experimental analysis of Cu K, rays. Dashed curve shows component due 
to variations of gas gain alone. 


It is possible that his assumptions are over-simple but it would be useful 
if the effect of the bandwidth of the amplifier on the pulse amplitude- 
distribution were examined in some detail. The relatively fast response 
of typical amplifiers could conceivably be the source of some of this 
disagreement with the theory. Certainly the observed value of y might 
be higher in this case than that obtained with an amplifier designed to 
reproduce faithfully the whole of the pulse profile even though such an 
amplifier is not desirable in a practical system. 
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It is clear that the proportional counter is a rather low resolution 
instrument but in practice this proves very little disadvantage in its 
special province. For example, it can generally distinguish between the 
K or L x-radiations of neighbouring elements and on the whole such 
resolution proves adequate for most purposes. 


3.4. Voltage Stabilization 

Closely linked with the question of resolving power and accuracy of 
energy measurement of the tube is that of voltage stabilization. With the 
gas mixture frequently employed, argon plus ~10° of methane, and 
with nitrogen, the gas gain increases exponentially over a wide range with 
increasing voltage. Typical tubes show a two-fold rise in gain for an 
increase of between 100 and 140 volts. This suggests that for very good 
stability of pulse amplitude (variation <1°% say) the voltage should be 
held constant to about 1 part in 10‘ if possible. An alternative method 
of maintaining constant conditions is to use a source of homogeneous 
radiation to provide pulses which are monitored and used to control the 
voltage to give constant amplitude (Wilkinson 1950). 


3.5. Amplitude Analysis 

Here again we have an important experimental limitation to the 
resolving power. The possibility of slow changes in the counter perform- 
ance (among these voltage drift) makes it necessary, for optimum 
results, to obtain the maximum amount of information on the pulse 
spectrum in as short a time as possible. Multi-channel electronic kick- 
sorters are generally sufficiently fast to handle the output pulses but 
100-channel instruments (for the range 0 to 100 kev say) are apt to be 
bulky and expensive. Recently Hutchinson and Scarrott (1951) have 
introduced an ultrasonic delay line method of amplitude analysis operating 
on an amplitude-to-time conversion principle and such an instrument, 
using mercury column or nickel wire as the storage house for the 
information, seems to possess many of the properties ideally suitable for 
the proportional tube. The display of the information on a cathode-ray 
tube is not the least of the attractive features of this system. In actual 
fact a rather satisfactory method of analysis is secured by the straight- 
forward display of all the pulses on the cathode-ray tube and the recording 
of these on moving film when required for detailed study. The energy 
intervals can be set at any desired value and the only drawbacks are the 
laborious nature of the film reading and the tendency for personal 
factors to appear in such reading. 


3.6. General Remarks 
Variations of the electric field strength and shape in the regions near the 
ends of the wire are to be avoided if possible since they lead to non- 
uniformity of gas gain and uncertainty of definition of the counting 
volume. The end-correction technique introduced by Cockroft and 
Curran (1951) eliminates these disadvantageous features and makes the 
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whole of the counter uniform in its response. A typical counter assembly 
with field-correcting tubes is shown in fig. 2, in schematic form. 

The excitation of fluorescence x-rays at the cathode, at the windows or 
at the source support is to be avoided. Polished aluminium makes a 
suitable lining or cover for the surfaces involved. It appears that cast 
aluminium tubes in general do not behave well as the main cylinder of the 
counter. Frequently the photoelectric absorption of the radiation occurs 
almost entirely in the gas and there may be an appreciable possibility of 
‘escape ’ of the K-radiation of the gas depending on the pressure and the 
size of the vessel. For argon this radiation has an energy of about 
2-6 kev and it is rather seldom that this soft radiation will escape 
from the gas without further interaction. When the effect is present it 
leads to the appearance of an escape peak at an energy lower than that of 
the quanta under investigation (by ~2-6 kev). The effect undoubtedly 


Fig. 2 


Schematic representation of proportional counter with different source positions 
indicated at a, b and c. T=t.eld-adjusting tube, G=grounded guard 
tube. 


contributes to the breadth of the observed lines even if a separate peak is 
not noted. The phenomenon appears more markedly in the heavier rare 
gases such as krypton and xenon (West and Rothwell 1950). 


3.7. Arrangement of Sources 


The spectrometer is fairly efficient in the range up to ~100 kev when 
high pressure argon (7 atmospheres is readily used) is the filling gas but 
better efficiencies are obtainable with krypton or xenon. These remarks 
apply when photoelectric interaction with the gas is involved. When 
internal conversion of the y-radiation occurs the photoelectrons emitted 
from a thin source (or a gaseous source) may be measured and the quantum 
energies deduced. Generally it is satisfactory to mount the radioactive 
material as a thin layer on a suitable surface within the counter. A number 
of arrangements have proved practicable :— 

(1) With very low intensity sources the material may be spread over 
the whole of the inner wall (generally over polished aluminium sheet). 
This procedure allows one to use a large mass of material in the form of a 
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relatively thin layer, e.g. 1 g at a thickness of 1 mg/em?. More is said 
later about the effect of the support on the resolution. 

(2) With sources of higher specific intensity the source may be mounted 
on a restricted area of the cathode surface or even better on a very light 
foil of say nylon as indicated in fig. 2 at a. The source can be made as 


thin as desired say 10-® g/cm?. 


Intensity/energy interval. 


(3) The source can be mounted on a probe as indicated at b in fig. 2. 
The probe may take the form of a small loop on the end of a stout wire, the 
loop supporting a nylon foil on which the thin source is deposited. This 
arrangement gives a 47 geometry but if the absorption in nylon is undesir- 
able the supporting material may be thick (27 geometry). 


Fig. 3 
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Electron lines of 1°Te*. The lower curve in each is for thin nylon support, no 
backing. The effects of thick aluminium and silver backings are shown 
by the upper curves in (a) and (6) respectively (from unpublished results 
of J. G. Balfour at Glasgow). 


(4) The source may be supported as in (3) above but maintained 
(fig. 2 at c) on the plane passing through the end of the wire and perpen- 
dicular to it. 

In (3) and (4) the supporting wire or thin rod is maintained at the 
required potential to give no distortion of the electric field. In a number 


_ of experiments of the author and co-workers no effect of charging of the 


source has been observed. The difference here as compared with the 
magnetic spectrograph is probably accounted for by the presence of gas in 
the counter and the ionization normally occurring within it to a greater or 
lesser degree. 
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As an illustration of the technique we show in fig. 3 the spectrum of 
internal conversion electrons released in the decay of the isomer }*°Te*. 
The source was a few times 10-® g/cm? in thickness and mounted on thin 
nylon foil at an aperture (arrangement (2) above). Thick sheets of alum- 
inium silver or gold could be brought behind the source and close to the 
nylon. The photoconversion peaks are in all cases clear but the scattering 
effects of the thick and heavy backings are fairly obvious. At the same 
time it is apparent that at the comparatively low energy of about 50 kev 
the distortion introduced by a thick aluminium support is not serious and 
it is easily possible to estimate it. This result for electronic radiation is 
reassuring. It was found useful to apply a magnetic field at right angles 
to the plane of the source and the wire of the counter, especially when the 
source emitted electrons of relatively high energy (~100 kev). Clearly 
the proportional tube is an extremely efficient spectrograph for photo- 
electrons produced by internal (or external) conversion. In many cases a 
rapid and accurate assessment of conversion coefficients in K and L shells 
can be made. When a source emits f-particles and de-excitation of the 
daughter element by internal conversion follows the spectrum of pulses 
consists of several distinct groups, e.g. B-rays, B-rays-+K—conversion 
electrons, B-rays+L conversion electrons and so on. The areas under 
various portions of the spectral curve can be used to estimate conversion 
coefficients. This procedure was adopted by Wilson and Curran (1951) in 
analysing the radiations of ?%Hg. Similar procedures are possible with 
the scintillation spectrometer. 


3.8. Analysis of Harder Radiations 


The application of magnetic fields allows us to extend the useful range 
of operation of a high pressure proportional tube. A variety of field and 
counter geometries has been explored. A fringing field was used to restrict 
the path of the primary particles in two dimensions and indeed restriction 
in three dimensions is feasible. The method generally employed at 
Glasgow is indicated in fig. 4 inset. The magnetic field is parallel to the 
wire and the source is mounted at the centre or end of the counter, as close 
to the wire as is reasonably practicable. With a counter wire 8 in. in 
length the end-position of the source seems to give satisfactory results, at 
least when the gas pressure is 7 atmospheres. When internal conversion is 
effectively absent external conversion in a thin radiator of high atomic 
number can be employed. In fig. 4 we show the peak due to external 
conversion of the y-radiation of "Cs in a thin gold radiator (16 mg/em?). 
The radiator covered the end face only of a small brass container con- 
taining a very weak source of the isotope ®’Cs. With greater area of 
radiator and substitution of uranium metal for gold the performance 
would be more impressive. The example is due to A. L. Cockroft and it 
demonstrates that it is possible to use the proportional tube as a sensitive 
spectrometer for y-quanta up to energies exceeding 500 key. It is 
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interesting to note that the main part of the spectrum, produced by 
Compton electrons emitted into a solid angle of 27, has a flat-topped shape 
rather similar to that obtained by Svartholm and Siegbahn (1946) in 
their studies of y-radiation with the magnetic spectrograph using external 
conversion and a relatively small solid angle. 

Extremely low intensity sources could in theory be analysed with the 
help of external photoeffect by spreading the source over the cathode 
of the counter and covering with a light absorber (to remove f-rays) 
coated with a suitable thin converter. This technique has not been 
employed largely on account of the introduction of the scintillation 
spectrometer. 


Fig. 4 
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Pulse Size (arbitrary scale) 
(a) Compton electrons from source container. 


(6) Compton electrons plus photopeak due to gold converter excited by gamma- 
rays of 187Cs (energy 662 Kev). 


The resolving power obtained at energies in the region of 500 kev does 
not appear to be in agreement with predictions based on observations at 
low energies. A primary electron of energy 500 kev releases about 
n=1-7X< 104 ion-pairs and we would expect a fluctuation of (1-67n)1/2 
i.e. 1:7 102. A line-width of-about 1° might be anticipated. The 
results fall far short of this quality and it is still not clear where the difficulty 
lies. More work in this direction is required. 


424 S. C. Curran on 


3.9. Coincidences 


The use of proportional tube spectrometers in coincidence with each 
other is a very suitable method of investigating the softer radiations 
from a variety of sources. J. G. Balfour has examined the radiations 
from RaD (mounted as a very thin source on a thin nylon foil) by means 
of two proportional tubes placed side by side within a common container. 
The source was placed at an aperture on the cathodes on the line of 
contact. He was able to record on a two-dimensional display system the 
coincidences of B-rays with electrons, B-rays and electrons with gamma- 
or X-rays and gamma- or X-rays with gamma- or x-rays (suitable absorbers 
over the source made these different combinations possible). All the 
known energy range for both particles and quanta was covered. The 
resolving time was relatively poor (10sec) and the observations were 
consistent with a simple decay scheme for this source, that is, a S-particle 
of energy about 16 to 18 kev followed by de-excitation from an excited 
state at 46-5 kev. Improved resolving time for the coincidence circuit, 
readily possible, should help to throw light on any further complexities 
in the decay scheme. The work is described here mainly to serve as an 
example of the possible application of the method. More generally, 
the proportional tube can be used in conjunction with the scintillation or 
magnetic spectrometers as the analyser of the softer radiations and so 
prove very valuable in elucidating decay schemes. 

An arrangement possibly useful for coincidence studies (and which 
has not yet been tested) may be suggested here. If the central wire of a 
counting vessel is in the form of two portions held together near the 
centre by a small insulator (glass bead for example) we have a pair of 
counters and at the common boundary there is not much ‘ fringing’ of 
the electric lines of force. This deduction is based on the results of Curran 
and Rae (1947) with beadel Geiger tubes. The ‘ fringes’ extend over 
the cathode for at the most a distance equal to several bead diameters. 
If a source is mounted on or within such a small insulator we have a 
useful arrangement for coincidence studies. Equally good performance 
is of course possible in this split counter if the source is held near the 
axis on the boundary plane. 


3.10. Photosensitive Proportional Tube 


A proportional tube with a semi-transparent cathode of sensitivity 
equal to that of a typical modern photocell or photomultiplier would 
serve many useful purposes. The difficulty in producing such a device 
is presumably due to the fact that only the rare gases or mixtures of the 
same would not react with surfaces of the antimony-caesium type. 
In these pure rare gases we are troubled with metastable atoms and 
instability of gain seems to be manifest at a multiplication factor of 
10° or less. If more suitable gases such as nitrogen could be used (or 
mixtures with argon) and a stable gas gain of 10° or 108 obtained a 
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photomultiplier of great simplicity would result. Ay an example of the 
application of such a device to y-spectroscopy we could use it as an efficient 
detector of the light emitted from a large volume of liquid used as a 
scintillator. Such a large volume of liquid could absorb high energy 
y-quanta and the system would serve as a simple and efficient spectro- 
meter in the high energy region. An arrangement incorporating such 
a large volume of liquid has been described but the light collection 
required the operation of sixteen photomultipliers in parallel. 


§4. THe ScrnTinLaTion Mrraop 


The second low resolution technique to be discussed is that employing 
scintillation counters. In contrast to the proportional tube method we 
have to consider all three important mechanisms of interaction of 
gamma-rays and matter—photoelectric, Compton. and pair production. 
This extension has given rise to three rather different types of instrument 
each of which was designed to operate primarily by virtue of one of the 
three mechanisms. This does not include many other types of spectro- 
meter based on the same principles but the three will suffice to illustrate 
the main trends in the present development. To some extent the three 
types can be regarded as applying to different parts of the energy range— 
the single crystal method to the range up to about 2 Mev, the Compton 
coincidence or Hofstadter type to the range from about 1 to 3 Mev say and 
the Johansson or triple coincidence type from about 2 Mev upwards. 
As we shall see the types may be applied over the whole of the energy 
range (except that there is a threshold for the third at 1 Mev) if necessary. 


4.1. The Single Crystal Spectrometer 

The arrangement in this case is sketched in fig. 5 (a), and is seen to 
consist of scintillating crystal-+photomultiplier. In the range up to 
~1mev the photoelectric interaction predominates at least when the 
scintillator is a crystal of high density such as sodium iodide, thallium 
activated, NalI(Tl). The conversion efficiency for particle energy to 
light energy is rather low (about 10%) even for relatively good 
scintillators such as NalI(T1) and with a normal light-collecting arrange- 
ment round the crystal it is found that between one and two photoelectrons 
are released from the cathode of the photomultiplier per kilo-electron-volt 
of energy dissipated within the crystal. At quantum energies in the 
region hv=100 kev we have at most about 200 electrons accelerated to 
the first dynode of the photomultiplier and a statistical fluctuation in 
pulse amplitude of +7% arising from this limitation alone. The 
instrument is therefore low in resolving power and inferior to the 
proportional counter in this region. As is to be expected the resolution 
improves with increasing energy (approximately as v/2) but, as in the 
proportional tube case, the improvement is not maintained to energies 
of several mev. The reasons for this failure will be discussed below. 
A set of sample spectrograms is given in fig. 6. They show the spectra 


P.M. SUPPL.—OCTOBER 1953 26 


426 S. C. Curran on 


of the gamma-rays of RaD, energy 46-5 kev, of ®°Co, energy 1:15 and 
1-33 Mev, and of 24Na, energy 1:38 and 2-75 Mev. They are taken with 
similar instruments under somewhat different conditions and by different 
observers and hence they are not strictly comparable. Nevertheless 
they show clearly the possibilities and limitations of this simple 
spectroscope. Clearly it is a very powerful tool for the rapid analysis 
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(a) Single crystal spectrometer. (b) Double crystal. (c) Triple crystal. 
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of y-spectra at low resolving power and it has played and will continue 
to play a major part in the elucidation of decay schemes. It can operate 
very successfully in coincidence arrangements, for example in separating 
a particular y-transition from a number of such transitions and so 
permitting the analysis of an associated radiation, say of the B-type in 
; second spectrometer. Obviously it is an extremely suitable instrument 
or investigating y~y coincidences and its size makes it well adapted to 


such tasks as observing the angular distribution of y-radiations of 
particular energies. 
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4.2. The Double Crystal Spectrometer (Hofstadter and McIntyre 1950) 


The simplest arrangement is shown in fig. 5 (b). The gamma-rays are 
collimated and allowed to fall on a crystal. The Compton scattering of 
the quanta in this crystal can be examined by virtue of the forward-going 
recoil electrons and the scattered quanta which leave this crystal can be 
detected in a second crystal placed to receive those scattered backwards 
(say between 120 and 150°). The recording of coincidences serves to 
separate the Compton process from other types of interaction and if the 
pulse amplitudes observed in the first simple spectrometer are examined, 
the signals being ‘ gated’ by the coincidences, we have a rather clean 
spectrum for the forward-going recoils. A single peak appears on this 
spectrum for each energy value of the y-ray components of the beam and no 
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Gamma ray spectrum of 14°La (Bannerman, Lewis and Curran 1951). 


confusion arises in the interpretation of the results. The backward scat- 
tered quanta are reasonably homogeneous and in practice the resolving 
power of the instrument depends mainly on the statistical fluctuations in 
the pulses occurring in the first crystal, assuming of course that the 
radiation is fairly well collimated. The pulses in the second crystal 
serve merely to select the process and their amplitudes need only exceed a 
minimum value. The detecting scintillators are usually organic crystals 
such as anthracene or liquids such as terphenyl in xylene. The colli- 
“mation reduces the overall efficiency very considerably but no serious 
effort to raise it to high values has been reported. Clearly the geometries 
are such as to make a large increase possible, and it might be as much as a 
few per cent. An example of the analysis with this instrument is given in 
five 


4.3. The Johansson Spectrometer 


The arrangement of three crystals used by Johansson (1950 b, 1952) is 
shown in fig. 5 (c). Gamma-rays from the source fall on a scintillator of 
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relatively dense material, e.g. NaI (Tl) and at energies hy>2m,c? where 
my=rest mass of an electron and c=velocity of light, the pair production 
gives rise to a negative and positive electron each of which dissipates its 
kinetic energy within this scintillator. At the instant of annihilation of 
the positron two oppositely directed quanta of energy hyv=0-51 Mev are 
emitted and these may be detected in a pair of scintillators suitably 
placed near the first. Triple coincidences arising in this fashion are used 
to ‘gate’ the pulses observed at the output of the first photomultiplier 
and so we separate the pair process from the other interactions. The peaks 
observed in the spectrum of these output pulses have an energy correspond- 
ing to the total energy of the electron-pair and this equals the y-ray 
energy less 27,¢?. 
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(a) Gamma spectrum for boron+ proton reaction observed with three-crystal 
system. 

(6) Gamma spectrum for fluorine real observed with single large erystal 
Both by pair process in NalI(T). i Pa a 
Some use of this type of spectrometer has been made at Glasgow by 
J . G, Rutherglen and P. Grant. At the same time they used the 
single crystal method (§4.1. above) on the same class of y-radiations— 
those emitted in the capture process for protons in light nuclei. Two 
examples of this work are shown in fig. 8 and two conclusions can be made 
on inspection. At quantum energies ~5 Mev we expect ~10‘ photo- 
electrons should leave the photocathode. Thus we might Ae ity 
line width if the major part of the width was contributed by senieeeal 
fluctuations in thisnumber. The line widths in both examples very much 
exceed this value. It appears probable that there are several contributor 
effects :—(a) non-uniformity of conversion efficiency of the scintillater 
perhaps associated with difficulties in obtaining equal activation by 
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thallium throughout the volume, which exceeded a cubic inch, (6) 
difficulty in collecting the fluorescent light equally well from all regions 
of the scintillator for directing onto the photocathode, and (c) non- 
uniformity of photosensitivity over the cathode area. Of these perhaps 
the second is the most intractible at present and further work on all should 
lead to considerable improvement in resolving power of the scintillation 
spectrometers, at least in this high energy region. 

In the operation of the proportional counter certain factors which had 
direct bearing on its accuracy, resolving power and so on were considered. 
Similar factors are involved with the scintillation spectrometers and we 
have discussed resolution briefly above. We now list a few of the more 
important points. 

Innearity. Here we are concerned with the constancy of the energy 
required to produce a useful photon of fluorescent radiation. This problem 
has been investigated in some detail recently and it appears that most of 
the widely used crystals, sodium iodide, anthracene and stilbene for 
example give a good degree of satisfaction for electrons as the primary 
exciting radiation. The best linearity is observed with sodium iodide 
and indeed it is possible that the output pulse amplitude is directly 
proportional to the energy dissipated, over the range from say 2 kev to 
3 Mev. The departures recorded by some observers are slight and they 
could be due to experimental difficulties of various kinds. On the other 
hand all the more detailed studies of the organic crystals, in particular 
those with anthracene and stilbene, seem to agree in showing that the 
linearity is indeed an approximation to a slow variation. The energy 
dissipation per photon decreases fairly rapidly as the energy of the 
primary electron increases from a few kev to about 20 kev and then more 
slowly. Above about 100 kev the response curve is nearly linear. Very 
similar behaviour has been noted with the liquid scintillators such as 
terphenyl in toluene, terphenyl in xylene and others. For the heavily 
ionizing particles the situation is much less satisfactory and indeed the 
non-linearity is a major impediment in such work as fast neutron spectro- 
scopy. The organic materials show a very marked variation in the 
amount of energy dissipated per photon over the whole range of energy 
values involved in most experiments and the performance of the inorganic 
materials such as Nal(TI) calls for careful examination. None of the 
better-established methods of y-spectrometry involves the excitation of 
crystals or liquids by protons or alpha-particles so we do not consider the 
matter in any detail. The subject has previously been reviewed by the 
writer (Curran 1953). 

Stability of performance. The voltage applied to the photomultiplier 
of a spectrometer determines the multiplication factor and hence must be 
stabilized rather well. For best results 1 part in 10¢ is a good target but 
on the whole satisfactory work can be done with the supply held constant 
to 1 in 10%. Fatigue effects may be exhibited in some tubes and it is wise 
to check their performance from time to time. It is likely that most of 
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such a long-term change occurs at the photocathode. Another source of 
variation in overall response is found in slow chemical changes at the 
surface of and within the scintillator. Such changes are accelerated by 
bombardment with radiations, particularly protons and alpha-rays but 
they do occur, although to a much smaller extent, with electrons (Black 
1953). The scintillator has in fact a very long but none the less finite 
life. More rapid deterioration can frequently be traced to direct chemical 
change such as occurs when sodium iodide, which is hygroscopic, is 
exposed to a humid atmosphere, or when certain materials which act as 
quenching agents go into solution. Most of the effects noted here are 
fairly readily eliminated and the stability of the instrument can be made 
really excellent. 

Amplification. The multiplier is effectively an amplifier of very 
great bandwidth and care must be exercised if post-multiplier amplifica- 
tion is required. The E.M.I.-type photomultipliers are obtainable with 
14 dynodes and they may be used in most experiments without further 
amplification. An R.C.A. tube (type 4646) in the developmental stage 
contains 16 dynodes and a saturating signal can be produced with a single 
photoelectron from the cathode. It has the additional advantage of 
feeding the output pulse directly into a low-impedance line (125 ohms). 
With such high gain tubes saturation is a potential source of error in 
spectroscopic work and it may be necessary to take some of the signals 
from one of the dynodes rather than from the collecting anode. The 
tube has been described by Greenblatt et al. (1952). 

Intensity. There does not appear to have been any comprehensive 
experiments on the problem of assessing the relative intensities of 
y-transitions with the various types of scintillation spectrometers. The 
single crystal method operating by photoelectric absorption is the most 
direct but the discontinuities at absorption edges leads to uncertainty 
of the cross section and the nature of its variation at low energies 
(<200 kev say). At somewhat higher energies the cross section changes 
smoothly and more reliable estimates can be made. Another difficulty 
arises from the presence of Compton background on the spectra. These 
questions have been studied very recently by Prescott (1953). He finds 
a useful criterion for the performance of the spectrometer is the value 
of the peak to trough ratio (height of line due to photoelectrons to depth 
of trough on the low energy side of the peak). His semi-empirical 
treatment leads to consistent results and shows that fairly accurate 
intensity measurements are possible with the system. 

The cross section for Compton interaction is known to be given 
accurately by the Klein—Nishina formula and clearly the double crystal 
spectrometer can be used to assess relative intensities accurately. The 
scattered quanta may indeed be detected with nearly 100% efficiency 
so that the absorption of the radiations in the measuring crystal is the 
main variable. The spectrometer is certainly well suited to the intensity 
evaluation. The situation regarding the three crystal (J ohansson) 
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spectrometer is rather similar but the more complicated geometry of the 
assembly leads to some uncertainties. Thus the probability of detection 
in the monitoring crystals (for annihilation quanta) depends on the 
point of origin within the central crystal and the average position of 
this origin depends in rather complicated fashion on the energy of the 
incident quanta, the size of the central crystal and other parameters. 
It would appear necessary to deal with the problem in an empirical 
manner. 

High energy quanta. The discussion above has referred mainly to the 
quantum energy range up to about 10 Mev. A bold attempt to use the 
scintillation method at still higher energies has been made by Cleland 
and Koch (1952). They examined the absorption of gamma-rays (of 
energy up to 35 Mev) in a large tank of terphenyl in xylene. They 
found a variation of about +10% in the efficiency of light collection 
from different parts of the large volume, four feet in length, and this 
variation was the limiting factor in the energy resolution. Sixteen 
photomultipliers were used in conjunction with the tank. The use of 
large photocathodes as in the recent R.C.A. Type H-5037 photo- 
multiplier which has a 4-in. cathode would make this instrument much 
more practical. Nevertheless the volume-absorption of the fluorescent 
light seems likely to remain a real difficulty. At best such instruments 
operating on the principle of complete absorption of the energy of the 
quanta, here by multiple acts of Compton scattering, will prove useful 
in low resolution spectrometry of the very high energy y-radiations. 


§5. Tot MaGnetic SPECTROMETER 


Case for high resolution. Hedgran (1951) has given a useful summary 
of the advantages afforded by the study of y-radiations at high resolution. 
Some of these are as follows : 

1. Accurate knowledge of the energies of levels may yield new 
information on the structure of nuclei. 

2. Accurate determinations make the construction of level systems 
much more meaningful. Deutsch (1952) has pointed out that an extra 
significant figure may lead to the confident rejection of a large number 
of possible schemes, as in the case of 1511 for example. 

3. Increased resolving power may reveal the presence of new levels 
or give information on natural line broadening, e.g. the Doppler effect in 
positron annihilation. 

4, Fundamental data for nuclear physics such as the binding energies 
of neutrons in nuclei may be found accurately. 

The need for high resolution is clearly well established and the two 
instruments satisfying this need are the magnetic spectrometer (in various 
forms) and the bent-crystal spectrograph. A number of low resolution 
magnetic spectrometers giving high luminosity have been described but 
it is perhaps fair to say that most of these, so far as y-radiations are 
involved, can be advantageously eliminated in favour of the scintillation 
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and proportional counter devices. We propose to discuss here the 
medium and high resolution magnetic instruments. 

Persico and Geoffrion (1950) have given an excellent account of the 
theory of B-spectrometers of the magnetic type and they consider the 
parameters controlling luminosity and resolving power. On the 
experimental side there has been much activity since 1946 when the 
principle of double-focusing was introduced by Svartholm and Siegbahn 
(1946). New types of instruments have been designed in many 
laboratories but till relatively recently they were mainly of rather 
limited accuracy and very moderate in resolving power. Mitchell (1950) 
has reviewed some of the conventional methods and he shows that the 
energy determinations are generally accurate to between 0-5 and 2% and 
the resolving power is usually ~5%. Some of the instruments depended 
on photoelectric interaction, internal and external and some used 
Compton recoils. Representative of the Compton method we have the 
measurements of Latyshev (1947). The recoil electrons have an energy 
E, given by 

2u(hv)? cos? 6 


o~ 2+ Qu (hv) + (hv)? sin? 6 (8) 


where p=m,c?, the rest mass of the electron and @ is the angle between 
the primary quanta (energy=/v) and the electron. The variation of 
E, with 6 led Latyshev to collimate the y-rays which were incident on 
the thin radiator, so that he confined 6 to small values. In this way he 
secured ‘Compton lines’ rather than a broad rather flat distribution in 
energy but even so the resolving power was rather low. Further work 
on the same lines has been reported by Martin and Parry (1952) and by 
Mladjenovic and Hedgran (1952). It should be noted however that the 
Compton interaction is very suitable in estimating the intensities of 
radiations. This follows from the well-established agreement between 
theory and experiment in the energy dependence of the cross section for 
Compton scattering (Klein—Nishina formula). 

Hedgran had discussed the reasons for preferring photoconversion 
(internal and external) in accurate y-spectroscopy, at least over the range 
extending from several kev to a few Mey. The quantum of y-radiation 
gives up the whole of its energy to a bound electron and its energy hv 
can be found with the accuracy that Z,,, the kinetic energy of the photo- 
electron, can be assessed since the binding energies are known precisely. 
A series of lines is observed, K singlet and L and M multiplets. The 
energy difference of the shells can be used in finding for instance the 
nucleus to which the y-radiation belongs or in checking the accuracy of 
magnetic field measurements but on the other hand the interpretation 
of the results may sometimes be rendered more difficult. 

Three requirements, which were satisfied in the early work of Ellis 
(1932), greatly facilitate the accurate study of y-radiations, viz. : 


(a) strong activities were available ; 
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(6) internal conversion was high ; 

(c) carrier-free sources could be prepared making it possible to have 
them small. 

These unusually favourable circumstances allowed determinations of 
energy to be made to within about 0-1% error. The majority of the 
artificially radioactive elements do not satisfy these three criteria. 
Although they can frequently be strong they show, in general, low values 
of internal conversion coefficient and the most usual method of prepara- 
tion, by (n, y) reaction, does not lend itself readily to the production of 
carrier-free radioelements. The use of (d,p), (d,n), (n, p), and other 
reactions readily induced by cyclotron bombardment, helps in respect of 
(c) above and may indeed serve the further purpose of identifying the 
particular isotope emitting the radiations. The combination of these 
reactions with electromagnetic separation gives us a carrier-free source 
isolated from its isotopes. It is interesting to note that the same method 
of electromagnetic separation is useful in preparing what might be termed 
pseudo-solid sources from gaseous elements. 

Condition (b) above introduces the most severely limiting factor in 
accurate spectroscopy. When the conversion coefficient is extremely 
small the technique applied by Turner and Cavanagh (1951) may greatly 
assist (see also Zeldes et al. 1951). In effect they used a scintillation 
spectrometer as the detector in the magnetic spectrometer. With 
sources giving weak conversion lines great difficulty is experienced due 
to the ‘ background ’ of radiation arising from reflection of the primary 
f-rays for example or from the secondary electrons released from various 
parts of the equipment by y-radiation. By using a scintillation counter 
as the detector those pulses whose amplitude corresponded at least 
approximately to the particle energy as defined by the Hp value in 
question could be selected from the total and considerable improvement 
effected. The use of spectrometers in series is a neat extension of magnetic 
spectrometry. Incidentally it is a point in favour of the ‘ wedge ’ type of 
field for the magnetic spectrometer that it lends itself readily to the use 
of a scintillation system as detector. The detector can be well outside 
the main field. 

Considerable problems arise when external conversion must be 
employed. An account of the main features of the external conversion 
technique has been given by Deutsch, Elliott and Evans (1944). Among 
the difficulties we have the loss of energy of photoelectrons in the 
converter and the presence of a relatively intense Compton continuum. 
Hornyak et al. (1949) have considered theoretically the effects of the 
converter on the energy distribution of the photoelectrons. Recently 
Lind and Hedgran (1951) have given the results of a detailed study 
of the effects of the converter. Little systematic experimental work 
has as yet been done at high resolution. Clearly it is difficult to obtain 
results with a converter approaching in thickness that of a carrier-free 
source with internal conversion lines (~several g/cm’). 
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5.1. Large Double-focusing Spectrometer 


The instrument described by Hedgran satisfies many of the stringent 
requirements. This type of spectrometer was designed and tested by 
Hedgran, Siegbahn and Svartholm (1950) and Hedgran describes in 
more detail how to operate it under optimum conditions. Svartholm 
and Siegbahn showed that with an inhomogeneous magnetic field of 
axial symmetry and specified radial field gradient a semi-circular spectro- 
meter possessing high resolution and high transmission was possible. 
Svartholm (1946, 1950) has examined the electron-optical relations 
involved. In the instrument in question the component H, of the field 
in the median plane takes the form 


Hfr)=Hrw) [1-3("2) +05") +....} ec) 


There is simultaneous first-order focusing in both the radial and axial 
directions after an angle 6=7V 2 (254°). The factor a, may be varied 
to improve the second-order focusing and here a,=3/8. The pole gap is 
28 cm X 28 cm and 7 is 50 cm. 


Fig. 9. 


Source arrangement for external conversion. 


It is found that, if the detector slit is small compared with the image | 
width B, the total relative line width expressed as 4H/H is 


AH/H=Bj(4r,). can oe ae one MEETS 


Here B consists of three components, viz. the sample width By, the 
aberration width By, and the astigmatic line width B,,. For optimum 
results rough equality of the three widths should be arranged. B, does 
not depend much on the height of the sample which is usually 15 mm. 
With a source width of 4mm the contribution to 4H/H is about 0:2%, 
and for a source 2x15 mm? this drops to 0-1%. The component B,), 
depends mainly on the radial aperture y, and the value is (4y,27,)/3. 
Radial shutters are used to adjust the aperture and at the maximum 
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value y;=y, (axial aperture)—0-2 radians, corresponding to a solid angle 
of 477/ 100. At this high luminosity the aberration contribution to line 
width is 13%. At a solid angle of 7/100 however the contribution is 
only 0-:1%. The component B,, is due to deviation of the field gradient 
from the prescribed value and for a relative deviation « is given by 
B=" enV 2. Obviously B,, does not decrease as rapidly with 
aperture as B,). 

Hedgran discusses practical details of procedure which can considerably 
affect the quality of the measurements. At low field strengths there is 
difficulty in reproducing results with accuracy and the range from 
0-3 to 3-0 Mev is considered as well within the compass of the instrument. 
Since most of the accurate work is comparative in nature the design of 
the source holder and support requires special care. The holder can be 
moved on a Wilson seal and the actual form of the holder itself is 
indicated in fig: 9. For external conversion the uranium converter is 
formed by evaporation on to an aluminium plate while for internally- 
converted sources the active material is mounted on a thin aluminium 
foil which replaces the converter. Interchange of sources introduces 
errors of less than 1 part in 104. 

The values of efficiency quoted by Hedgran are very interesting. The 
converters of uranium generally have thicknesses of 0-7 mg/cm? for 
hv<0-5 Mev and 3 mg/cm? for hy>0-5Mev. At these thicknesses the 
multiple scattering is insufficient to make the angular distribution of the 
photoelectrons isotropic. At hy=1 Mev with a thickness of 3 mg/cm? 
about 5x 10-* photoelectrons are released per incident y-quantum. The 
‘ overall efficiency ’, that is, the ratio of the counting rate at the peak on a 
photo-line to the quanta emitted by the source, is about 5 x 10-® when the 
solid angle is 47/100 and the resolution is 0:2%. This indicates how 
valuable the large amounts of shielding that can be used in this big 
instrument prove to be. 

Field measurements. The accurate determination of the field strength is 
critical to the success of magnetic spectrometry and Hedgran states that 
the centre of a line can be located to about 1/10 of the half-width which 
means the energy can be found to some 2 parts in 104. The field value 
should therefore be found with less than this uncertainty, or to about 1 
part in 10. Proton resonance cannot readily be employed at the fields 
required for this big instrument (10 to 200 gauss) and more conventional 
equipment has been brought to a very high level of performance. Cork 
et al. (1947) and Langer and Scott (1950) described very acceptable 
procedures and that employed in Stockholm resembles that of the latter 
workers. The system is shown diagrammatically in fig. 10. Two flat 
coils mounted on a common shaft are rotated with the angular velocity w 
in the spectrometer field, H,, and in the field of a set of Helmholtz coils, 
H,. The field H, must be known accurately in terms of the current [ 
through the coils, ic. H,=F(I). The coils are set so that the induced. 
voltages which are taken off via brushes are as completely as possible in 
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anti-phase. The fields are balanced in this null method by changing I 
to adjust H,. With H, including stray effects and taking into account 
the earth’s field we have 

Hy= 3 H_=const FLY” 3) 2 eee 
where S, and S, denote the effective areas of the coils. The accuracy of 
the measurement of H, thus depends on the precision balancing oi HH, 
against H, and the accuracy of measurement of J. 

Analysis shows that for an accuracy of 1 part in 2x 10* the two coils 
must be adjusted relative to each other to an accuracy of 0-003° and the 
torsional vibration in the rotating system must be less than this in ampli- 
tude. The same battery may be used to produce both fields by having the 
circuits in series and balancing is achieved by adjusting a resistor in 
parallel with one of them. The null position is found with the help of a 
selective audio amplifier (20 cycles/sec) similar to that described by 
Sturtevant (1947). The minimum is defined to about 1 part in 30 000. 


Fig. 10 
Spectrometer Helmholtz coils 
=—| poy 
Ss — 


Potentiometer 


Circuit of magnetic spectrometer. 


The current measurement is effected by comparing with a potentiometer 
the voltage drop produced in an oil-immersed manganin resistance with 
the E.M.F. of a standard cell. All of these points indicate clearly that the 
mechanical and electrical design of the spectrometer reaches a high 
standard of precision. 

Calibration energy standards. Lind and Hedgran (1951) discuss in 
detail the important question of setting up energy standards with the 
large spectrometer. They point out that the success of high-accuracy 
comparative work depending on photo-conversion, internal or external, 
involves at least three stages :— 

(1) Establishment of an absolute energy standard. 

(2) Measurement of standard gamma-rays over the whole energy region. 

(3) Comparison of a standard with the unknown radiation. 
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Initially some difficulty with regard to requirement (1) was experienced 
since the original standard assumed was the energy of the gamma- 
radiation of Au at 411-24+0-10 kev as measured by the bent crystal 
technique. Comparison of the radiation of !8Au and that produced in 
annihilation of positrons with the magnetic spectrometer threw doubt on 
the crystal determination of the first. The energy of annihilation y- 
radiation could be calculated with the help of known atomic constants, 
viz. 


hv=c?(e/m)-1 x 10-4=0-5109689-+0-0000102 ev. 
But in the magnetic spectrometer the K-conversion of the annihilation 
radiation and the Ly, conversion of the 18Au rays with uranium as 


converter gave rise to a doublet of energy difference 4 V=—0-73-L0-04 kev. 
We have 


(Loy —L x) — (Eyy—F,, = A V 


where the # terms are self-explanatory and where H,=115-595, 
Ey, = 17-16, each in kev. Hence 

Ep Ey + (Ex—Hy,,,) + 4V_ or 

Ey=E sy +(99-16+0-04) kev. Roce ice Kip ar Nig ah (a 
With H,,, given the value 411-24 the value for H,,, was 510-4 kev which 
differed too much from the value calculated above, 510-97 kev. Subse- 
quent work by Muller eé al. (1952) revealed certain non-linearities in the 
operation of the bent crystal and the value for #9°Au was modified upwards 
from 411-24 kev to 411-770+0-036 kev. This value substituted in (12) 
above gives H,,—510-93 kev which agrees with, 510-97 kev within the 
experimental errors. Previous to this remarkable demonstration of the 
high accuracy of both bent crystal and magnetic spectrometer, work at 
Stockholm had shown that the latter instrument could derive its own 
absolute standard. The method employed was suggested by K. Siegbahn 
(1944, 1952) and it very neatly avoids the absolute measurement of 
magnetic field strength. A pair of electron lines from the K and L shells 
must be measured but only the ratio of the fields corresponding to equiva- 
lent points on the lines is sufficient. The formulae are 

(H,p)2=A+VA2+B; A=(a2+1)(a2—1)-2(smoce)? ; 
B=(4—s)?(a2—1)-*(smgc/e)*s—? : s= (Ey, —Ex)/(m ec") (¢=H/H,. 

Lindstrém measured the I line of ThB and the L line of ThC" in his 
semi-circular spectrometer and plotted them to coincide. With the 
help of the necessary X-ray data the absolute value for the L line was found 
to be Hp=2607:17-0-30. With this line as the absolute standard the 
values for the annihilation radiation and 1*Au were found to be 
510-96--0-15 Mev and 411-75-L0-10 kev respectively. Clearly at least 
two gamma-radiations are now very accurately known and self-consistent 
go that the second requirement—building up a series of standards—can 
proceed. We can anticipate a big increase in both the amount and 
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accuracy of the data in y-spectroscopy in the near future. As examples of 
the recent work we show in fig. 11 the spectra for #9°Au and the annihilation 
radiation. The first was obtained with uranium converter of thickness 
0-7 mg/cm?, in which only 20% of the photoelectrons suffer inelastic 
collisions within the converter so that the effect on line width is negligible. 
The resolving power is 1-7 parts in 10% and the transmission about 0-1%. 
Other studies by Lind and Hedgran reach resolving powers of about 3 parts 
in 104. The Doppler effect in the annihilation of positrons gives rise to an 
energy spread of 975 ev in the case of copper and the effect is obvious in 
the spectrum. Siegbahn (1952) has suggested the study of this pheno- 
menon in other elements as an informative programme. Some of his 
Fig. 11 


Cc? @ eTt+o” U, 


Intensity 


8650 8700 8750 8650 8700 8750 R 
(2) The Ly; line from uranium for 1°*Au and (b) the K line for annihilation 
radiation. Curves displaced vertically to permit comparison of the three 
representative sets of data. R is potentiometer reading, proportional to 
Hp except for correction for earth’s field. 
general remarks on the future application of such spectrometers are useful. 
He notes that Slitis has examined the width of K-conversion lines (~80 ev) 
at resolutions of 1-2 in 104 and he considers that increased accuracy (by a 
factor ~10) in determining the binding energies of electrons in atoms could 
be achieved with iron-free double-focusing spectrometers. Such a device 
will be described later. 

To conclude the discussion of the Stockholm work in obtaining standard 
energy values for spectroscopy we give in table 1 the results collected by 
Siegbahn (1952). He notes that remarkable consistency in the value of 
the energy of the line of 1°’Cs was obtained in four different measurements. 


5.2. Iron-free Flux Ball Spectrometer. 
We have noted above that the presence of iron is one of the limiting 
factors in obtaining very accurate reproducibility of results in magnetic 
spectrometry. An instrument has been designed recently at Pasadena 
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along rather novel lines and it seems certain that this spectrometer 
possesses a number of very advantageous features among which is the 
absence ofiron. It is described as an axial-focusing beta-ray spectrometer 
of high luminosity, resolving power and precision with proton-resonance- 


Table 1. Calibration Standards 


: Hp E 
Line (ae ae Method Reference 


Th-F 1388-56 +0-15|] 238-63 -+0-04 | absolute Lindstrom 
(1951a, 1951b) 
Th-l 1754-01 40:20] 238-62 -+0-04 | absolute 


Th-L 2607-17 £0-30| 510-85 -£0-80 | absolute 

193Au, K-conv. | 2222-4 +0-04| 411:75+0-10| against Th-L | Lind and 
Hedgran 
(1952) 

13703, K-conv. | 33381:28+0-5 | 661-65 -+0-15 | absolute Siegbahn 
against Th-L | (1952) 
against annih. 

Th-X 9986-7 41-5 2614-25 +0-50 | absolute Lindstrom 

(1952) 
Na, U; 10363 +4 2753-5+1-0 | against Th—X | Hedgran and 
Lind (1952) 


Rays from source through slit and parabolic behaviour of energy ¢ of rays through 


stabilized homogeneous field (Du Mond e¢ al. 1952). Du Mond (1949) 
showed that the design of the spectrometer was optimized as regards 
highest luminosity combined with highest resolution when the quotient 
(Ac/e)/(44)? was a minimum. Here de and « are respectively the instru- 
mental line width and the energy of the primary electron and 4¢ is the 
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utilized fraction of the total sphere about the source. It was shown that 
the mean colatitude angle of the utilized electrons should be about 45° 
relative to the axis of the instrument and that the radiation could readily 
cover the range 45-+10°. This amounted to showing that most lens 
spectrometers were too ‘thin’ and indeed a six-fold improvement is 
effected by increasing the colatitude from 15° to 45°, In _ practice 
Du Mond’s theory suggests that a big wide lensisrequired. He considered 
in detail the distribution in energy of the rays passing through a defining 


slit P as shown in fig. 12 and found that it was parabolic. The effects on’ 


the line profile of (1) spherical aberration, (2) size of the annular resolving 


Fig. 13 


Conical source 
semi-apex angle 23%48’ 


| ‘2 5 4 
Component pulse profiles obtained with conical source and modified slit (above), 
due to spherical aberration, 1; main resolving slit, 2; finite size of cone 
source, 3; beta population with retardation at source, 4. 


slit, (3) size of the disc or cone-shaped source, (4) retardation in the source 
were studied. All of these were found to possess the very important 
property of maintaining a first-order discontinuity on the high energy 
edge of the profile with the exception that the effect of the disc source 
was to make the profile continuous. The cone-shaped source was 
suggested by D. A. Muller and is such that the semi-apex angle is chosen 
equal to the obliquity of the sinusoidal ray traces for the ultimate ray 
in contact with the inner jaw of the principal annular resolving slit: 
Figure 13 illustrated the idea and shows schematically the predicted 
shapes of the component line profiles (for cone-shaped source). 


ee 
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Resolving power. Th-B. was deposited on a cone-shaped source of 
semi-apex angle 23° 48’ and the predicted linear decline on the high energy 
edge observed as seen in fig. 14. The point Q refers to the field value at 

Fig. 14 
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Nominal oscillator setting. 
Internal conversion beta-ray profile for electro-deposited thoron coating. 


which the ultimate ray of the pencil just clears the inner jaw of the 
defining slit. Location of Q to about 1 part in 104* seems possible. It is 
claimed that Q can be located accurately without using thin sources and 
that converters in the required conical form can be fabricated and the 
y-ray source placed within them. The technique is regarded as a good 


Table 2. Matched Parameters for Spectrometer 


Fraction of sphere, 4 & 0-005 0-020 
e>1: Ae/e 0-000126 0-00202 
e=1: Ae/e 0-000189 0-00303 
eX1: Ac/e 0-000252 0-00404 

Angular opening, 40 deg. 48-5/ 3° 15' 


Resolving slit width, 4a cm. 0-0101 0-1620 
Source radius, cm. 0-00449 0-0719 
Source disc area, cm? 4-99 x 10-5 1-62 x 10-? 
Source cone area, cm? 1-235 * L041 4-01. 10-2 


one for analysis of Compton recoils. Theoretical considerations show that 
three factors, viz. the utilized angular aperture at the source, the annular 
resolving slit width and the diameter of the cone base, must be matched for 
the best results. Some calculated values are entered in table 2. Here 
the resolving power is given in terms of Ae/e (the relative line width at 
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half-maximum height=1-62 4e/e). With ¢« in mc? units the cases 
«<1, e=1 and «>1 have been computed. Preliminary tests show that 
the calculations are realistic. A maximum value of 20cm for the 
departure of the B-rays from the axis applies here. 

The tabulated data are very impressive and lend strong support to 
the belief that the device will prove of great value in spectroscopy. Some 
of the design features deserve comment. The homogeneous field. is set 
up within an assembly of forty coaxial circular coils of diameters 
graduated so that the mean turns of the coils lie on an ellipsoid of 
revolution. The spacing of the coils is such that there is everywhere the 
same number of ampere turns per unit axial length. The design is 
regarded as a variant of the Maxwell ‘flux ball’ and appropriate 
theoretical treatment of ellipsoidal windings has been given by Blewett 
(1947). A special Geiger counter is used as the detector. Its perforated 
wall gives about 70°% transmission (mica } mil thick acts as a seal over 
the cylinder). The magnetic field control is very interesting. It is 
automatically stabilized against rapid fluctuations by an electronic 
stabilizer and against slower drift by three proton resonance heads 
permanently installed in the assembly and arranged to give an overall 
range of control from 35 to 1200 gauss. Special small Helmholtz coils 
(2 pairs) are provided around each head to supply the a.c. sweeping field, 
which is of sufficient amplitude to cover the entire width of the proton 
resonance signal. A solution of MnSOQ, in water (1/100 molar) is used and 
the R.F. field is set up by a coil of wire wound round each bulb. The 
power is supplied by a d.c. generator driven by a synchronous motor. 
The proton resonance device forms an element in a feedback (servo) 
system giving an easily measured reference, i.e. the frequency of the 
R.F. oscillator. The resonance serves to measure the strength of the 
magnetic field and to stabilize it with precision. The d.c. generator is 
made into a power source of low noise and controllable output by 
incorporating a subsidiary inverse feedback system as described by 
Sommers et al. (1949). Proton resonance as used here is based on the 
method of Pound and Knight (1950) and it is an interesting example 
of the application of basic studies of nuclear magnetic resonance and 
nuclear induction (see Bloembergen et al. 1948, Bloch 1946). The 
proton has two energy states corresponding to the two possible 
directions for its spin axis. In each the spin precesses about the 
direction of the magnetic field with the Larmor frequency w=ypH 
where yp/27=4-2578 ke/gauss. Absorption of energy from a radio- 
frequency source can occur each time such a source produces a field 
passing through the resonance value and demodulation gives an audio 
signal which can be amplified and used for control. 

Du Mond notes the very good agreement among some recent measure- 
ments on the I and F lines of thorium. The ratio of Hp values for these 
two lines for three different experiments (the first refers to the new 
helical instrument just discussed) are (Hp),/(Hp),»=1-26329-+0-00006, 
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=1-26317+-0-00012 (Craig 1952)—1-26319-+-0-00005 (Lindstrém 1951). 
These results amply demonstrate the power of the new magnetic 
techniques and show conclusively that the bent crystal is required at 
most for supplying some standards of energy value or checking those 
secured by other means. It seems certain to be relegated to this minor 
role. 

Some comparisons of magnetic instruments. Du Mond et al. have 
tried to compare their new helical instrument with a number of fairly 
recent types. For this purpose they use the dimensionless quantities 
® and Y suggested by Persico and Geoffrion as figures of merit. For 
helical instruments with ring exit slit these are given by 


Cae ane ne lag i o— (18) 
where »=Ap/p is the spread in the momentum scale at half-maximum, 


w=the gathering power, L=ow with o=source area in cm? and 
/=distance in em from the source to the detector. 


Table 3. Comparison Data for Magnetic Spectrometers 


3 : 
Reference Iron ee vy o eee @ x 104 Polo? 
Deutsch, Elliott and no | 100 1-7 0-18 0-22 5:8 1:4 
Evans (1944) 
Siegbahn (1946) yes | 50 5-0 2-00 10 72 9-0 
Hornyak, Lauritsen and} no | 114 1:5 0-50 3°55 100 4-1 
Rasmussen (1949) 
Jensen, Laslett and no | 102 2:3 0:34 1-22 14-6 2-2 
Pratt (1949) 
ee rae | ela OF ee Ox NOR es Sc 1102 
% em? 
Agnew and Anderson no | 96 9-0 11-0 585 260 
(1949) 
Slatis and Siegbahn yes 4-0 8-0 16 160 
ae (1950) 34 9-0 0-31 10 35 
Feldman and Wu (1 no : : 
Du Mond e¢ al. (1952) no | 62-8 | 10-5 5:63 | 15-0 x 108 1065 11-0 


\ v 1 
0-032 | 0-338 |7-89x 10-7 


The values ¥* =T'max/n1/? are such as to take instrument and counter losses 
into account and ©*=(o/7/?) Y refers to a conical source with base diameter 
matched to the other parameters. A wide variation of 7, @ Or T nax and L or L* 
is possible by varying source size, acceptance angle slit jaws and annular | 
resolving slit jaws. 


$6. Tae Bent CrySTAL SPECTROMETER 


The second of the high accuracy instruments, the bent crystal, has 
been: discussed recently by Sandstrém (1951), with reference to soft 
X-ray measurements. It is comparatively recently that it has been 
used by Du Mond for the investigation of hard y-radiation although its 
high resolution combined with its ability to operate at low intensities 


2H2 
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sts its suitability for this purpose. The early names associated 
hee serail ates Du Mond and Kirkpatrick (1930), Cauchois 
(1932), Johann (1931) and Johansson (1933). The Pasadena school has 
employed the instrument in a number of interesting investigations 
(Du Mond et al. 1947, 1948, Du Mond 1947, 1951, Lind et al. 1949, 1950, 
Muller et al. 1952). The principle of their equipment is indicated in 
fig. 15. It operates by transmission with R and V the real and virtual 


Fig. 15 
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Collimator 


Defector 


Arrangement of bent crystal spectrometer. 


foci. The Bragg angle is denoted by @ and B is the point at which the 
reflecting atomic planes converge if produced. The instrument has a 
radius of about 1 metre for y-ray studies. The crystal is a plate of 
monocrystalline quartz 2mm thick which in the unstressed state has 
accurately optically flat plane parallel faces, cut so that the 310 planes 
are normal to the flat faces. The sheet crystal is bent to the required 
curvature by compressing in a vice. 

The source is placed at the focus and arranged so as to maintain the 
selectively reflected beam stationary while the range of wavelengths is 
investigated by moving the source, A lead collimating system on the 
side of the crystal remote from the source serves to shield the detector 
from the direct beam of y-radiation while transmitting the reflected beam. 
The displacement of the source-supporting carriage must ‘be accurately 
proportional to the wavelength. In the Pasadena equipment a neutron- 
activated source of 'T'a—strongly emitting the x-rays of its daughter 
element 16*W, was found to perform well as a calibrator. 
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The most satisfactory detector takes the form of a large crystal of 
Nal (Tl) acting as a simple scintillation spectrometer. It is viewed with 
a pair of photomultipliers acting in coincidence and it selects the required 
radiation from the background and makes the device much more sensitive. 

Among the sources successfully investigated were 1°8Au, Cu (annihila- 
tion y-rays), °1I—3 lines, ®°Co—2 lines, 18?Ta—19 lines, radon and radio- 
thorium. The upper limit of the energy scale was estimated to be at 
1-5Mev. With this instrument D. Muller was successful in measuring 
each of four y-ray lines in the first, second and third orders of Bragg 
reflection. These measurements, given in table 4, show a small systematic 
departure from linearity in the spectrometer and a correction curve is 
required. 


Table 4. Higher Order Reflection 


A in 10-1! em 


Energy 

k 

oe lst order 2nd order 3rd_ order 
295-79 41-909 41-894 41-898 
308-26 40-214 40-205 40-191 
316-28 39-197 39-183 39-174 
467-53 26-515 26-503 26-496 

191-025 190-993 oe 


§7. Tue Parr SPECTROMETER AND GENERAL REMARKS 


No precision method of analysis such as those described above has been 
applied to y-radiation of energy exceeding 3 Mev. It is unlikely that the 
methods described above, and which depend on photoelectric interaction, 
will prove useful at higher energies. As we have seen the pair process is 
clearly suitable in this higher range but it does not lend itself readily 
to high accuracy work due to the continuous nature of the energy 
distribution of any one of the electrons. A method of measuring the 
energy of both particles simultaneously was introduced by Walker and 
McDaniel (1948) but it does not seem capable of acting as a high resolution 
device. Two sets of counters are set on opposite sides of a thin radiator 
of heavy material and associated electronically so that the sum of the 
energies of the electrons is effectively measured as a function of the 
magnetic field. In this way the pair process is used to give rise to peaks 
on the spectrum but the slits of the spectrometer are essentially the 
counters themselves. This means that at best we have a spectrometer of 
intermediate resolution. An excellent illustration of the application of 
this type of instrument is found in the work of Kinsey, Bartholomew and. 
Walker (1951, 1952) on the process of capture of neutrons. The sample is 
placed in a high neutron flux near the reacting core of the Chalk River 
pile. A block of bismuth, 5 in. thick, is placed between the sample and 
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the hole to reduce the intensity of the pile y-rays. Lead collimators in 
the hole limit the y-rays from the sample to a beam which is allowed to 
fall on the radiator of the spectrometer. The radiator is generally gold or 
other heavy element about 5 mg/cm? in thickness. The line width in 
this arrangement is usually 130 kev but may be as low as 65 kev. The 
efficiency rises roughly as the fifth power of the quantum energy above 
3mey. The instrument performs rather satisfactorily in the range 3 to 
about 10 mev and it is suitable for the study of nuclear reactions. A high 
precision spectrometer for this kind of investigation is still required. 
An example of the work is shown in fig. 16. 


Fig. 16 
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Corrected gamma-ray spectrum from chlorine. Line width, 130 kev. 


General Remarks. The bulk of the above survey has been concerned 
with a limited number of new instruments which seem to the writer to 
demonstrate adequately the main trends in y-spectrometry. Many 
matters have been omitted or considered at best very briefly. Fortun- 
ately an excellent and rather comprehensive treatment of most of this 
material has been supplied by C. S. Cook (1953). He considers in turn 
the source, the vacuum system, the detector and the magnetic field and 
each section contains many references to original papers. Here we 
will limit ourselves to a few remarks on the counter. Until recently it 
was extremely difficult to employ a Geiger counter for detecting electrons 
of energy less than about 25 kev. The new materials like nylon and 
formvar, which can be prepared in thin layers (see Richards et al. 1948, 
1949, Zworykin et al. 1945, Langer 1949) appear to act well as windows 
and transmission down to about 5 kev is achieved. Charging of such 
windows must be prevented and the rise in transmission between cut-off 
and complete detection gives rise to uncertainties. Some transmission 
studies have been made by Saxon (1951) and Heller et al. (1950). Good 
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progress has been made nevertheless in observing low energy photoelectric 
and Auger lines. Attempts have been made to solve this same difficulty 
at low energies by accelerating the particles both at the gap between the 
exit slit and the counter (see Langer and Cook 1948, Agnew and Anderson 
1949, Butt 1950) and at the source (Chang and Cook 1952). Townsend 
finds that some of the spurious electrons which cause difficulty in this 
technique are of thermal origin while Chang obtains evidence of the 
presence of low energy secondaries emitted from the source and its 
support. Another possible answer to this problem is to be found in the 
use of the scintillation detector. Although the detector is windowless 
there remains the difficulty of ensuring that it is 100° sensitive to low 
energy electrons and scattering complicates intensity measurements. 
Sometimes a light guide must be employed and this reduces the size of the 
signal. Butt (1950, 1952) has discussed his work with anthracene as 
scintillator and it may be possible with this material to operate down to 
about 5 kev. It would seem preferable in view of the better conversion 
efficiency of Nal(Tl) to use this phosphor as the detector even if the 
scattering is increased. A special geometry of the crystal or crystals could 
help considerably (see fig. 17) in eliminating this source of trouble. 


Big. 17 
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Slotted crystal C behind aperture slit S. 
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§1. INTRODUCTION 


THE present review is concerned with the physics of a class of substances 
which have in common certain features of bonding, a peculiar variability of 
structure in the solid state, and a liability to order—disorder phenomena 
of different kinds. The materials in question are normal long-chain 
paraffins and certain of their derivatives. They are closely related to 
plastics such as polythene, polyvinyl chloride, etc., but the present 
review is restricted to ‘ short’ long-chain substances of exact chemical 
formula and less than some 100 carbon atoms, which are crystalline 
when solid. 

It is an essential feature of long-chain substances that they form 
homologous series, whose members differ from one another by small and 
distinct steps. For instance, the chemical formula of any normal long- 
chain paraffins can be rendered as CH;(CH,),.»CH, (or, more simply, 
C,,H,,.), where » is the number of carbon atoms in the molecule. In 
some cases, a homologous series implies two or more serial numbers, as 
is the case with long-chain ketones whose formula may be written 
Cine ym1 - CO. Cr mH on—m)i1 Where n denotes again the number of 
carbon atoms and m denotes the position of the CO group on the chain. 
With regard to physical properties, the repetitive nature of the chemical 
constitution has the effect of introducing an unusual kind of independent 
variable. Within certain limits it is often possible to consider a homologous 
series as an entity whose properties can be treated as depending not only 
on pressure and temperature, but also on the serial numbers. 
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Interest in the physics of long-chain crystals arose when it was realized 
that solid long-chain substances exist in a variety of structural modifi- 
cations, or polymorphic forms, some of which are disordered, and that these 
depend in a more or less systematic way on serial numbers. The physics 
of these materials was first treated systematically by A. Miller. In a 
series of papers published between 1928 and 1941, he investigated the 
crystal structure as well as the bonding and physical properties of typical 
paraffins and their derivatives, and presented a general picture which is 
still valid in its essentials. 

In more recent years, a fair amount of work was done on long-chain 
crystals, the impulse for which came mainly from two directions : interest 
in order—disorder phenomena, and secondly, interest in dielectrics. 
From the latter point of view, long-chain crystals are of interest because 
dipoles within them may be mobile, so leading to comparatively high 
dielectric constants. 

The disorder occurring in long-chain crystals is of two kinds. Firstly, 
transitions occur from ordered forms, stable at comparatively low 
temperatures, to high temperature forms where molecules turn and twist 
while being part of a solid lattice. Secondly, cases occur where the repeat 
of the lattice is irregular ; for instance the lattice may contain holes or 
molecules may be ‘ upside down ’. 

Transitions in solids to forms with ‘rotational freedom’, i.e. with 
extensive mobility of molecules about their sites in the lattice have been 
reported for many materials. In particular, Smyth and co-workers 
and Perdok and Terpstra (1947) investigated numerous solids consisting 
of approximately spherical molecules, mainly methane derivatives. 
A good review of the fundamentals of ‘rotational transitions ’ in general is 
given by Eucken (1939) and this is relevant to such transitions in long- 
chain crystals as a special case. However, the repetitive nature of the 
long chains implies many features which justify the treatment of long- 
chain crystals as a special subject. 

In the following, the structure of long-chain molecules and the forces 
bonding them into crystals will be discussed, and it will be shown how the 
structural features lead to the different observed order—disorder phenomena. 
The exposition aims at a unified presentation, in terms of physical 
concepts, of a subject which belongs partly to the realm of chemistry. 
The diffuse boundaries of the field surveyed make it unprofitable to 
attempt a complete bibliography. 


§2. STRUCTURE AND BONDING 
2.1. The Structure and Flexibility of Long-Chain Molecules 


The formation of long-chain molecules is possible when atoms are able 
to form at least two essentially homopolar bonds. In paraffin deriva- 
tives, chains consist mainly of carbon atoms, and their properties are 
determined chiefly by the valency characteristics of carbon. 
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Carbon has four electron orbitals which take part in bonding and has 
a valency of four in almost all of its compounds. It is well known that 
carbon is capable of forming homopolar, single, double and triple bonds 
and systems of conjugated single and double bonds, as in benzene. The 
formation of these bonds can be understood in terms of quantum mechanics, 
although calculations have to rely on approximate methods which yield 
accurate results only in favourable circumstances. The quantum 
mechanics of homopolar bonding is outside the scope of the present 
paper. However, it will be useful to refer briefly to the properties of 
the single bonds formed by carbon, since these are the basis of the 
structure of paraffins. 

The properties of single carbon bonds in hydrocarbons have been 
reviewed recently by Lennard-Jones and Pople (1951). Of the four 
outer electrons of carbon three are in 2p orbitals and one is in a 2s orbital, 
but the carbon atom forms four equivalent single bonds of hybrid charac- 
ter. When these four equivalent bonds are formed with four hydrogen 
atoms, in methane, the resultant electron configuration resembles that of 
a noble gas. The four electron pairs formed tend to be as far from one 
another as possible and in consequence, the carbon atom in methane is 
in the centre of a regular tetrahedron, the hydrogen atoms being at the 
four corners and the angle between a bond and any of its three neigh- 
bours being 109° 28’. 

The tetrahedral bond angle appears not only in methane, but, with 
slight variations, also in all other carbon compounds which contain 
carbon bonded to four other atoms. The bonds in question are very 
strong, and the bond angles are maintained by powerful forces. The 
dissociation energy for carbon bonds is very high as compared with RT’ 
at room temperature, of the order 50-100 kcal/mol. As for the forces 
that maintain the bond angles, these are most conveniently characterized 
by the vibration frequencies appertaining to vibrations where the bond 
angles are slightly bent out of their stable position, since the frequencies 
of harmonic vibrations are related to the restoring force acting against a 
small deformation. 

The vibration frequencies characteristic for the stretching and bending 
of carbon bonds are known with considerable accuracy from spectro- 
scopic measurements. For instance, according to Pitzer (1937) the 
frequency corresponding to the stretching of a C-C bond is around 
1000 cm~!, while the bending frequencies of H—C-H bonds are around 
1440 cm™ in different hydrocarbons, other bending frequencies being 
mostly lower. The vibration frequency is not unambiguously related 
to the dissociation energy, since the former characterizes the resistance 
of a bond to a small deformation, while the latter is the energy required 
to break it completely. Nevertheless, the high dissociation energy of 
carbon bonds and the similarity of frequencies for bending and stretching 
vibrations indicate that not only the bond lengths, but also bond angles 
are maintained by strong forces. 
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A complete knowledge of all the intramolecular vibrations and rotations 
of a polyatomic molecule implies a knowledge of the partition function of 
the gaseous material and thus of all the thermodynamic data of the gas. 
The spectroscopic data available and the fact that the paraffins consist of 
simple repeating units has been utilized to calculate partition functions 
for them. It is not the object of the present review to discuss this work in 
detail, but only to summarize it briefly ; a general review is given by Fowler 
and Guggenheim (1939). 

The bending and stretching of bonds is not the only way in which paraffin 
molecules may be distorted, but they may also rotate about single bonds. 
It was realized at an early date that this movement is likely to be opposed 
by weak forces only, since it does not involve a change of either bond 
lengths or angles. In constructing partition functions, assumptions have 


g - 0° (trans) 


Three successive bonds of a long-chain molecule, illustrating 
rotation about the central bond. 


to be made about the energy conditions involved in rotation about bonds, 
and various models have been verified by comparison with the experimental 
data furnished by spectroscopy and calorimetry. The main features of a 
model which fits the facts are by now well established largely due to the 
work of Pitzer, but details are still under discussion. The subject has 
been reviewed by Pitzer (1951) and Aston (1951). 

The model evolved by the study of gaseous paraffins gives very clear 
data on the flexibility of the molecule. It will be described with the help 
of fig. 1 (according to Taylor) which shows a chain of four carbon atoms, as 
in normal butane, CH,-CH,-CH,-CH;. The four carbon atoms in 
fig. 1 are connected by bonds of constant length o=1-54 A, forming 
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tetrahedral angles with one another. The constancy of bond lengths and 
angles permits all configurations for which the two outermost carbons are 
placed along the two circles drawn. 

The energy of the different possible configurations shown in fig. 1 is 
known to have minimum values and the relative energies may in first 
approximation be described according to Pitzer by an equation for the 
potential energy V of a configuration with the angle ¢ in the figure 


V(¢)=4 Vy (1—cos n¢) ee ey sone ON 


where n for ethane or long-chain paraffins is taken as equal to three. The 
height V, of the barrier is around 2750-3000 cal/mol., and thus appreci- 
able as compared with k7’. Equation (1) is on the whole satisfactory in 
accounting for the specific heat and spectroscopic data, but it has to be 
refined to account for the average shape of long-chain paraffins. Under 
the effect of heat motion a long flexible chain in the gas phase will clearly 
be bent or coiled, to an extent that is calculable from eqn. (1). However, 
the result of the calculation leads to too strongly coiled configurations, so 
that the equation has been modified by Taylor (1948) to 


V(¢)=4 V,[X (1—cos ¢)+(1—X)(1—cos 34)] . . . . (2) 


where X is a constant between zero and unity. For X=—0 the formula 
reduces to (1), while for X 40 it implies that the configuration for which the 
chain is plane is more stable than the two alternate ones. Taylor quotes a 
difference of 800 cal/mole. between the stable and the subsidiary minima, 
but Ubbelohde and Woodward (1952) consider this too large. However, a 
small difference seems generally accepted. 

Equation (2) implies that the configuration of lowest potential energy of 
a paraffin molecule is uniquely determined as a flat zig-zag, as shown 
in fig. 2, which represents the crystal structure of a solid paraffin. In the 
gas phase, the molecules deviate from this configuration, but according to 
Taylor (1948) the crumpling of short chains is not very marked at 
room temperature. 

While data for normal paraffins are remarkably complete, those for 
derivatives are less so. In the following, a brief description will be given 
of the derivatives which come within the scope of the present review. 

Chain molecules containing double bonds in even one of the skeletal 
links differ greatly from normal paraffins in their flexibility. Such 
molecules, and others with too different a flexibility will be excluded. 
This leaves substituted paraffins, where the carbon skeleton is unchanged, 
but one or more hydrogens are replaced by other atoms or groups, and 
esters, where a carbon in the chain is replaced by an oxygen. In the 
first case the flexibility may be slightly affected by the bonding or size of 
the substituted group. In the second case, the oxygen link in the skeletal 
chain probably does not greatly affect the flexibility, since the bond angle 
between the two bonds of oxygen is very close to the tetrahedral angle. 
For each paraffin of a given chainlength there exists a large number of 


a 
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derivatives of the kind considered here, according to the kind of substituted 


group, the number of groups on the chain, and the link or links at which 
the groups are attached. 


2.2. The Crystal Structure of Long-Chain Solids 


The crystal structures of all crystalline paraffins and derivatives have a 
great deal in common. In all long-chain crystals the chains pack as 
straight rods, in layers. However, within the bounds of this similarity, a 
great many structural variants exist, only a fraction of which has been 
investigated in any detail. 

The determination of the structure of a long-chain crystal by x-ray 
(or electron or neutron) diffraction is not an easy task. It is obviously best 
carried out with the help of single crystals, but the growing of good single 
crystals of paraffins is a difficult art. On the other hand, powder 
photographs, which are easy to obtain, are difficult to interpret with 
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certainty because of their complexity. For these reasons, most structure 
determinations are more or less incomplete. Some authors, notably 
Piper and co-workers, concentrate on rapid methods whereby they 
obtain surveys of homologous series, without attempting a detailed 
description of any particular structure. Other authors determine in- 
dividual structures in more or less detail. So far, no structure determination 
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has been published of the high accuracy now attainable in x-ray analysis, 
but such a determination on a normal paraffin (C3,H,4) is in progress 
(private communication by Dr. Vand). 

The fundamental features of long-chain crystals may conveniently 
be described for the simplest and most important structural variant, 
namely the orthorhombic structure observed for many normal paraffins. 
This was determined first by Miiller for CygH,o, and his results are shown 
in fig. 2. The unit cell contains four molecules in two layers, in the most 
stable, flat configuration of the chains. Within any layer the chains are 
parallel to one another, their cross-sections being arranged in a face centred 
pattern. 

The structure of C.,H¢. is comparatively complicated because of the 
existence of two layers (drawn in full and dotted lines respectively) 
each layer being staggered with regard to the layers above and below. 
This complication is due to the finite length of the chains and the need 
for most convenient packing of the chain ends. The simplest structure 
would be for infinitely long straight chains, and this is realized, approxi- 
mately, for the crystalline regions in polythene, which consists of paraffin 
chains of some 1000 carbon atoms. Bunn (1939) finds that the structure of 
polythene closely resembles that of a layer shown in fig. 2. 

The electron density maps for polythene, according to Bunn, show loose 
bonding between chains in the basal plane (fig. 3 (a)) and strong bonding 
within the chain itself (fig. 3 (b)). This is characteristic for all long-chain 
crystals. Bunn finds anomalies of the shape of the CH, group, but his 
material, a polymer which is only imperfectly crystalline, is probably not 
quite suitable to be used for a decision on a subtle point of structure 
of this kind. According to Vand (private communication) the electron 
density maps for C,,H,, show a normal shape of the CH, groups. 

The different structural variants of long-chain crystals may conveniently 
by described with reference to the orthorhombic structure discussed 
above. In particular, Miiller’s nomenclature of crystallographic axes is 
always used to describe long-chain crystals, c being the long-chain axis 
and a b the ‘ basal plane ’. 

In classifying the structures observed, it is convenient to distinguish 
differences in the arrangement of layers, and in the arrangement within 
: sub-cells ’,a term due to Vand (1951). As for the arrangement of layers, 
unit cells may contain one, two or, occasionally, more than two layers. 
Chains within any one layer are almost always parallel to one another, but 
chains within alternate layers need not be. Chain axes may be either 
normal to the basal plane, so that the distance between layers (described ; 
erystallographically as 1/c*) is the full length of the chains, plus the width 
of the gap between subsequent layers, or else they may be inclined to the 
basal plane, when the distance between layers is shorter than corresponds 
to the length of a chain. Here it is necessary to distinguish between the 
inclination of chains and the inclination of the crystallographic c axis. 
The latter, which may be chosen in different ways, for a given structure, 
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need not be parallel to the chains in any one layer, and its inclination 
with regard to the basal plane has little physical significance. 

The details of the arrangement of chains in layers essentially depend on 
the end groups of the chains, and are of less intrinsic importance than the 
mutual arrangements of chains within a layer. If we disregard the end 
groups, a layer of chains can be regarded as composed of repeating units, 
‘ sub-cells ’ according to Vand (1951). For simple structures, one sub-cell 
is sufficient to describe the arrangement of matter within a layer. For 
more complicated ones, more than one kind of sub-cell may be present. 

The structures of long-chain crystals investigated so far may be des- 
cribed in terms of three types of sub-cells, although other types may yet be 
discovered. By far the most important type of sub-cell is that common 
to the structure of Cy H¢o (fig. 2) and polythene (fig. 3). This sub-cell is 


Fig. 3 
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(a) 
(a) Electron densities on plane perpendicular to c-axis, at 3/12 c and 
(rotated 180°) 9/12 c. 
_(6) Section through a molecule. 


characterized by a rectangular basal figure a b, to which the chain axes are 
normal. This figure contains two molecules (at the corners and centre) 
the planes of zig-zag of which make angles of +-¢ and —¢ respectively with 
the a-axis. 

The Miiller type of sub-cell has a number of variants which may be 
described as distortions of the basal figure. For some variants, the figure 
a b is rectangular, but the angle 6 between its diagonals varies. If 6=60, 
the basal figure is a regular hexagon, a most important special case. 
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The basal figure may in some cases be also distorted so that the rectangle 
a b deforms into a parallelogram with unequal diagonals. 

A second type of sub-cell is described by Vand and Bell (1951) for 
trilaurin, and probably occurs in the triclinic form of the normal paraffin 
C,sH3, (Miiller and Lonsdale 1948). This is similar to the sub-cell already 
described, but all the planes of zig zag are parallel to one another, while 
making an angle 4 with the a axis. This sub-cell occurs also for C.,H;5, 
and other paraffins. 

A third type of sub-cell was found by Vand, Lomer and Lang (1951) 
for potassium caprate. This differs from all others in that here chain 
axes (as distinct from zig-zag planes) are not all parallel to one another, 
but that they cross one another within the same layer, so that the zig-zags 
of neighbouring chains fit into one another. This form was found only in 
potassium soaps e.g. potassium caprate (Vand, Lomer and Lang 1949) 
but may occur also in other long-chain crystals. 

The nomenclature used in the above survey, as in terms of chain tilt 
and the different sub-cells, is not the same as that used in large parts of 
the original literature, where structures are given according to their 
crystal class. However, little ambiguity exists, as a rule, with regard to 
the most symmetrical structures, the orthorhombic one shown in figs. 2 and 
3, and its hexagonal, still more symmetrical variant. Monoclinic structures 
correspond often to forms with Miiller’s sub-cell with rectangular basal 
figure, and a more or less slight tilt of the chains, and need not be very 
different from orthorhombic forms. However, in general monoclinic and 
triclinic forms may have sub-cells of all types with various angles @ and 
angles of tilt of the chains. 

As regards the distribution of the different structural variants among 
different long-chain crystals, data are not complete, even for normal 
paraffins. However, the study of homologous series, initiated by Piper 
and co-workers, show that for any such chemical series the same variants 
reappear for different chainlengths. In particular, Piper and co-workers 
find that the long-chain paraffins, acids, ketones, alcohols, etc., can be 
divided into crystallographically homologous series, for which the distance 
of layers is 1/c*=k,n+k., k, and ky being constants. For instance, 
for paraffins, Piper, Chibnall and others (1931) find series 4, B and 0 
where /, and k, are constant over wide ranges of chainlength, and Francis, 
Collins and Piper (1937) distinguish series of acids which run from n—16 
to n=46. Vand, Aitken and Campbell (1949) point out that small 
variations of sub-cell or tilt of chains may occur within series of this kind, 
but their careful investigation of seven silver salts of fatty acids shows 
that the concept of a homologous series of crystals is yet valid. 

Of homologous series investigated, normal paraffins will be discussed 
later. Acids and soaps tend to have complicated structures because of 
their ionic endgroup, while ketones are simpler (McArthur 1944) since 
they all appear to crystall’ze in the orthorhombic form. In this review, it 
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is not possible even to attempt a systematic survey of the crystallography 
of all long-chain crystals investigated. The reader is referred to the 
original literature, and to ‘ Structure Reports ’ which review all structure 
determinations year by year. 

All the different variants of sub-cells and of arrangements of layers, 
combined, make up a great variety of structural forms. Nevertheless, 
if we neglect the rare type of sub-cell with crossed chains, all the structures 
have a great deal in common. The chains pack in layers as straight 
rods, and the closest approach of neighbouring chains varies within 
fairly narrow limits, being around 4-5 4. The simplest structure which 
exemplifies the basic type, is the hexagonal structure obtained from 
Miiller’s sub-cell with equal diagonals and 6—60°, where chains are normal 
to the basal plane defined by the endgroups. This hexagonal form 
corresponds to the best packing of cylindrical rods. 


2.3. The Bonding of Long-Chain Crystals 


The basic features of the packing of chains in long-chain crystals can 
be explained by van der Waals bonding between CH, groups. Miiller 
(1936) succeeds in calculating the van der Waals potentials, for the general 
case of a chain of m links, and finds fair quantitative agreement with 
experiment. The van der Waals potential, i.e. the energy to remove a 
molecule out of the range of the attractive forces of its neighbours, is 
between — 2-58 x 10-* erg and — 1-52 x 10~18 erg per CH, group, according 
to different approximations used, while the heat of sublimation (if cor- 
rections for changes of specific heat are neglected) is —1-62 x 10-8 erg 
on the basis of experimental data. (More recent data, quoted in table 2, 
give only —1-2x 10-} erg for the heat of sublimation. However, this is 
not of major importance.) 

The success of Miiller’s calculation shows that van der Waals forces are 
almost exclusively responsible for the bonding of the solid normal paraffins. 
This is not surprising if one bears in mind that paraffins have electron 
shells analogous to those of noble gases. 

De Boer and Vand (1947) show that not only the bonding of normal 
paraffins, but also that of fairly complicated derivatives, namely soaps, 
may be treated quantitatively, as a combination of van der Waals 
bonding between the chains and ionic bonding between endgroups. 
Although the forces of the endgroups are large compared with the van der 
Waals forces, the packing of the chains is not much influenced by them. 
This is due to the fact that the endgroups are arranged in two dimensional 
arrays, so that van der Waals and ionic forces act in separate layers of the 
structure. De Boer and Vand show that for paraffin, where van der Waals 
forces only operate, the cross-section of the chains increases with increasing 
chainlength, but that in silver and sodium soaps the ionic forces cause a 
contraction which is the less the longer the chains, so that here cross- 
sections increase with chain lengths. 

212 
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2.4. Macroscopic Properties 


The layer structure of long-chain crystals and their bonding by van der 
Waals forces determines a number of their macroscopic properties. The 
forces holding layers together are very weak, even compared with the 
weak forces within layers (unless layers are cemented by ionic forces, as is 
the case with every second layer in soaps). In consequence crystals 
cleave very easily along layers and crystal growth occurs predominantly 
parallel to the layers, in the a 6 planes. Paraffin type crystals are almost 
always flat plates or needles, and the crystals often have a rhombic shape. 
These characteristics are strikingly demonstrated in fig. 4 (Plate VII), an 
electron micrograph of a rhombic crystal whose growth spirals show the 
comparatively slow growth in the ¢ direction (Dawson and Vand 1951). 
Mechanically, paraffin type crystals are soft and deform particularly easily 
by glide along layer planes, that is, they are ‘ waxy’. 

Many physical properties are very anisotropic, as has been shown by 
Miiller (1941) for the compressibility. X-ray measurements on crystals of 
C,,;H,, show that the compressibility 

Lick 

= x — 

Peete 
(p being the pressure and / a length) in the a direction is 10 x 10~!cm?/dyne, 
that in the b direction 11 x 10~!* cm?/dyne, while that in the c direction is 
only 3x 10-8 cem?/dyne. Miiller is able to calculate the compressibilities 
theoretically, on the basis of van der Waals bonding. Qualitatively, the 
an‘sotropy in question is due to the fact that compression in the chain 
direction implies a deformation of the carbon bonds, while compression 
normal to the chains is opposed only by the weak van der Waals forces. 
However, it should be noted that Miiller’s observations refer to paraffins 
near their melting point, where their thermal expansion is abnormally 
large (see later). In temperature regions of normal thermal expansion the 
anisotropy of the compressibility may be considerably less. 


2.5. Long-Chain Materials in the Liquid State 


Liquid long-chain paraffins are structurally remarkably similar to solids. | 
Various aspects of this characteristic have recently been discussed by 
Moore, Gibbs and Eyring (1953), and will be quoted here only briefly. 
The chains in the liquid appear to be approximately straight, though not 
arranged in regular layers. The packing is similar to that of cylindrical 
rods, and it has been suggested that in the liquid the circular cross- 
sections tend to pack in a square pattern, as distinct from the closer 
approximately hexagonal packing in solids. Electron diffraction mensuees 
ments by Lufcy, Palubinskas and Maxwell (1951), support this contention, 
Moore, Gibbs and Eyring, on the basis of the measured densities of liquid 
paraffins and the space requirements of the chains, suggest that for normal 
paraffins with »=5-10 the chains in the liquid at the melting point are 
packed with very little room to spare. For Cz,H,,, Vand (paper to be 
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Paraffin n—C3,H,, crystallized from petroleum ether on formvar substrate. 
Gold-palladium shadowed. x 24 000 


lectron micrograph by H. F. Kay and Mrs. B. J. Appelbe after Dawson and Vand (1951). 
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published in Acta Crystallographica) finds that at the melting point the 
densities are consistent with hexagonal packing in the solid; and with the 
square packing of circular chains in the liquid, distances between chains 
being the same for solid and liquid. Thus the available space in the 
liquid at the melting point is larger for long than for short chains. 

An experiment by Miiller (1940) on two diketones, with n=10 and 
n=11, one with both ketogroups on the same side of the chain, and the 
other with two ketogroups on opposite sides of the chains, show that 
paraffins in the liquid state are very flexible : both liquid diketones have 
nearly equal dielectric constants, which indicates that the two dipoles 
attached to the chains move practically independently even though they 
are only three or four links apart. 


2.6. Mono- and Multi-Layers 


Before concluding the discussion of different states of aggregation of 
long-chain materials, it is worth while to mention briefly mono- and multi- 
layers. Monolayers are, as it were, two-dimensional crystals formed by 
paraffin derivatives with polar endgroups at or near the end of the chains. 
They occur when a small amount of long-chain substance is allowed to 
spread on the surface of water, while the long chains stick outwards. 
Several monolayers may be deposited above one another by a technique 
developed by Blodgett (1935) to form multi-layers, which are similar to 
very thin crystals formed in the usual way (Knott, Schulman and Wells). 
This confirms the similarity of the packing of chains in all states of aggre- 
gation. More recently, interesting work on monolayers was done by 
Stenhagen (1951). 


§3. DisoRDER IN LoNG-CHAIN CRYSTALS 


The crystal structure of a solid represents a solution of the problem how 
to satisfy the bonding forces of atoms or molecules of given shapes and 
forcefields. Usually, one solution of this problem is much better than any 
alternative ones, that is, a certain structure is stable over a fairly wide 
range of temperatures. Long-chain crystals are unusual in that numerous 
alternatives exist of similar stability. These alternatives differ from one 
another partly in minor details of the packing of molecules, and partly in 
the extent to which the bonding forces restrain the twisting and turning 
motions of the long chains. In the following, the experimental evidence 
regarding alternative structures and rotational disorder in the crystals 
will be presented, with particular reference to the latter. Disorder due 
to irregular repeat of the lattice will be discussed later. 


3.1. Huperimental 
Experiments on long-chain crystals employ three main techniques : 
(1) Crystallography, including x-ray and optical methods. 


(2) Calorimetry. 
(3) Dielectric measurements. 
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Of these methods, the first two give information about bulk properties, 
while dielectric measurements are of particular interest because they give 
information on the movements of individual molecules. However, 
while the first two methods are universally applicable, dielectric measure- 
ments yield their information only for polar molecules. Since unsub- 
stituted normal paraffins are non-polar, dielectric measurements can be 
used only for derivatives with polar groups. 

The three methods in question show evidence of rotational disorder for a 
variety of long-chain solids. However, it would not be profitable for the 
purpose of the present review to discuss all the homologous series that have 
been investigated. It will be sufficient to consider in detail only the normal 
long-chain paraffins, since these are theoretically the simplest, and to 
consider derivatives mainly in connection with dielectric properties. 

The calorimetric and x-ray evidence on long-chain paraffins has been 
reviewed by Gray (1943) and by Cines (1950); the latter author is con- 
cerned in the first. place with liquid-solid equilibria. The melting points 
of normal paraffins have been determined up to n=70. For n<18, a 
systematic alternation is found between paraffins with n odd and n even, 
the latter having higher melting points. For longer chains, this alternation 
disappears, and the melting points as a function of chainlengths lie on a 
smooth curve (see fig. 5). Garner and co-workers show that the melting 
points can in general be represented by equations of the type 


Ty=(an+p)\mt8) . 2... - (4) 


(«, 8, y and 6 being constant) which hold over wide ranges of chainlengths. 
This will be discussed later. Formulae of this kind imply a convergence of 
the melting point to a finite value for chains of infinite length. 

As regards solid structure, the normal paraffins may be divided into four 
groups according to chainlengths (Gray 1943) and the different groups are 
not of equal interest here. If we disregard the shortest chains, with n<11, 
three groups can be discerned. For 11<n<18, a systematic difference 
exists for paraffins with n odd and n even. For even chainlengths, the 
forms observed are generally of low symmetry ; for n odd, more symmetri- 
cal forms are found, some of which seem to be hexagonal. However, the 
evidence on this group is comparatively scanty, except for the very 
interesting calorimetric work by Ubbelohde for paraffins with n=14, 15, 
16 and 18, and by Parks and Light for n=14. Lonsdale and Miiller 
determine the structure of C,,H,, with the help of single crystals. This 
paraffin has a triclinic form stable up to one degree below its melting point, 
but changes into a more symmetrical form just before it melts. The 
triclinic form, with the subcell discussed above, occurs also in other paraffins 
with n even. 

While paraffins with »>36 are not very interesting, since they all 
crystallize in the orthorhombic form, those with 18<n<36 show a variety 
of forms and are of the greatest interest as regards rotational disorder. 
This group is most accessible to measurement since it is solid at room 
temperature, and records for it are extensive. 
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The variety of forms observed for the paraffins in question is well 
illustrated by C,,H;, which has recently been examined in detail by 
Mazee (1948). The notation used by this author is unconventional and 
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partly in error and has been rectified in ‘Structure Report’ 11. The 

latter notation is used here, the observed forms being as follows : 

Hexagonal high temperature form: a=4-77, c=33-3, at 46-5°c. The 
melt congeals into this form at 50-7°c, and on cooling the hexagonal 
form is stable down to 42°c. 
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Monoclinic formI: a=7-94, b=4:95, c=33-2, B=94°, at 42°c. This 
form, which is similar to an orthorhombic one arises on cooling the 
hexagonal form to below 42°c. On heating, however, it is retained 
up to 46-5 ©. 

Monoclinic form II: a=7-50, b=4-:99, c=32-7, B=94°. This form arises 
on further cooling of the Monoclinic I form, but it is metastable, 
though long-lived. 

Triclinic form: a—=7-42, b=5-35, c=32-5, a=102°, B=99-5°, y=99°, at 
room temperature. This is the stable form at room temperature, 
and may be formed either via Monoclinic II, or directly via Mono- 
clinic I. On heating, it converts into Monoclinic I at 46°c, but not 
below. 

The behaviour of C,,H;) shows that a number of transitions occur, 
most of which are badly reproducible. The only reasonably well repro- 
ducible transition is that which occurs at 42—46-5 c, from the monoclinic 
to the hexagonal form. Transitions of this kind have been recorded for 
all the paraffins of chainlength between n=20 and 36 that have been 
examined and occur at temperatures of up to about 12° below the melting 
point. The dependence of the transition temperature on chainlength is 
systematic. 

The reproducible transitions in question are found to be connected with 
large latent heats, as shown in table 1, for paraffins of chainlengths between 
20 and 43. Data by two authors are quoted to show the discrepancies 
which are found even for carefully purified materials. These may partly 
be due to traces of impurities, and partly to the ambiguities introduced by 
the ill reproducible transitions between different low temperature forms. 
Data on these can only be inferred and are not well substantiated. 

The evidence of the latent heats of the reversible transitions in question 
shows them to be transitions where the order of the solid state is greatly 
changed. The latent heat of a first order transition 4H is related to the 
corresponding change of entropy 4S by an equation 


AH=TAS* ied a, ee 


where 7’ is the transition temperature. As Ubbelohde and Oldham (1938) 
point out, the entropy of transition can be defined as 


AS=RlogW,/W, ..... . (6) 


where W, and W, are the respective numbers of ways in which the two 
states can be realized. A latent heat of 6 kcal/mole at 50°c corresponds to 
W,/W,=104, a great change in the number of configurations if one bears in 
mind that the provision of two configurations in the high temperature form 
for each one in the low temperature form means only W,/W,=2. 

The calorimetric measurements on paraffins of chainlengths n=20 
to 36 indicate that the transitions observed occur within fairly narrow 
temperature intervals, that is, that they are approximately of the first 
order. However, eqn. (6) is also applicable to a gradual (second order) 
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transition, when AS has to be calculated from specific heat measurements, 
according to 


CoC 4 
AS= | >-"ar LY Cae ae ieee 


where C,, is the measured specific heat and C,, is the specific heat that 
would be measured in absence of a transition. Relevant specific heat 
measurements have been made by Ubbelohde (1938) for a transition in 
C,5;H39. The measurements give a A curve typical of a so-called second 
order transition. More recently, work by Mazee (1948) and Hoffman 
(1951) has cast doubt on the occurrence of gradual transitions in paraffins. 


Table 1 
Total heat of 
fusion of low Heat of 
n Tu T transition | temperature form | transition 
(kcal/mole) (kcal/mole) 
21M 40-3-40-4 32°8 16-9 3:6 
22G 43-8 40-2 18-60 6-90 
23 M 47-3-47-4 40-6 19-4 5:4 
24M 50-7-50-8 47-0 20-5 71 
26 G 56-1 50-1 22-41 8-37 
28M | 61-25-61:3 54-2 22-7 6-1 
30 M | 65-8-65-95 59-2 23-9 6-1 
G 65-5 59-0 25:17 8-72 
31M 67-2-67°3 61-8 25°2 6-7 
| 34 G 72-4 67-9 30°59 11-48 
35 M 74:4-74:5 71-6 28-3 7:3 
36M | 75-8-75-85 73:5 28:1 7:2 
40 M | 81-35-81-45 . 31-5 
43 M | 85-25-85:35 34:3 


Note: M=figures by Mazee (1948); G=figures by Garner, van Bibber 
and King (1931). 


The recent work appears to show that paraffins have sharp first order 
transitions when they are quite pure. Nevertheless, the results of 
Ubbelohde are of interest, since they show that a quite minor disturbance 
of the lattice (Ubbelohde’s material was fairly pure) is able to change a 
sharp transition into a Atransition. In this context it is worth mentioning 
that impurities have appreciable effects not only on the sharpness of 
transitions but also on their temperature. Piper, Chibnall, and co-workers 
(1931) investigate transitions in a number of mixed paraffins, and find that 
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the temperature difference between transition and melting points is 
generally greater for mixtures than for pure substances. 

What is the structural significance of the large change of entropy at the 
transitions in question ? Miiller succeeded in answering this question in 
broad terms, by an ingenious interpretation of the x-ray evidence. If the 
transition in question occurs between forms with chains normal to the 
basal plane, as in CygHgy (and approximately in C,,H;9, since the tilt 
of chains in the monoclinic I form is slight) the structural transition 
corresponds to a change in the basal plane a 6 only. This can be 
described by two parameters, namely the area a b and the angle @ between 
the diagonals of the basal figure (see figs. 2 and 3). 

Miiller investigates the changes of a b and @ with temperature for several 
of the paraffins which show the reversible transitions discussed. He finds 
in some cases sudden jumps at the transition temperature, while in others 
changes are more or less gradual. Figure 6 shows a transition which is 


Angle 6 (obs.) between diagonals of cross-section. 


gradual from beginning to end, for C.;H4,, and this is of particular interest 
for the mechanism of the transition (Mazee claims that the gradualness 
observed for Cy;H,4, is due to impurities, but this is not of major importance 
in the present context). The characteristic feature of the gradual transi- 
tions in question is the behaviour of @, which drops with increasing 
temperature to 60°. The thermal expansion near the transition is aniso- 
tropic and consists of a gradual distortion of the basal figure until it reaches 
hexagonal shape, while the volume expansion is anomalously large (see 
Miller 1930, and Vand’s paper on the density of paraffins (in the press)). 
In many cases this process is sudden, in others it occurs in two stages, and 
sometimes hexagonal character is not quite reached, @ at the melting 
point being higher than 60°. Nevertheless, the behaviour of 6 gives a 
vital clue to the nature of the high temperature forms. 

Miiller points out that the hexagonal form is the closest packing possible 
for straight rods of circular cross-section, while the less symmetrical forms 
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give a more compact packing, possible for the flat cross-section of the 
paraffin molecule in its most stable form. He suggests that the transitions 
in question are due to increased mobility of the chains. If the flexible 
chains twist and vibrate, their cross-section becomes statistically circular 
rather than flat, and they need to be packed in the hexagonal arrangement. 
This explanation fits also the calorimetric evidence, since extensive 
twisting and turning can account for a large change of entropy at the 
transition. 

At this point it is worth drawing attention to the fact that Miiller’s 
elegant interpretation of the crystallographic data makes use of only a 
small part of the information contained in x-ray diffraction patterns. 
A good photograph of a single crystal contains a number of Bragg 
_ reflections, as well as pattern of so-called diffuse or ‘ dynamic ’ reflections 
which are due to the modulation of the periodicity of the lattice by thermal 
vibrations. In principle, such a photograph contains information of the 
movements of molecules which is much more specific than that used by 
Miller. However, this information is difficult to extract, partly because 
of the experimental difficulties of obtaining good single crystals, and partly 
because of the great theoretical difficulty of interpretating diffuse 
reflections. Nevertheless, a very interesting observation has been recorded 
by Lonsdale (1942), who states that diffuse reflections indicate that thermal 
vibrations in long-chain crystals are much more intense in the a b plane 
than normal to it. This confirms Miiller’s ideas. 

While the calorimetric and crystallographic data provide information 
on bulk properties, dielectric measurements are of particular interest 
because they provide information on the mobility of individual polar 
molecules. However, the relationship between dielectric properties and 
molecular mobility is complicated, so that it will be useful to summarize 
the relevant dielectric theory before reporting on the experimental data. 
Tn discussing the theory, it is convenient to treat separately the behaviour 
of dipoles in a static and in an alternating field. The static case will be 
considered first, following Fréhlich’s (1949) treatment in broad outline. 

The dielectric constant of a material measures its polarization in an 
applied electric field, that is, the displacement of the centres of gravity of 
positive and negative charges due to a field. The polarization can be 
divided into four parts : 


P=Pee PsP reP) ae eke te ca (8) 


where P, is due to displacement of electron clouds, P, to that of ions mobile 
with respect to the structure, P, is a contribution due to the distortion of 
molecules and P, is due to the orientation of permanent dipoles. In the 
present paper we are concerned only with the latter contribution, which is of 
interest in that the orientation of dipoles rigidly attached to chains 
implies movements of the chains. 
In long-chain materials, P, is widely independent of time and tempera- 
ture. The non-dipolar part of the dielectric constant can be characterized 
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by «.,, the dielectric constant measured for very high frequencies. This 
is in general not very different from the square of the refractive index, n?. 
However, a difference of the order ¢,.—n?=0-1 or more is observed in 
some cases. 

For a gas of sufficiently low density, containing polar molecules, the 
static dielectric constant «, is given by 

4au,2N 
fo a Re te os 

where N is the number of dipoles per unit volume and p, an average dipole 
moment which for spherical molecules is equal to p. 

For solids the relationship between dielectric constant and dipole 


moment is much more complex than for gases. In general «, can be 
defined 

Se, 47N mm* 
5a, deg eee 
where the quantity mm*, though of the dimension of a dipole moment, 
squared, is different according to the surroundings of the dipoles. 

Frohlich shows that mm* for solids can be calculated if sufficient 
information is available about the structure. The calculation demands a 
knowledge of both short and long range forces acting on a dipole. ~The 
former signify the restraints on the direction of a dipole enforced by the 
forcefields acting in its neighbourhood. They have the effect of reducing 
the contribution of a dipole to P;, or in other words, they reduce mm* in 
the average over all directions. Long range forces signify the effect of 
dipolar interaction on the polarization, and may increase or decrease mm*. 
They depend critically on the configuration of dipoles with regard to one 
another. 

In most solids the dielectric properties are dominated by the short range 
forces, and appreciable dipolar contributions to the dielectric constant 
occur only where these forces permit alternate positions of dipoles. This 
will be briefly discussed with the help of a simple model. 

The positions permitted to a dipole within a solid may be described by 
the potential energy of the dipole as a function of its direction with respect 
to the rest of the structure. In general, it may be a complicated function 
with minima in various directions in space. In the simplest case, however, 
we may assume that the potential energy of a dipole has two minima only, 
when it is pointing in two opposite directions. Fréhlich shows that this 
model is reasonable for a dipole attached to the centre portion of a long 
chain in a crystal; the dipole then has a minimum of potential energy in 
its normal position in the lattice, and we may assume another minimum, 
when the chain has turned round its axis by 180°. 

A dipole with two minima of potential energy in opposite directions can 
contribute to the polarization of the solid of which it forms part only when 
the applied electric field has a component in one of the permitted 
directions. The dipole then contributes to the polarization according to, 


(10) 
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the relative magnitude of the two minima of potential energy. Ifthe two 
minima are of equal depth (fig. 7 (a)) the dipole contributes to the polariz- 
ation according to its dipole moment jy. If the potential energy of the 
two minima differs by a term V (fig. 7 (6)), the contribution of the dipole is 
reduced by a factor 


fb V 

et Oe! opr ee he ee ee LL.) 
(Daniel and Stark 1952). This implies that for VS RT the dipole 
contributes only a negligible amount to the polarization of the solid. 

The model shown in fig. 7 is analogous to that described by eqns. (1) and 
(2), which applies to the rotation of paraffin molecules about single bonds, 
and the reduction of polarizability given by eqn. (11) is analogous to the 
stiffening of chains when rotation is restricted by a difference of potential 
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Schematic model of ‘ double potential well’. 


energy characterized by V or X respectively. Accordingly, the model in 
question is described by the terms of hindered or restricted rotation. 
However, this description is somewhat misleading since no motion 
analogous to the spinning of a top is involved. The term ‘ rotation ’ merely 
implies that the alternate positions accessible to a dipole (or other entity) 
are symmetrical about an axis. 

With regard to the static dielectric constant, only the permitted relative 
directions of dipoles and their potential energies are of interest, but not the 
potential maxima separating the minima. This is not so for alternating 
applied fields, because a potential hill of height W implies a time constant 
7 according to an equation 


Tes CODSts, Cli tie. Wp cpu nse nena we Cle) 
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The significance of 7 is that a large number of dipoles, when disturbed 
in their equilibrium configuration, return to equilibrium with a relaxation 
time r. In an alternating field, the behaviour of dipoles depends on 
whether 7 is large or small compared with the frequency of the field. 

The model used for the discussion of relaxation times is a special one. 
However, the concept of a relaxation time is general, and the dielectric 
behaviour of a material with a given relaxation time 7 and frequency w can 
be described by two equations, derived by Debye : 


Es — Een 
€1(w)—€0= ——s 53 5 Yt Ee ee ae 


and 


(14) 


where ¢, is the real part of the dielectric constant and tan 6 the loss angle. 

Equations (13) and (14) apply to a material with a single relaxation time 
7, but not many materials can be described thus. In most cases, the 
explanation of measured loss angles demands a distribution of relaxation 
times, which implies that different relaxation times apply to different 
groups of dipoles. 

Experimental work on dipolar long-chain materials has produced 
evidence on the static dielectric constant which shows that in some 
cases mm* is negligible, ie. dipoles are immobilized by short range 
forces, while in others mm* is as high as in a liquid or even higher. This 
links up with the crystallographic and calorimetric evidence on ordered and 
disordered forms, and indicates that dipoles are free to move within the 
solid material. The relaxation times give a certain amount of information 
as to the nature of the movements of chains by which the reversal of 
dipoles takes place. 

There exist numerous paraffin derivatives which have a dipole moment, 
but not all of these yield information on the mobility of the paraffin 
chains, as distinct from the mobility of the dipoles themselves. In this 
context, long-chain molecules may be divided into two groups. Firstly, 
molecules where the dipole is rigidly attached to the chain and where it 
cannot perform movements unless the chain itself moves with it. Secondly, 
molecules where a dipolar group is attached so that the dipole can appreci- 
ably alter its direction while the chain remains rigid. To the first group 
belong ketones, halogen substituted paraffins, and esters. The second 
group includes secondary alcohols. Only the first group is of interest 
where the investigation of the mobility of paraffin chains is concerned ; 
the most useful members of this class, as far as theory is concerned, are 
the ketones, since here a large dipole is attached in the plane of the chain, 
at right angles to the chain axis, and the chains themselves are identical 
with those of normal paraffins. 

A great deal of work has been done on the dielectric properties of solid 
dipolar paraffin derivatives, both when they were dissolved in paraffins 
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(Jackson, Sillars, Daniel and Stark) and in the pure state (Buckingham, 
Miiller, and numerous papers by Smyth and co-workers). This work 
establishes beyond doubt that dipoles attached to paraffin chains are in 
many cases free to reverse in the solid state. The dipolar part of the 
dielectric constant is very anisotropic, which confirms the idea that dipoles 
reverse by some twisting or turning motion of the chains and that they are 
able to point in diametrically opposite directions (Daniel and Turner 
1953), 

Dipolar rotation, evidenced by dielectric properties, occurs in some 
long-chain materials within some temperature ranges, but not in others. 
In most cases dielectric and other physical measurements have not been 
correlated, but satisfactory evidence exists to prove that the forms with 
dipolar rotation are at the same time forms with high entropy and volume 
(e.g. Buckingham (1934) and Garner and King (1934) for the ester ethyl 
behenate). This proves that the high entropies in question are connected 
with the turning or twisting of chains. 

The dielectric constant is in some cases as high as in the liquid, in others 
so low that the contribution due to mobile dipoles is negligible, while in still 
others it is intermediate. It may be considered in terms of dipolar rotation 
as will be discussed later. In most cases observed (for stable forms) the 
dielectric constant rises with rising temperature. In some, the rise 
proceeds in several discrete steps, or is gradual, being interrupted by 
steps (Buckingham). The increase of dielectric constant is closely 
correlated with volume, small volume changes being connected with large 
changes of dielectric properties (Daniel and Stark 1952). 

The evidence of static dielectric constants shows that dipole carrying 
chains are sufficiently mobile in the solid state to allow the reversal of their 
dipoles in a plane normal to the chain axes. Dielectric relaxation times 
indicate that the movements responsible for this consist of a twisting, 
rather than a turning of the chains. 

Experiments on ketones and esters dissolved in paraffin waxes of chain- 
lengths larger than that of the polar molecules show that relaxation times 
as function of chain lengths lie on a regular curve. Frdhlich (1942) 
interprets this relationship theoretically, by assuming the chains to be 
elastic rods which are twisted while their axes remain straight. He finds 
a relationship between the relaxation time 7 and the chainlength n given 
b 
a log ant A+ 7 tanh % coceal oad eeatene me ces (LD) 
where A, B and C are constants to be determined by experiment. Figure 
8 shows that this relationship fits the data. The constants give best 
fit for A——50:4, B=13800 and C=26. This formula also gives 
reasonable agreement for materials other than those used to determine Ae 
Band C originally (Meakins 1949, Daniel and Stark 1952). Hoffman 
(1952) gives a slight modification of the formula, which is however, not 
very important. The distribution of relaxation times tends to be rather 
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wide, and in some cases relaxation times were found to cluster about two 
values (Hoffman 1952). 

Before concluding the discussion of the experimental evidence it will be 
useful to draw attention to some very deceptive spurious effects due to 
impurities. These apply particularly to dielectric measurements, and to 
measurements near the melting point of a material. 

Dielectric measurements can be falsified by very small amounts of 
ionized impurities, because one mobile ion introduces a polarization 
(see eqn. (8)) that equals the polarization due to many dipoles. There 
are several mechanisms by which ionized impurities can produce high 
dielectric constants and loss angles, and the spurious nature of the results 
is by no means always easy to detect (Daniel and Turner 1952). 
Impurities are most easily ionized in materials with high dielectric 
constant, and dielectric data for such materials have to be treated with a 
certain amount of caution. 


Fig. 8 
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Dependence of log 7 on chain length. 


Near the melting point not only dielectric but also calorimetric measure- 
ments are subject to errors due to partial melting caused by impurities. 
This occurs if an impurity is soluble in the liquid but insoluble in the solid, 
so that for an appreciable temperature range below the melting point the 
material is two phase consisting of solid, and a liquid solution of the bulk 
substance in the impurity. Calorimetric as well as dielectric measurements 
on such an impure substance near its melting point give results indicating 
a gradual order—disorder transition which culminates at the melting point, 
even though the pure material may have no such transition. 

Three otherwise very interesting experiments by Miiller appear to be 
vitiated by experimental errors of the kind discussed. Dielectric measure- 
ments for the ketone (CsH,,),.CO indicate a gradual rise of the dielectric 
constant towards the melting point, and Miiller (1937) interprets this as 
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an order—disorder transition. He follows this experiment by a further one 
on the same material (1938) where the volume of the specimen under 
test is kept constant. The results are of interest theoretically, since 
conditions of constant volume are more simple than conditions of constant 
pressure. However, later work showed the experimental results to be due 
toimpurities. It should be noted that the impurities in question amounted 
to only 2-3 /, and that their presence was only detected when carefully 
purified material was found to behave differently, the rise of the dielectric 
constant in the solid being the less the purer the material (Daniel and 
Stark 1951). Impurity affects may atfect also Miiller’s interesting 
comparison of two diketones (1940). The rise of the dielectric constant of 
both these substances observed below the melting point is small, and may 
easily be due to impurities. 


_ 3.2. Theoretical Considerations 

The many structural transitions, and the molecular mobility observed in 
paraffins and derivatives appear capricious when viewed on the level of 
morphology. However, the different forms may be viewed also as variants 
of a basic structure whose variability is one of its essential features. 
Attempts to understand the properties of long-chain crystals in general 
have been made at different levels. The most fundamental level is that of 
the Physics of the Solid State, when the solid is considered in terms of the 
properties of individual molecules, and of their interaction. On a less 
fundamental level, a general understanding may be sought in terms of a 
_ phenomenological treatment of molar properties. The simplest, and so far 
most rewarding, approach to the theory of long-chain crystals consists in a 
systematic correlation of the thermodynamic information for different 
materials and structural forms. This approach, which has been initiated 
by Garner and co-workers and owes much to Ubbelohde, makes use of the 
fact that long-chain molecules are members of homologous series and that 
substituted groups of the kind considered here alter the molecular structure 
_ and the bonding forces only slightly. With this in view, the systematic 
appreciation of the thermodynamic data leads to an understanding of the 
variability of the structure of paraffins, and, in favourable circumstances, 
to a method for predicting the thermodynamic stability of ordered or dis- 
ordered forms. 

Table 2 presents a conspectus of the thermodynamic data for normal 
paraffins with 6<n<25. Data for longer chains have already been quoted 
in table 1. 

The quantities listed are defined as follows : 


Pe 
Hy—Hy=| 0, a0 ir homage hee (16) 


that is, H,—H, is the heat content of the solid or gas due to its thermal 
agitation, H, being the heat content at absolute zero. S is the so-called 


‘ third law entropy ’, TO 
s=| yar ee te iaires Tat) 
0 
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the value of S at 7’=0 being zero for any state of aggregation in inner 
equilibrium (Fowler and Guggenheim, p. 223 ff.). 

The data for H,—H, and S for gaseous paraffins of chainlengths up to 
n—8 have been determined with great accuracy, but data for solids are 
less readily available. The figures for S are given in the original papers by 
Parks and co-workers, but those for H,—H, are not evaluated, and all 
data given in the table have been estimated by the author of the present 
review by means of graphical integration of the specific heat data. These 


Table 2. Thermodynamic Data for Normal Paraffins 


gaseous solid 
H-H,| S |H-H,| 8S AH,| 4H, | 4Hz | ASr 
cal/g | cal/g/° | cal/g | cal/g/° | cal/g| cal/g | cal/g | cal/g/° 
for 298-6°K Jat melting point 


Sources 


TTB), 1-072 121 | 35 
TG. SEOLS 121 | 33-8 
76-6 | 0-970 45 0-397 | 130 


- 
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4+HhAw 
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Limit for n—0o 
73:2 | 0-656 
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61 0-450 | 123 
70 133 
120 


eo 
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P=Pitzer 1944; T=Timmermans 1950; U=Ubbelohde 1938: P+H+-T= 
Parks, Huffman. and Thomas 1930; P+-L=Parks and Light 1934: oe 
1948 ; G+B-+K=Garner, van Bibber and King 1931, 


reach down to about 90°K for the measurements of Parks and co-workers 
and only to 200° for those of Ubbelohde. 

The quantities 4H,,, 4H, and 4H, and AS, are the latent heats of 
melting, reproducible transitions, and sublimation respectively and the 
entropies of the transitions quoted. Heats of melting for n>11 refer to 
low temperature forms, the value for high temperature forms being less 


by the heat of transition. For n<11 some figures appear to refer to 
low, others to high temperature forms. 
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Heats of sublimation are more difficult to obtain than heats of melting. 
The figures listed for n=6, 7 and 8 have been obtained by adding the 
heats of melting and vaporization given by Timmermans (1950). 
Ubbelohde’s data for n=14, 15, 16 and 18 are based on vapour pressure 
measurements which are difficult for the compounds investigated. The 
author gives more than one estimate for some of the substances, and the 
figures in the table are arithmetical means. None of the latent heats 
of sublimation quoted have been corrected for differences of specific 
heat between the different states of aggregation. 

In considering the table, we shall first discuss the orders of magnitude 
of the quantities listed, noting that most of these quantities per gram of 
substance reach a limiting value for long chains. The most remarkable 
point is the magnitude of heats of transition as compared with those of 
melting. The heats of transition are appreciable, about 4 to 4 of the 
heat of melting of the low temperature forms and up to 4 of the high 
temperature forms, so that the disordered solid appears almost as a kind 
of missing link between ordered solid and liquid. This condition is typical 
of the so-called rotational transitions discussed by Eucken (1939). 
Eucken’s data show that in some cases, e.g. for O, and for cyclohexane, 
a heat of transition in the solid is larger than the latent heat of melting 
for the disordered solid. 

Apart from the information on transitions, the table is also of interest 
with regard to the statistical thermodynamics of the solid forms. In 
assessing the nature of packing in paraffins, perhaps the most striking 
feature of table 2 is the relative smallness of the heat of sublimation as 
compared with the heat content of the solid at its melting point. To 
appreciate this feature, the heat contents and sublimation energies of 
paraffins are compared with corresponding data for other solids (see 
table 3), estimated from figures given by Parks and co-workers, and in 
the 28th edition of the Handbook of Chemistry and Physics. Since 
specific heat data down to low temperatures and heats of sublimation 
are available only for a few solids, the table cannot claim to be repre- 
sentative of the conditions in all types of solids. However, it seems 
significant that the fraction 7 given in the last column increases with 
chain length, and that for the longest chains investigated it is large not 
only as compared with simple strongly bonded solids, such as a metal or 
ice, but also with a noble gas, and even with benzene, which is also a 
molecular compound bonded by van der Waals forces. 

The significance of the large value of 7 may be appreciated in the light 
of Fowler and Guggenheim’s discussion of the partition functions of 
solids consisting of polyatomic molecules (p. 151). They suggest that for 
a solid containing N molecules of s atoms each, the degrees of freedom 
may be divided into 3N modes of vibration of the solid as a whole, and 
3N(s—1) modes of internal vibrations, in the first approximation unaffect- 
ed by the packing of molecules in the solid. If such modes are indeed 
independent they can be activated without shaking the solid to pieces, 
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so that 7 may be large, the larger the more atoms s the molecules contain, 
However, as a rule one would not expect internal vibrations to be quite 
irrelevant to the order of the solid and so it seems remarkable that in 
solid paraffins H,—H, can reach as high a value as three quarters of the 
lattice energy. This appears to indicate that the internal vibrations of 
chains are independent of one another to an extent which may well be 
uncommon in solids and significant for the theory of long-chain crystals. 

The statistical thermodynamics of solids with complicated molecules 
is not sufficiently advanced to allow more than a qualitative appreciation 
of the order of magnitude of the data given in table 2. However, the 


Table 3 


1 AH” \i Dylon | Hee 
Substance °G cal/g °q cal/g cal/g 4H,+AH,, 


Argon —189 6:7 +185 37-6 14 44-3 0-315 
Mercury 39 2-82 357 65 6:3 67:8 0-093 
Water 0 79-71 100 | 539-55 78 619-2 0-126 
Benzene 5:5 30°3 80 94-3 55 124-6 0-442 
C,H, —56:5 45 130 0-35 
C1 .4H 59 5:5 | 53-3 61 123 0-50 
CoH. 36 78 0-63 
Ore, Do 87 0-71 


68 


33 


4H,,,=Heat of melting. 
AH, =Heat of vaporization of liquid. 
T'm =Melting point. 
T', =Boiling point. The heat of vaporization is measured at 7’y. 


dependence of thermodynamic data or chainlength can be used to derive 
results of quantitative value. 

It is a striking feature of long-chain crystals, apparent also in table 2, 
that for a homologous series the heats and entropies of melting, per gram 
(or per CH, group) are nearly independent of chain length, provided that 
m is larger than 8 or 10. Garner and co-workers show that this circum- 
stance accounts for the character of the melting points in a homologous 
series as a function of chain length. At the melting point the Gibbs 
free energy 

GATS, . ee ites eee 


is equal for the solid and liquid, so that the melting point the latent 
heat, and the entropy of melting are connected by the equation 


H,—H,=4H,,=T ,(S,—S8,)=T 48m. « . (19) 
If 4H,, and AS,, per CH, group is a constant for a crystallographically 
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homologous series in the limit for infinitely long chains, the values for 
finite chains will be 


4H,,=on+B, AS Os ees Aes (20) 


where f and 6 are characteristic for the chain ends. This leads to eqn. (4), 
quoted earlier, for 7',, as a function of n. Garner and co-workers find 
that B, the energy contribution due to chain ends, is negative: the chain 
ends tend to keep the material liquid. 

At this point, it will be useful to digress a little from thermodynamics 
to structural considerations, to explain the alternation of melting points 
observed for paraffins up to n=18 and for odd and even acids. Malkin 
(1931) suggests that this is due to a tendency of these short chains to 
pack with chain inclined to the basal plane. For even chains this is 
favourable, but less so for odd ones and it means a negative contribution 
by the chain ends to the energy H which is numerically larger for odd 
than for even chains: hence the alternation of 4H,, and T',,, while 
AS,, remains constant. The latter means that the freedom of movement 
of the chains is unaffected by the different behaviour of the chain ends. 
For long chains the packing with inclined chains appears to be less 
advantageous, so that the difference between odd and even numbered 
chains disappears. 

Garner and King (1934) propose a tentative theory for the constancy 
of the energy and entropy of melting per CH, group with a homologous 
series. At the melting point, the probability of a molecule leaving the 
crystal must be the same as the probability of a molecule settling on the 
erystal surface. The former depends exponentially on the heat of melting 
according to . 
. aay OxD (a Allele loa tet ven oes, 5 (ZL) 


where y is a frequency, while the latter depends on the probability of a 
chain settling on the surface. This probability depends on the chain- 
length, and on the extent to which chains in the liquid are bent. Garner 
an King find that a probability for a molecule to leave the liquid 


p"=(kln)r™1, 2. 2. ee et ww (22) 


can account for the constancy of latent heat per link, where 7 is the 
probability that a molecule of chainlength n is straight in the liquid. 
The quantity 7 can be calculated from the melting point data. 

~The equations of Garner and King are interesting because the quantities 
involved should in principle be measurable; in particular the shape of 
molecules in the liquid could be assessed with the help of x-ray data, and 
the data on coiling in the gas. However, the assumptions made are 
difficult to check because they involve not only the structure of the 
liquid, but also the kinetics of crystal growth. In principle, since the 
melting point is an equilibrium property of a substance, the constancy 
of energies and entropies of melting per link need not involve kinetic 


assumptions. 
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Garner and co-workers, and later on Ubbelohde and Oldham (1938) 
by their study of several homologous series, find that constant terms 
AH and AS can be ascribed not only to CH, groups within a chain, but 
also to substituted groups, such as the COOH group of an acid or the 
C=O group of a ketone, or to features of bonding, such as a configuration 
of chain ends. That is, the energy and entropy of melting of a given 
substance may be calculated by adding the right number of the character- 
istic contributions of its constituents, and allowing certain terms for 
features of bonding. 

The calculation of heats and entropies of melting from contributions 
of unit parts means a calculation of the structural stability of the solid, 
as compared with the liquid. A similar method can be used also to 
estimate the stability of different forms in the solid. If, in first approxi- 
mation, we suppose that an ordered and a disordered form are the only 
alternatives to be considered, we can again use eqn. (19) which now 
defines a transition temperature 


T= AHAS.:. \-.-eoee ae 


If 7, is lower than the melting point the material will have a disordered 
solid form. If 7’. is above the melting point, the solid material will 
always be ordered. 

The quotient of 4H and AS may be interpreted as characterizing a 
kind of balance between these two terms: if 4H is comparatively large, 
the material is ordered, if 4S is comparatively large the material is 
disordered. This viewpoint leads to an understanding of the structural 
variability of long-chain crystals if we take for granted that in unsub- 
stituted paraffins 7’, is near the melting point. Since 4H and AS consist 
of many small additive terms the change of a single feature of the molecule 
or structure may tip the balance so that of two similar solids one is always 
ordered, while the other becomes disordered at some temperature below 
the melting point. 

The explanation advanced for the variability of structure in long- 
chain crystals is supported by quantitative evidence in a special case, 
investigated by Daniel (1950), where the entropy term AS could be 
varied by an artifice so that the balance between order and disorder could 
be shifted in discrete steps. The materials used in this experiment were 
ketones of equal chainlengths (n=18) but different position of the keto 
group. For the six ketones of this kind for which the ketogroup is more 
than two links distant from the chain ends, all physical properties are 
so similar that these ‘equivalent ketones’ are reminiscent of isotopes. 
This justifies their free energies per mole being taken as equal and for 
an ordered form it may be written as 


Go=H\—TS,.. “et see 


For a disordered form, both H and S are larger, for reasons illustrated 
by fig. 9 which shows the role played by the keto groups in the structure. 
As far as the chains are concerned there are two possible positions for 
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the keto groups in any layer of chains. If they are all on one level, the 
dipoles are arranged so that their interaction contributes a considerable 
term hp to the energy bonding the lattice. On the other hand, if they 
are distributed at random, the entropy of the resulting arrangement is 
larger by a term Rlog 2. If we assume that for the ordered form all 
dipoles are on one level, while for the disordered form they are arranged 
at random the entropy of the disordered form may be written as 
Gg=Hyt+hpt+h,—T(S)+R log 2+s), . . . . (25) 
where h, and s are an energy and entropy corresponding to disorder in 
the lattice due to increase in intermolecular mobility. The balance 
eg the ordered and disordered forms is defined, according to eqn. (25) 
af 
hpth,=T (RB log 2+s). Be apa coe eit bieaineeg 2 3) 
This balance can be shifted by increasing the difference of configurational 
entropy between the ordered and disordered forms by mixing several 
equivalent ketones in equal proportions. This is possible because their 
equivalence implies that such mixing alters the free energy per mole of 


Fig. 9 


Schematic representation of a layer ina crystalline ketone. The positions which 
may be occupied by oxygen atoms are marked by crosses. 


any one ketone only by virtue of a configurational term. If p ketones 
are mixed, the entropy term in eqns. (25) and (26) is increased by R log p, 
while all other terms remain unchanged. This means that for a sufficiently 
large p the disordered form can be made stable at low temperatures. 
Experiment confirms the predictions made on the above assumptions. 
Mixtures of less than the 5 ketones always segregate into the five constituent 
ketones, while for five or more ketones the disordered, random solid 
solution is stable at the melting point. The random form has a hexagonal 
crystal structure and a transition to an orthorhombic form occurs some 
10-15° below the melting point. The measured transition temperature 
and melting points of the different mixtures permit the calculation of 
the latent heat of melting of an equivalent ketone and of hp, the contri- 
bution to the bonding energy due to electrostatic interaction of keto 
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groups. Both values agree with data otherwise obtained, the electro- 
static energy in question being obtained by Frohlich (1946) by direct 
calculation of the basis of the structure. 

The quantitative agreement between the calculated and observed 
shifts of the balance between order and disorder confirms the general 
validity of the concept of such a balance. Besides, the agreement 
suggests a method of quantitatively predicting the stability of ordered 
and disordered forms on the basis of the known energy and entropy 
terms that correspond to individual groups or features. Obviously, such 
a quantitative prediction is only possible if these terms are known fairly 
accurately, together with their dependence on chain length. So far, the 
data available are sufficient only in materials where it is known that 
structural forms of low symmetry do not occur. Thermodynamic data 
of such forms are scanty since transitions between them are ill reproducible. 

Probably the most important consequence of the concept of a delicate 
balance between an entropy and an energy term is that it accounts, to some 
extent, for the variability of structure observed in long-chain crystals. 
The argument outlined suggests that the reproducible transitions of large 
latent heat, observed in some materials, are a feature of long-chain 
crystals in general, and that the differences between materials are more or 
less fortuitous, being due to the shifting of a delicate balance by relatively 
minor structural features. This justifies the attempts made to explain 
the transitions in question in terms of a general theory, on the basis of 
structural considerations. 

The experimental evidence leaves little doubt that the transitions in 
question consist of a change from a state where molecules carry out 
limited vibrations about their position in the lattice, to a state where they 
perform large turns or twists. This means that these transitions are 
analogous to the so-called rotational transitions observed in other materials 
and discussed by Eucken (1939). 

Frohlich (1946) treats a transition of this kind as a co-operative 
phenomenon, an approach suggested in more general terms by Fowler 
(1935) for rotational transitions. Fréhlich approaches the problem 
from the point of view of dielectric theory, and thus applies his argument 
to a dipolar derivative, a ketone. This provides a clear justification for 
using the model of a double potential well discussed earlier (see fig. 7). 
The dipolar ketogroup, rigidly attached to the chain acts, as it were, as 
a pointer which indicates the state of the chain. In the low temperature 
form it has only one stable position, in which it points in the direction 
determined by the structure, while in the high temperature form the dis- 
order of chains is characterized by two opposite directions of equal 
potential energy for the dipole. 

Frohlich’s consideration of a ketone instead of a paraffin has the 
disadvantage that pure ketones show no transition (Daniel and Stark 
1951) and that his results can thus not be verified. Besides it also 
introduces into the theory of long-chain crystals the complications of 
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dipolar long-range forces. Nevertheless, his use of a dipole as a ‘ pointer ” 
allows a very instructive description of a co-operative transition in a 
long-chain crystal. 

In a ketone the keto groups are arranged in planes, in a face centred 
arrangement (see figs. 2 and 3). Frdhlich shows that in such a ‘ dipolar 
plane’ the electrostatic interaction tends to make the dipoles parallel, 
' thus causing a large difference V of potential energy between the right 
(parallel) and wrong (antiparallel) direction of an individual dipole. 
In consequence, the structure is ordered, and a molar parameter can 
conveniently be defined to characterize the order: the degree of order 
being the fraction of ‘right’ less the fraction of ‘wrong’ dipoles. The 
energy difference V is due to the co-operation of neighbours and thus 
depends on the degree of order, and if this decreases, due to thermal 
agitation, V decreases also. This has the effect that the degree of order 
decreases with increasing temperature, first slowly and then catastro- 
phically. The order, and V, vanish completely at a transition tempera- 
ture, which is a simple function of the low temperature value of V. 

Fréhlich’s treatment leads to a gradual transition and an entropy of 
only R log 2, which is much too small to account for the measured data 
even though it refers to constant volume. Besides, his model implies 
ferro electric properties, and thus raises complicated issues not relevant 
to the theory of long-chain solids as such. However, this model is only 
one of many possible ones that can be used for the treatment of a 
transition in a paraffin as a co-operative phenomenon. 

Hoffman (1951) uses a different model, applicable to paraffins and 
derivatives with dipoles at chain ends. He assumes one position of the 
chain of minimum potential energy and a number 2 of subsidiary 
minima, the difference of potential energy between the stable minimum 
and any of the subsidiary ones being V. This leads to an entropy of 
transition of R log 2 and by a suitable choice of 2 the measured entropies 
of transition can be accounted for. This model leads to fair agreement 
with experiment. It is interesting that for a value of 2 larger than about 
6, the transition encountered is of the first order with some gradual change 
before and after the transition, a satisfactory result in view of the experi- 
mental evidence on pure substances. Hoffman also finds agreement 
between experimental dielectric constants and those calculated according 
to his model. Furthermore, he explains cases when materials have two 
relaxation times as being due to energy barriers W between subsidiary 
minima, and a different barrier V-++-W between a subsidiary and the 
stable minimum. Comparison of the dielectric and thermal data leads to 
the conclusion that a value of about Q@=12 can account for all the experi- 
mental data for C,.H,,Br. 

The treatment of a transition as a co-operative phenomenon implies 
that it is possible to define a molar parameter, a degree of order 7, in terms 
of some model which is based on the solid structure. For any such 
model, the co-operative nature of interaction in a solid implies that order, 
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as a function of temperature, drops to zero more or less suddenly. A co- 
operative transition may be described also in terms of thermodynamics, 
a treatment useful for further reference. The Gibbs free energy may 
then be written as a function of r 


G(r) =H (r)—TS(r), os. 2 fe (27) 


and the methods whereby the degree of order 7 as a function of tempera- 
ture may be derived from such an equation are discussed, for instance, 
by Hoffman. 

In spite of the success of Hoffman’s treatment the models used so far 
seem somewhat arbitrary since they bear no close relationship to the 
molecular structure of paraffins. Essentially they treat the molecule 
as a rigid rotator, having no degrees of freedom except for a rotation 
about the long axis. In reality, long chains twist as well as turn, as is 
shown by the success of Fréhlich’s treatment of relaxation times (eqn. 
(15)). In view of this, Hoffman’s assumption of twelve potential minima 
seems artificial, since the high entropy of transition may well be due to 
chain twisting (see Frohlich 1944, 1945) rather than to twelve alternate 
favourable positions of the rigid chain in the forcefield of its neighbours. 

In the last resort, the theoretical treatment of transitions in long-chain 
crystals should be based on the information on molecular structure and 
flexibility which is well established by Pitzer’s work on gases, and on the 
forcefields of the molecules, treated by Miiller. It remains to be seen 
whether this approach may lead to a realistic model of a hindered rotator. 
The problem is difficult, even for the simplest solids discussed by Eucken 
(for instance solid H,, CH, and CD,) and the chain character of the 
paraffins, which proves an advantage for the phenomenological approach, 
is of much less help for the theory on the basis of molecular parameters. 

While transitions in paraffins have been treated only by. the semi- 
phenomenological method of co-operative phenomena, some problems 
of the solid state have been treated in terms of the parameter of molecules. 
These are problems where the interaction of molecules may be neglected. 
Such a neglect seems more justifiable in paraffins than in most other 
solids, since evidence of heat contents quoted earlier, indicate that 
internal vibrations of chains do not greatly interfere with one another. 
A successful example of a treatment of the kind in question, namely — 
Fréhlich’s calculation of relaxation times, has been quoted earlier 
(eqn. (15)). This calculation is based on the simple model of a rod or 
ribbon of a certain rigidity, and Frdhlich (1945) uses this model also to 
calculate the dielectric constant of long chains with several polar groups 
attached at intervals. The latter calculation has so far not been verified 
experimentally. 

Szigeti (1952) uses a model closely related to that of Pitzer (1940) to 
investigate theoretically the modes of internal vibrations of paraffins 
and ketones. He assumes that the chains are in their straight, flat, 
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configuration and neglects, in the first place, the interaction with neigh- 
bours which keep them so. However, the author id 

: : promises to consider 
these interactions in a later paper. The published papers contain the 
interesting result that the amplitudes of torsional vibrations of the 
paraffin chain are appreciably influenced by the end groups even for long 
chains. For ketones Szigeti calculates a comparatively large contribution 
of the torsional vibrations to the dipolar part of the dielectric constant. 

Szigeti’s work represents an approach to a rigorous treatment of the 
statistical thermodynamics of long-chain crystals, and is thus of great 
interest. So far, such a treatment is so much in its infancy that it may be 
worth while to discuss briefly some of the problems that await a solution. 

Many problems of the disordered form seem to revolve around the 
mutual interaction of the flexible molecules. For instance, it may be 
asked : Which of the twisted configurations occurring in the gas phase are 
suppressed in the solid, ordered and disordered ? In other words : How far 
is the increased entropy of disordered forms due to twisting, and how far 
to turning of the chains ? Furthermore: By what means do chains inter- 
lock to give the mutual interaction that leads to the co-operative nature 
of transitions ? Do neighbouring chains rotate co-operatively in clusters ? 
So far, none of these questions can be answered. 

Apart from the questions mentioned above, there remains unsolved also 
the most important problem regarding long-chain crystals, namely the 
structural reason for their variability. Experiment shows that for a 
given substance, a wide variety of structures have very similar free energy 
at room temperature. This means that when two structures differ in 
entropy by a term AS, the energy of thermal vibrations is increased by 
a term 7'4S8, where 7 is always near the melting point. This is valid 
not only for the rotational transitions discussed above, but also for the 
changes of structure involving forms of low symmetry. This feature, 
which makes the morphology of long-chain crystals so very complicated, 
must in the last resort be explicable by the way in which the packing in 
the solid restrains the different movements of the flexible chains. 


3.3. The Problem of Premelting 


The similarity of the packing of paraffin chains in the solid and liquid 
states, and the molecular mobility in disordered forms, tend to convey 
the impression that for long-chain crystals the difference between solid 
and liquid is less than it is in general, for pure substances. This, essentially 
qualitative, argument suggests that a study of the disordering and 
melting of long-chain crystals is of interest for the understanding of the 
liquid state and of melting, and has led to investigations of ‘ premelting ; 
in paraffins, particularly by Ubbelohde. However, this concept needs to 
be clearly defined if it is not to be misleading. 

An exact definition of premelting is not practicable because of the 
lack of an accepted theory of the liquid state and of melting. (Problems 
relating to this subject have been recently discussed in the Royal 
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Society, in March 1952.) For the purpose of a discussion of relevant 
concepts and variables, it seems convenient to start from the theory of 
Lennard-Jones and Devonshire (1939). With this theory, melting, at 
any rate for spherical atoms, is defined as a co-operative increase of the 
number of flaws (holes and interstitial atoms) in the solid, culminating 
at the melting point in a break up of long-range order. This implies a 
clear distinction between a rotational transition, which affects only the 
freedom of molecules to rotate (or twist) at their appointed lattice sites, 
and between melting, which affects the number of flaws which disturb 
the long-range order. A rotational transition need not imply a marked 
increase of flaws, and may thus differ in kind from melting. 

In general any crystal contains in equilibrium some holes and/or 
molecules in interstitial positions, but the number of these is very small. 
Fowler and Guggenheim (p. 542 ff.) give approximate formulae for the 
average number NV” of holes and NW‘ of interstitial atoms for a crystal 
of one component. These are 


N*=N . exp (—CYRT) 4 ee eee 
Ni=aN .exp(—CYET), . . a} 


where « is a constant of the order of unity and V the number of molecules. 
¢” is the energy to remove a molecule from an internal lattice site to the 
surface, and ¢' the energy to take it from a surface point to an inter- 
lattice site. 

While the energies €* and ¢' are not very accessible, and besides, 
somewhat ambiguous, since not all.surface sites are equivalent, it is yet 
clear that the number of holes in equilibrium in a crystalline paraffin 
will be very small if ¢” is as large as even a half or a quarter of the subli- 
mation energy. Thus a large number of holes is to be expected only if 
eqns. (28) and (29) are not valid, e.g. because the number of holes is 
increased by a co-operative process. 

In general, pure solids melt very sharply, instead of more or less 
gradually as predicted by the theory of Lennard-Jones and Devonshire. 
It would thus be very interesting if long-chain crystals were found to 
be an exception. This seems plausible to some extent, because the 
difference between the packing of chains in the solid and liquid states may 
be reduced largely to a difference in their packing in the basal plane since 
the chains pack more or less as straight rods in the liquid as well as the 
solid. Vand, in a paper on the density of long-chain compounds (to be 
published in Acta Crystallographica) suggests that chains in the 
liquid form a square mesh, while in the solid their packing is approxi- 
mately hexagonal. If cross sections are regarded as gears that mesh, 
a square array is not resistant to shear, and is thus consistent with a 
liquid structure, while in the hexagonal array the molecule in the centre 
of the hexagon causes a resistance to shear. This theory may not apply 
to the short chains discussed by Moore, Gibbs and Eyring, but in any case 
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melting seems to consist largely of a re-arrangement in the basal plane, — 
and might be more gradual than the co-operative. increase of holes in a 
three dimensional array. 

Oldham and Ubbelohde (1940) present experimental evidence which 
suggests that near the melting point of paraffins the number of holes in 
equilibrium in the lattice is abnormally large and that it may increase 
progressively with rising temperature. Their argument for premelting 
is based on three kinds of experimental evidence: measurements of 
specific heat (Ubbelohde 1938), cryoscopic measurements near the melting 
point, and experiments with holes introduced into the lattice by means of 
mixing. 

The significance of specific heat measurements near the melting point 
depends, as has been discussed earlier, on the purity of the material used. 
Ubbelohde (1938) claims that impurities should be negligible for his 
measurements for C,,H3, which shows a gradual increase of C,, as the 
melting point is approached. However, as far as the question of holes 
is concerned, this evidence is not quite relevant, since the increase of 
C’, may be due to the beginnings of a rotational transition. 

The cryoscopic method used by Oldham and Ubbelohde (1940) is 
based on the melting point depression of a substance due to admixture of 
a second substance that is completely soluble in the liquid and completely 
insoluble in the solid. The melting point depression gives the latent 
heat of melting, and Oldham and Ubbelohde find that for C,,H,, and 
O,H,,.CO.C,H,, this decreases, as the melting point is approached, by 
9% and 34% respectively. This is an interesting confirmation of the 
results of the specific heat measurements, because for the cryoscopic 
Measurements impurities act so as to make the latent heat of melting 
spuriously too large near the melting point. However, the decrease of 
H near the melting point may still be due to the onset of a rotational 
transition. 

The mixing experiments of the same authors give more direct evidence 
on holes. The character of the holes in question is shown in fig. 9, a 
schematic representation of a mixture of the ketones C,H, .CO.C,,H 9, 
and C,;H,.CO.C,;H,,, the holes in question corresponding to missing 
CH, groups rather than whole chains. (Incidentally the number of holes 
is likely to be less than that shown, because some of the layers may be 
upside down, so that their ragged ends interlace. However, this does not 
invalidate the argument.) 

Oldham and Ubbelohde (1940) on the basis of experiments on defect 
lattices of different kinds, and the other data, suggest that near the 
melting point the number of holes increases co-operatively, ie. that 
genuine premelting takes place. However, the evidence is not conclusive. 

In investigating premelting in the sense of agglomeration of holes near 
the melting point, the existence of the order—disorder due to twisting 
is rather an obstacle, not only because the two kinds of phenomena are 
difficult to separate, but also because they must of necessity influence one 
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The existence of holes reduces the bonding energy which restrains 
and thus reduces the temperature at which a form disordered by 


another. 


having a lower energy ¢”. 


s the most direct test of premelting would be the accurate 
number of holes increases co- 


, both by x-rays and with a method 
the x-ray and macroscopic densities 


Conversely, it must be easier to create holes in 


Perhap 
measurement of density for a paraffin 


a structure disordered by twisting and thus 
giving the macroscopic density. If the 
operatively near the melting point, 


twisting becomes stable. 


twisting, 


Fig. 10 
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3) in CygH,,O (m 


(m= 


Composite lattice with random holes obtained by packing C,,H;,0 


This method would 


would diverge progressively with rising temperature. 


distinguish between hole-formation and a rotational transition, even 
though it would not eliminate the difficulties due to impurities. 


3.4. Irregular Structures and Unstable Forms 
Long-chain crystals may contain a great variety of irregularities of 


different kinds, of which the holes, discussed earlier, are a special case. 
Irregular structures may be in equilibrium, or they may be 


Such irregularities may be present in crystals of pure substances, but 
even more so in mixed crystals, and they may be present in large concen- 


trations. 
unstable, but so long lived that they may, in somewhat loose usage of 
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the word, be called metastable. The long life of unstable structures is 
also characteristic for structures which are not irregular. 

It is not possible in this review to discuss all the irregular and metastable 
structures investigated, and a selection has to be made of structures 
where the irregularity can easily be interpreted, or where physical 
properties are interesting. 

A simple kind of irregularity has been studied for ketones and esters. 
If a substituted dipolar group is not in the centre of the chain to which it is 
attached, an ordered crystal structure implies that within a layer all the 
groups should be on one level (see fig. 9). However, irregularities may 
arise if some chains are ‘ upside down’ as compared with the others. If 
this is the case for many chains the structure may be described as one where 
substituted groups are distributed at random between two levels in a layer. 
For this type of irregularity, the number of flaws can be calculated fairly 
simply as a function of the thermodynamic data of the structure. For a 
pure ketone, Oldham and Ubbelohde (1940) point out that the number of 
upside down chains is approximately 


N,=N exp (—A,/kT) . . . . . . (30) 


where h, is the contribution of the ketogroup to the heat of sublimation. 
They estimate V,/N=2x10~4, a small number. For mixed equivalent 
ketones of equal chainlength and different position of the keto group 
Daniel (1950) shows that the number of upside down chains increases with 
the number p of ketones mixed, until a random form becomes stable for 
p=5. Random forms are unstable for lower values of p, being the more 
long lived the larger » is, i.e. the less the difference of free energy between 
the ordered and random forms. 

Experiments by Crowe, Hoffman and Smyth (1951) for esters, show a 
similar effect, except that a random form is long lived for the pure 
materials. This can be explained by the fact that esters have a smaller 
dipole moment than ketones, so that the regular arrangement of ester 
groups contributes less to the bonding energy than is the case for ketones. 

The irregularities due to upside down chains are probably the best 
clarified example of structure with irregularities, but many other such 
structures exist with interesting properties. Such are, for instance, the 
mixtures of paraffins of different chainlength investigated by Piper, 
Chibnall, Malkin and others (1931), and the mixtures of ketones and para- 
ffins investigated by Oldham and Ubbelohde (1940). Oldham and 
Ubbelohde point out that the co-operative nature of ordered lattices implies 
that the number of flaws is usually either large or small, intermediate 
numbers being restricted to narrow ranges of concentration. Irregular 
structures exist also for pure compounds where no ‘ upside down’ chains 
occur, as found by Vand, Lomer and Lang for potassium caprate. 

One of the most interesting examples of unstable long-chain crystals are 
the alcohols investigated by Meakins, Sack and co-workers in a series of 
papers. The static dielectric constant of secondary long-chain alcohols, 
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when freshly cooled from the melt, was found by Meakins and Walsh 
(1951) to reach values up to «=15. This value is considerably higher 
than can be accounted for on the basis of dipole reversal, as observed in 
ketones, since the dipole moment of the COH group is only about 
1-6 10-38 e.s.u. as against »~=2-8 x 10-18 e.s.u. for the CO group. The 
measurements of loss angle show a wide distribution of relaxation times 
clustering around 7=10~4 seconds, and this value does not significantly 
vary with chainlength. However, according to Dryden the relaxation 
time decreases with rise of temperature. Ageing of the high dielectric 
constant forms decreases the dielectric constant while the relaxation time 
does not change significantly. 

Meakins and Sack (1951) explain the high dielectric constants found as 
being due to the formation of chains of linked hydrogen bonds, in analogy 
to Ubbelohde and Woodward’s (1946) theory of the ferroelectric properties 
of Rochelle salt. The model envisaged implies that the polarity of a long 
chain of linked hydrogen bonds can be reversed when each proton in the 
chain transfers its allegiance from the oxygen on its left to the oxygen on 
its right, or vice versa. Sack (1952) treats this model quantitatively, and 
derives both dielectric constant and relaxation time as functions of the 
number of linked hydrogen bonds. The polarizability should increase 
with the cube of the number of linked bonds, and the relaxation time with 
the square of this number. Dryden (1952) in a later paper, quotes 
calculations by Sack according to which the relaxation time for very 
large numbers of links reaches a saturation value. Sack’s theory is 
able to account qualitatively for the observed dielectric properties, but a 
quantitative verification, which would demand an independent measure- 
ment of the number of linked bonds, has so far not been possible. 

The high dielectric constant forms are unstable, the stable forms having 
low dielectric constants. Experiments by Meakins and Sack indicate that 
the difference of stability between the two forms is less at high than at 
low temperatures. The authors tentatively suggest that the formation 
of chains of linked hydrogen bonds demands a crystal structure which is 
unfavourable with regard to the packing of paraffin chains and that these 
two competing factors account for the difference between the high and 
low dielectric constant forms. At high temperatures the hydrogen bonding 
predominates, while at low temperature a structure is formed similar to 
that of paraffins, the long chains of hydrogen bonds being broken up. 
This explanation seems difficult to reconcile with the evidence of the 
melting points. These are higher than those of ketones of equal chain- 
length, and this implies high bonding forces. Such forces can be explained 
without difficulty on the basis of hydrogen bonding, which contributes 
some 5 kcal/mole to the bonding energy. However, if we assume that in 
the low dielectric constant form the hydrogen bonds are broken up, it 
seems difficult to explain its melting point, which is even somewhat higher 
than that of the unstable high dielectric constant form. Further work is 
needed on this problem. 
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In view of the prevalence of long-lived unstable forms, and of the 
importance of thermal movements in paraffins, it seems curious that no 
systematic work had been done on diffusion in long-chain crystals. In 
principle, the existence of long lived unstable forms may be due to two 
causes : molecules may change place only rarely, or the difference of free 
energy between stable and unstable forms may be small. 

The number of place changes of molecules in a solid is characterized by 
the coefficient of diffusion D, and this obeys in general an equation 


D= A exp (—E/kT) 3 Se Sopeeteren Ge 


where A and E are approximately constant. The energy H represents the 
barrier to be surmounted when two molecules change place. For self- 
diffusion, this is characteristic for one kind of molecule alone, and of 
course, for the crystal structure. The magnitude of # and its dependence 
on chainlengths would be of interest with regard to the formation of holes 
in the lattice. Besides, it would be interesting to compare H with the 
activation energy of dipole reversal. 

Although no work has been done so far to measure coefficient of diffusion 
in long-chain crystals, or to assess the role of thermodynamic conditions 
with regard to the long life of unstable forms, recent experiments by 
Daniel and Turner (1953) suggest possibilities of approaching this problem. 
Mixtures of p equivalent ketones which in equilibrium consist of p separate 
phases, can be retained for a time as unstable solid solutions which 
subsequently break up. This involves the development of thermal 
fluctuations into nuclei of the p phases, and subsequent growth of the 
nuclei, by diffusion, a process which can be followed at or near room 
temperature. Conditions of this kind have been used by Daniel (1948) 
to measure coefficients of diffusion in an alloy. If similar work, which 
demands difficult x-ray techniques, were to prove possible for equivalent 
ketones or other equivalent paraffin derivatives it would be of interest, 
since ‘equivalent’ substances may be assumed to have very nearly 
equal coefficients of diffusion, so that the coefficient measured would be 
practically a coefficient of self-diffusion. Besides, the differences of thermo- 
dynamic stability in such systems can be assessed fairly accurately and 
can be varied in small steps. 

Another method of measuring coefficients of diffusion, using multilayers, 
has been suggested by Professor Piper in a private communication. Multi- 
layers of different substances may be deposited above one another, and 
their interdiffusion might be followed with the help of x-rays. A method 
of this kind would yield coefficients of diffusion in the direction of the 
long-chain axis, and thus provide data on the anisotropy of diffusion. 


3.5. Long-Chain Plastics 


From the practical point of view the greatest interest of long-chain 
crystals lies in the help they may give to the understanding of plastics. 
Numerous important plastics are very closely related to paraffins. In 
particular, polythene is a mixture of normal long-chain paraffins of 
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chain lengths of the order n= 1000. Experiments shows that such plastics 
show order—disorder transitions, and that admixture of other materials, 
so called plasticizers, affects order and disorder. The rather complicated 
conditions in question are analysed by Boyer (1949) who finds conditions 
rather similar to those in long-chain crystals : for instance, dipolar groups 
which bond by dipolar interaction, increase the transition temperature 
for order to disorder. 

The connection between disordering in plastics and crystals can most 
conveniently be stated in terms of thermodynamics. Firth and Tuckett 
(1944) explain why polythene changes gradually from order to disorder, by 
assuming that the material consists of a fraction x of ordered regions 
where the chains are aligned as in crystals, and of a fraction 1—a of 
amorphous regions where the chain ends become entangled. ‘They assume 
that the entropy of disordered regions depends on a because entanglements 
are the more confused the more chains take partinthem. This leads to an 
equation for the free energy which may be written 

G= (1 x ) Gea ceases 1s TS (x) 
where the first term applies to the average of crystalline regions, the 
second to the average of amorphous regions, and S(x) is characteristic of 
entanglements. 

The free energy equation for the plastic is more complicated than that 
for the crystalline material, because of the inhomogeneity of order. 
However, there seems no reason why some of the individual terms in the 
equation should not be derived from a knowledge of long-chain crystals. 
The experience obtained on such crystals shows that it is possible to 
subdivide energies and entropies into contributions from individual mole- 
cular groups or structural configuration, and that these contributions 
remain constant for numerous different materials. There seems to be no 
reason why the same should not apply to long-chain plastics in thermo- 
dynamic equilibrium. It would thus seem that the disordering in plastics 
such as polythene is due to similar causes as in long-chain crystals, but that 
in plastics additional causes of disorder are present, due to entanglements 
and differences in the length and constitution of the constituent chains. 

So far the systematic consideration of structure and bonding of long- 
chain materials has been mainly concerned with paraffin chains. However, 
recently two new classes of long-chain plastics have gained prominence, 
namely the silicones and fluorocarbons. The fundamental work on long- 
chain crystals may have practical applications for these new materials. 


$4. CONCLUSIONS 
The behaviour of long-chain crystals may be interpreted on different 
levels, in terms of different concepts and variables. On the level of 
morphology, the many structural transitions observed can only be 
catalogued. At the other extreme, on the level of theoretical physics, 
the structural variability of long-chain crystals is a problem of the 
statistical thermodynamics of solids consisting of polyatomic molecules. 
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a far the most fruitful interpretation of the experimental facts can be 
eG ivarrcact i duo be havnt 
ct that long-chain 
molecules consist of repeating unit parts, and that constant energies and 
ay terms can be ascribed to these units. The free energy of a given 
g-chain material and structure can be assembled by adding appropriate 
numbers of such terms, and this means that the relative stability of different 
structures can be estimated. This method is not very fundamental 
because it does not by itself imply a knowledge of the structural significance 
of the different unit terms. However, its application shows that order and 
disorder in long-chain crystals are delicately balanced, so that the crystal 
structure of a long-chain material may be decided by a single constituent 
group of its molecules. 

On a more fundamental level, the basic principles of the crystallography 
and the bonding forces in long-chain crystals are well enough understood. 
Paraffin chains are bonded by van der Waals forces, and these suffice to 
keep the chain axes straight without implying localized linkages from 
molecule to molecule. The packing approximates that of straight rods, 
so that the mutual arrangement of molecules reduces, in first approxi- 
mation, to a two dimensional problem, the packing of the cross sections of 
chains in a plane normal to them. 

The basic principles outlined are not specific enough to prevent the 
structure of every long-chain crystal from being an individual, and not an 
easy problem. However, from the point of view of the physicist, the 
problems of greatest interest are those which concern long-chain crystals as 
a class. Unsolved problems of this kind concern mainly the thermal 
vibrations which seem to contain the clue to much of the variability of 
long-chain crystals in general. 

The thermal vibrations of long-chain crystals appear to be largely 
internal vibrations, which do not disturb the arrangement of chain axes, 
so that the heat content of a crystal may reach three quarters of its heat 
of sublimation. In low temperature forms, the thermal vibrations are 
mainly limited torsional vibrations, but in high temperature forms the 
chains become twisted and/or turn around. It has been shown that 
transitions between the two kinds of forms are co-operative, but not enough 
is known about the extent to which the vibrations in high temperature 
forms include rotations of chains as a whole, and about the way in which 
movements of neighbouring molecules interact with one another. The 
mechanism of melting also involves interesting problems. 

Apart from general problems concerning the properties of a typical 
long-chain crystal, a great many problems remain regarding individual 
materials, and it should be borne in mind that comparatively little is 
known about the sluggish transitions between the various low-temperature 
forms. Apart from these, certain disordered and unstable structures 
present interesting problems, in particular where dielectric properties are 
unusual. 

DL2 
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- The value of a detailed theoretical investigation of a class of materials 
such as the crystalline paraffins and derivatives depends not only on its 
intrinsic interest, but also on the practical use to which the knowledge 
obtained can be put. The long-chain molecules in question are such that 
chemical features can be modified by mixing, so that solid structures, 
ordered or disordered, can in principle be tailormade to practical specific- 
ations. In this regard the long-chain crystals, as a class, are analogous to 
alloys or ionic semiconductors. Possibilities of this kind are valuable where 
materials have useful properties. So far, the long-chain crystals as such 
have no practical application, but their close relationship to important 
types of plastics gives to their theory a possible practical significance. 


Thanks are due to Professor Fréhlich and Mr. C. G. Garton for their 
interest and helpful criticism, and to Dr. Vand for valuable discussions 
and the communication of unpublished material. The article is published 
by permission of the Electrical Research Association. 
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§ 1. INTRODUCTION 


THE interplanetary dust swept up by the earth as it moves around the sun 
is divided into two categories, the sporadic and shower meteors, res- 
pectively. The former are distributed over a large volume of space and 
have a wide range of velocities and directions of motion. A meteor 
shower consists of a homogeneous velocity group of particles moving in 
a stream with a well defined orbit, so that the particles incident on the 
earth are moving. in approximately parallel paths and appear to an 
observer to radiate from a fixed position on the celestial sphere known as 
the radiant point of the shower. 

Tn recent years radio echo techniques have been applied intensively to 
the study of the ionization created by meteors in the region of the atmos- 
phere 80-120 km above the earth’s surface. Much of this work falls 
within the field of Meteor Astronomy in which important advances have 
been made concerning the delineation of the orbits of streams, the 
distribution of velocities of sporadic meteors (thereby establishing them 
as members of the Solar system), and the discovery of the daytime 
meteor streams, to mention but a few. It is not the purpose of this 
report to discuss these questions, but to present an account of current 
investigations, utilizing the radio echo methods, in the field of Meteor 
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Physics which is concerned with the physical properties of the meteors 
and the columns of ionization which they produce, and of the atmosphere 
in the meteor ionization region. 

The equipments used cover roughly the wavelength range 4-10 m and 
may radiate energy continuously or in short pulses. Except where 
otherwise indicated, the experiments referred to were made with pulsed 
equipments. 

The units in the following, unless specifically defined, are c.g.s. 


§2. Tue THrory or Mereor IoNnrIZATION 


A meteor, passing through the upper atmosphere, is subject to collisions 
with air molecules, a great proportion of which are trapped in its surface. 
The heat liberated causes the meteor to evaporate atoms which carry the 
meteor velocity relative to the surrounding atmosphere and subsequent 
collisions between these and air molecules are responsible for the luminosity 
and ionization. 

The modern theory describing the motion of a meteor in a rarified atmo- 
sphere (Whipple 1943, 1952; Herlofson 1948) is based on the following 
assumptions :—(a) That the meteoric particle should make individual 
impacts with air molecules. This requires the meteor dimensions to be 
small compared with the molecular mean free path in the atmosphere. 
We will see later (§11) that this will be so for all except the brightest 
visual meteors. (b) That energy loss due to direct radiation from the 
meteor should be negligible. According to Whipple (1952), particles with 
radii less than about 10~4cm (which he terms ‘ micrometeorites ’) will 
dissipate most of their heat energy by radiation and hence will not cause 
ionization. This limiting size is far below the range observed both 
telescopically and by radio techniques. 

Following the above authors we can derive the following expression for 
the rate of evaporation of meteor atoms : 

ee (uH) mv COS x (2/Pmax)(1—4-P/Pmax)” 3 iS ales (2.1) 
Pmax 18 the atmospheric pressure at the point of maximum evaporation and 
is given by 
Pmax=[(2lg)/(Av?A)]m" cos x. eRe) 
The quantities appearing in (2.1) and (2.2) are as follows: »=number of 
meteor atoms evaporated per second ; v=meteor velocity (geocentric) ; 
m=initial mass of the meteor ; «mass of an individual meteor atom ’ 
/=latent heat of evaporation of the meteor ; y=zenith angle of the meteor 
path ; p=atmospheric pressure ; H7=atmospheric scale height ; g=gravi- 
tational acceleration ; A—heat transfer coefficient defined as the ratio 
(heat energy liberated)/(air mass intercepted by the meteor x $v") ; 
A=shape factor defined by A=.%/m-*/3 where . is the cross sectional area 
normal to the line of flight. 

Herlofson takes A=1, i.e., he assumes that all impinging air molecules 

are trapped in the meteor surface. while Whipple suggests a value of 0-7. 


Radio Echo Studies of Meteor Ionization 497 


It should be noted that A is independent of m, being a function only of 
the shape of the meteor and its density, p. For a spherical particle, 
A=1-21p-2/, 

In deriving (2.1) and (2.2) the meteor deceleration is assumed to be small, 
and approximations have been made which require v?> 121. 

From (2.1) and (2.2) we obtain the following expressions for the electron 
line density, «, produced by the meteor : 


ee nas (Dice) ut Di Donan) os pete aS a (2.3) 
pa=e SUE \e? WM COBY. ae ek te eh (2A) 
h ob 2.2 2 
re ce ( ) Pmax= ue (Haas cos? Neo Oh Gp ore 5 (2.5) 
where L /p\13 
o= 7a (5) eS 


@ is a function of the properties of the meteor and of 8 which is defined as 
the probability that a single evaporated meteor atom will produce a free 
electron. 

The ionization curve derived from (2.3) is given in fig. 1, with abscissa 
—In (p/Pmax)=(h—Inax)/H, where h is the height of the meteor in the 
atmosphere and h,,,, refers to the point of maximum ionization. The 
maximum ionization produced is proportional to the initial mass and may 
depend on velocity through the probability of ionization, 6. ‘The meteor 
luminosity likewise follows a law of the same form as (2.3). 


Fig. 1 


~ln (p Perry) = (h-fmax)/ H 
~2 = ro) | 2 


3 4 
Ionization along a meteor trail as a function of height. 


§ 3. Disstpation OF METEOR IONIZATION 


We will consider here various processes which cause the decay of the 
volume density of meteor ionization. 


(a) Diffusion 


The most important of these is ambipolar diffusion which reduces the 
volume density without affecting the total electron line density, «. 
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A mathematical treatment has been given by Huxley (1952), which is 
essentially similar to the following. 

The diffusion equations for a cylindrical distribution of electrons and 
positive ions in the absence of a magnetic field are given by (3.1) and 


(3.2) respectively, 
D) sa # (re. En, 
Re ome) +Y Gran) 2 ee 


On; D; 0 On; eEn; 
Fe ES (rE) a—\/ (=) ° ° . ° (3.2) 


where 7=radius, t=time, n, and n, are the electronic and ionic concentra- 
tions, D, and D, are the respective diffusion coefficients given by the 
kinetic theory, m and m, are the electronic and ionic masses, v, and v; the 
collision frequencies, e the electronic charge and E is the space charge 
electric field ; eE/(m,v,;) and —eE/(my,) are the drift velocities of ions and 
electrons respectively in this field. 

Taking advantage of the fact that throughout the main part of the 
cylindrical column the difference between n, (and its derivatives) and 7; 
(and its derivatives) will be small, we can eliminate E from (3.1) and (3.2) 


and obtain 
on,” DO. / 50K, 

eC x). ow «bleh ueers (3.3) 
where the ambipolar diffusion coefficient, D, is given by 


my, mp, 
she (Det mo, »,) | Cee) 


From the kinetic theory, 


and 


where A, ;=mean free path, 7’, and 7’; are the electron and ion tempera- 
tures. Hence 
D=D, (+7, /7 )=2D; , (3.4) 
if'we put 7.=7 =—T. 
Using the kinetic theory expression for D;, assuming the ions and neutral 
molecules to be of approximately equal mass, j., we obtain 
8 /2 1/2 
D= = A; anal 
3 *\ ap 
A= (/27nQo 7)-} 
where k is Boltzman’s constant, n, is the atmospheric particle density and 
o;, is the mean classical diameter of the ions and neutral molecules. 
The initial radius of the meteor trail will be of the order of the molecular 
mean free path A,, which in the meteor height range is a few em. Since 
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we will be mainly concerned with the behaviour of a meteor trail when it 
has expanded to a radius considerably larger than this, we may assume an 
initial line distribution of ionization, and solution of (3.3) gives 


Me= TOP ROS Oa eae eG) 


Thus the radial electron density distribution is Gaussian and the 
quantity (4Dt)!/2 may be taken as a measure of the radius of the column 
at time t. 

The effect of the earth’s magnetic field will be to reduce D,, and hence 
the space charge electric field, E, in the direction transverse to the field ; 
its direct effect on the ions in the meteor ionization region will be negligible. 
We would thus expect the trail to become elliptical in cross section. At 
some height the coefficient, D,, for transverse diffusion will become equal 
to D,, leading to D=D, (i.e. a reduction by a factor of two), and at greater 
heights the space charge field will be reversed and we will get D<D,, i.e. 
the slowly diffusing electrons will retard the positive ions. Now for 
transverse diffusion, the diffusion coefficient, D,, is reduced in the ratio 
v,7/(v.2+w,") where w, is the angular velocity of the gyrating electrons. 
The condition D,=D,; is satisfied when (w,/v,)?=4(2m,/m)!2 which, 
putting w,~10? sec—1, m,/m~6 x 104, becomes v,~3 x 10°. The collision 
frequency will be of this order for a neutral particle density of 210 
which corresponds to a height of about 100 km. Thus, below this height, 
the earth’s magnetic field should have little effect. 

Owing to diffusion, the shape of the column of ionization will be that 
of a paraboloid of revolution extending back from the meteor. However, 
for many purposes it is sufficient to assume that the meteor instantaneously 
creates a narrow column which subsequently diffuses radially. 

The energy released by the meteor will cause heating of the column of 
air and may initially enhance the rate of the expansion. From (2.1) and 
(2.4) we find that the heat energy released by a meteor producing « 
electrons per cm is dyvx«/B erg cm—!. For a meteor with v=40 km sec}, 
B~1 (see § 11), <=10" cm-!, we find that heat energy is produced at the 
rate of about 400 watts. It is shown in §8 that this corresponds to 
the point of maximum evaporation of a 5th magnitude meteor. If the 
heat energy is distributed over a radius r (which may be of the order of the 
dimensions of the column of ionization), then the increase, 57’, in the 
temperature of the column will be 


ST’ (ura) /(37Br2n,k). eA ie on Pee (317) 


The possible effects of the increased temperature on the radio echoes will 
be discussed in § 5. 

Tn addition to the reduction in the volume density of ionization due to 
diffusion, electrons may be lost through recombination with positive ions 
and by attachment to neutral molecules, followed by mutual neutralization 
of the positive and negative ions. We will deal with these separately. 
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(b) Recombination 
Instead of (3.3) we now have 
én, Da an, 7 
~~=-—=|('=)- jets 4 oar 
Ga SW 5 (7 7) ht oa) 
where «, is the electronic recombination coefficient (and should not be 
confused with the electron line density, «). 
From (3.8) the rate of change in the electron line density, «, is 
da. a 
ine — 2H, |, rn,?dr. 
To a first approximation we may assume that the variation of n, with r is of 
the same form as (3.6) and thus obtain 
a t 


=l ge te ——— In = at) oe eee 


~~ pa le ee 
n= t[47D-+- daa, In (t/to)] 


where (4D¢,)1/2 and %» are the initial radius and line density respectively. 


exp[—r2/(4Dt)], . . . (3.10) 


(c) Attachment 
We now have 


eo = # e ( o —f,n, %»,+detachment terms, . . (3.11) 
where f, is the coefficient of electron attachment. 

We can obtain an upper limit to the effect by assuming that each 
electron attached to a neutral molecule is permanently removed ; hence, 
neglecting the detachment terms and integrating (3.11) with respect to 
radius, we have 

a=, exp (—f, n,, t), +o. fs <u, eee tL ogee) 
and since n, retains the same form as (3.6) 

% 
47 Dt 
where 7,,, is here the density of the negative ion forming molecules. The 
magnitude of the effects of the above: processes on the character of the 
radio echoes is discussed later. 


exp [—r?/(4Dt)—B,n,, #], il ehe dy fed Locka 
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§ 4. THEORY OF THE SCATTERING OF Rapio WAVES FROM METEOR TRAILS 


Pierce (1938) suggested that meteors form a cylindrical column of 
ionization in which the electron volume density is greater than the critical 
density out to a radius which is large compared with the wavelength, in 
which case the reflection is similar to that from the ionosphere. This 
assumption, however, leads to scattered energies considerably greater than 
that found for the majority of echoes obtained with equipment operating 
in the wavelength range 4-10 m. 


Radio Echo Studies of Meteor Ionization 501 


Blackett and Lovell (1941) derived the reflection coefficient for a 
cylindrical column of diameter small compared with the wavelength, 
assuming the electrons to scatter coherently and independently. This has 
been applied to the case of meteor trails by Lovell and Clegg (1948) who 
predict a variation of echo amplitude from the extending column with time 
as given in fig. 2. The echo exhibits a relatively slow rise, reaching 0-5 


Fig. 2 
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Echo amplitude vs. time for an extending narrow column of ionization. The 
meteor reaches the intersection of the normal from the observing station 
with the trail at the instant t=0. 


times its final amplitude the instant the meteor reaches the intersection of 
its path with the perpendicular from the observing site (we will term this 
the ‘reflection point’), and subsequently fluctuates about the final 
amplitude with a period dependent upon the meteor velocity and range and 
the wavelength. Measurements of the period of the rapid fluctuations 
enable the velocity of the meteor to be determined (Davies and Ellyett 
1949). Taking the backward scattering cross section of a free electron to 
be 47(e2/mc?)? where c is the velocity of light, the power delivered to the 
receiver terminals by a meteor trail of line density, «, at range R, is* 


PG23 e2 \2 
e= 3272h? (<3) Oo, yet atn SoG mee Aa (4.1) 


hence the echo amplitude is proportional to « P—=transmitter power, 
G=aerial gain in the direction of the reflection point (referred to an 
isotropic radiator), A=wavelength. 

The part of the complete meteor trail which is effective in reflection is 
restricted to a few Fresnel zones in the vicinity of the reflection point, 
and in evaluating the echo decay we may assume to a first approximation 
that the trail is formed instantaneously and subsequently diffuses radially, 
maintaining cylindrical symmetry. Assuming the incident wave to 


* The formula quoted by Lovell and Clegg (1948) is in error by a factor /2 
due to these authors taking GA2/87 instead of GA?/47 for the aerial absorption 


cross section. 
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penetrate throughout the column, and using the electron density distri- 
bution (3.6) we find that the echo amplitude decays from its initial value 
(given by the Lovell-Clegg formula) according to the law 


exp [—(16s2Dt){A2]/) 9 as ee ee 


Thus we expect an exponential decay (with the fluctuations of fig. 2 
superimposed) with duration 7, defined as the time for decay to exp (—1) 
times the maximum amplitude, given by 

t= A?/(1672D). 2 Od ot en een hee 


The radius of the trail is represented by (4Dt)/?, hence (4.3) shows that 
the echo amplitude is reduced by exp (—1) when the radius of the trail 
reaches A/27.. 

If collisions of electrons are taken into account, the received echo power 
is less than (4.1) by a factor [1+-4(v,/w)?]-1 where w=2z7 times the radio 
frequency. Substitution of numerical values shows that collisions have 
negligible effect at the frequencies commonly employed (30-80 Mc/s). 

In addition to the above theories, various authors (Herlofson 1951, 
Feinstein 1951, Kaiser and Closs 1952) have predicted that when the 
incident electric vector is normal to the meteor trail enhanced scattering 
due to a resonance phenomenon may occur. 

A comprehensive theory of scattering of radio waves by meteor trails 
has been given by Kaiser and Closs (1952) (see also Eshleman 1952) and 
is outlined in the following. 

We may treat the trail as a column of dielectric constant 
47n,e7 


aaah eviet (42) 


k= 1 
where k=27/A, and if « is at all times close to unity, the Lovell—Clege 
formula is valid. In most cases, however, the axial dielectric constant will 
be initially negative, and only after diffusion has gone on for some time 
will it approach unity. For a narrow trail with a large negative axial 
dielectric constant, the condition for penetration of the incident wave is 
that —x(kr)?<1 throughout the trail. Using (3.6) and (4.4) the limiting 
case is that for which 


9 


ta (5) exp (—1)=1 ran Sper 


or, substituting numerical values, «2-4 1022 em-1. 
Thus the echoes fall into two categories depending upon whether the 
line density is greater or less than this value. 


(a) Decay Type Echoes, «<2-4 1012 em-1 


The incident wave penetrates the column and when the incident electric 
vector is parallel to the trail (parallel polarization) the electrons scatter 
coherently and independently, and the Lovell—Clegg formula is valid. 
For the electron density distribution (3.6) the decay of the echo 
is exponential with the decay constant, 7, given by (4.3). 
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When the electric vector is transverse to the trail, a restoring force acts 
on the electrons, and resonant scattering may occur. This is most easily 
seen by first considering a homogeneous cylinder of ionization. If we 
polarize the column by displacing the ions and electrons by a small 
relative amount, x, then each electron will experience a restoring force 
towards its initial position equal to 27n,ex, and hence the resonance 
occurs at a frequency w/27 where w?=27n,e2/m, i.e., when the dielectric 
constant x=—1. The main factor which limits the magnitude of this 
resonance is radiation damping of the cylinder oscillations which becomes 
severe, and the resonance disappears, when the radius becomes of the order 
of A/27. For a homogeneous cylinder, with «=—1, r=A/2z, we find 
a=2x10". Thus «<2-4 x 10! is also the condition for the existence of 
the resonance. 


Fig. 3 
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Resonant scattering from a cylindrical column with ne=n, exp [—(r/79)°]. 
(a) (91/9u)res> (6) Kyes- 


The case of a diffuse cylinder of ionization, as defined by (3.6) is more 
complicated. In this case the resonance with transverse polarization is 
limited mainly by collision damping in the vicinity of the region where « 
goes to zero. This occurs because the continuity condition for the normal 
component of the electric displacement at the boundary «=0 requires the 
electric field to become large, and even if the collision frequency is small 
compared with w, the damping may be considerable. In the approxi- 
mation of a continuous dielectric medium, it remains finite even when 
y, > 0 since the electric field goes to infinity. In this case the magnitude: 
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of the resonance depends on the gradient of « near x0, except when « 
approaches 1012 cm~! when radiation damping becomes effective and the 


resonance disappears. 
Kaiser and Closs (1952) have derived the echo characteristics with 


transverse scattering for an electron density distribution 
n, « exp [—(?/"o)*] 

and have considered the case s=2, 7)2=4Dt (i.e. the density distribution 
of (3.6)) in detail. Figure 3 gives (g,/Gi)res ANA Kyes VS. 8 Subject to 
a< 2-4 101, where g, and g, are the reflection coefficients (defined below) 
for transverse and parallel scattering, (7,/9i);.; is the value at resonance, 
Kes is the axial dielectric constant for the resonant condition. It is clear 
that the sharper the boundaries of the column, the more pronounced is the 
_ resonance. 


Fig. 4 


i) iO722 s 2x1O7!2 Sxi0es 


Polarization ratio, g,/g,, as a function of 4k?Dt/« for various values of a. 


For the important case, s=2-0, the resonance occurs when the axial 
dielectric constant is —1-4 and its effect is to increase the maximum 
amplitude of an echo obtained with transverse scattering by a factor 2-0 
above that given by the Lovell-Clegg formula. Figure 4 gives g,/g, as a 
function of 4k?Dt/« for s=2, with various values of «. The effect of the 
radiation damping as « approaches 10” is clearly shown. Subsequent 
to the resonance the axial electron density will continue to fall and when it 
is somewhat less than the critical value, the parallel and transverse 
reflection coefficients will become identical and the echo will decay 
exponentially according to (4.2). 

Owing to the finite meteor velocity, the reflection coefficient will not, 
in general, be constant along the length of the meteor trail and if its 
variation is significant over the length occupied by the central Fresnel 
zone it will affect the echo amplitude. The reflection coefficient at a given 
point on the trail at a time ¢ after the passage of the meteor is as given by 
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the theory of Kaiser and Closs, hence the echo amplitude may be deter- 
mined by numerical integration (for transverse polarization it is necessary 
to take account of the phase shift in g, with t). This has been carried 
out for the case of parallel polarization and the results are given in figs. 5 
and 6 which show the echo envelopes for various values of the quantity 
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Echo envelopes (with parallel polarization and «<2-4 10}? cm-*) for 
various values of 4=167?DR1/?/(vA3/?), 


Fig. 6 
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Maximum echo amplitude (from fig. 5) vs. 4. 


A=1672DR1?/(vi3/2) and the variation of the maximum echo amplitude 
with 4. As an example, with R=150km, A=4m, v=40 km sec}, 
D=2~x 104, we find 4=0-4. It can be shown that the periods of the 
fluctuations are not significantly affected, and hence meteor velocity 
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Reflection coefficient as a function of 4k?D#, evaluated on the assumption that 
the meteor instantaneously forms a narrow ionized column which 
expands radially at a rate determined by the diffusion coefficient, D. 
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theoretical amplitude—time envelopes of the echoes. 
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measurements remain reliable, provided 4<2. The case of transverse 
polarization is more complex. Kaiser and Closs (1952) suggest that for 
the above parameters, the effect will be to reduce the observed resonance 
for electron line densities less than a few times 1011 per cm. 


(b) Persistent Echoes 


The incident wave does not penetrate throughout the column and the 
reflection coefficient is obtained by assuming total reflection from a 
cylinder of radius r, given by the larger of the two values satisfying 
K(r,)(kr,)2=—1. For the case (kr,)? >1, 7, is just the radius at which «=0. 
Reflection from the diffusing column of ionization is thus similar to that 
from an expanding metallic cylinder of radius 7,, and the echo persists until 
the axial electron density falls below the critical value. 

Figures 7 and 8 give the echo characteristics derived by Kaiser and Closs 
on the assumption that the meteor instantaneously forms a narrow 
column of electrons which expands radially due to diffusion, The curves 
of fig. 8 give the expected amplitude—time envelopes of the echoes. The 
reflection coefficient, g, is defined so that the received echo power is 


P7Gr 
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The duration and maximum reflection coefficient of the persistent echoes 
are given by 


Las (2 4.7 
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The theory of Kaiser and Closs thus reconciles the theories of Pierce 
and of Lovell and Clegg. Two echo types are predicted, decay and 
persistent, depending upon whether « is less or greater respectively than 
2-4 1012. In the former case the echo amplitude is proportional to «, 
whereas the duration is independent of «, being determined by the 
wavelength and diffusion coefficient. The persistent echo amplitude 
varies only as the fourth root of «, but the duration is proportional to «. 
In both cases the duration is proportional to the square of the wavelength 
and inversely to the diffusion coefficient. 


(c) The Head Echo 


In addition to the specularly reflected or ‘ body’ echo, the properties” 
of which are predicted above, a very short duration echo is sometimes 
seen, apparently originating at the head of a dense trail and moving with 
the velocity of the meteor (Hey ef al. 1947, Millman and McKinley 1948, 
McKinley 1951). 

2M2 
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Millman and McKinley (1949) have attributed the head echo to ioniza- 
tion instantaneously produced by ultraviolet radiation from the meteor. 
The difficulties inherent in this explanation are firstly that the ionization 
would need to disappear almost instantaneously to explain the short echo 
duration (usually too short to be measurable) and secondly an appreciable 
proportion of the meteor energy would have to be dissipated in this way in 
order to produce sufficient ionization to explain the observed head echo 
amplitudes. 

A simple theory has been propounded by Browne and Kaiser (1953) 
which is based on the well established properties of the meteor trail and 
which gives good agreement with observation. The head echo received as 
the meteor approaches the reflection point is easily explained. Its 
amplitude is given by the rising portion of fig. 2 and its range, of course, 
must be the instantaneous range of the meteor. The head echo received 
after the meteor has passed the reflection point (the receding head echo) 
is explained as follows. The incident pulse of radiation is in the form of a 
spherical shell which first illuminates the trail in the vicinity of the 
reflection point and the scattering from this region gives rise to the body 
echo. Subsequently the spherical shell intersects the trail in two short 
strips and the nett scattered energy towards the observer may be shown to 
be zero. Finally, however, the shell intersects the discontinuity at the 
head of the meteor trail and there is a nett scattering of energy giving rise 
to the head echo. The amplitude, A,, of the head echo from a meteor at 
range Ff, is related to the amplitude A of the body echo which would be 
received from a trail of the same line density at the same range by 


A 1 /A\12 
3-57) ’ . . . a tate ty © (4.9) 


where d= f,,—, R being the range of the specular reflection point on the 
meteor trail. Thus the intensity ratio of the head echo to the body echo 
should be proportional to the wavelength and inversely proportional to the 
range difference between them. The head echo amplitude will also vary 
with the electron line density «, at the head of the trail either as «, or 
o%,1/4 depending upon whether «, is less or greater than a few times 
10%cm~!. The latter case will usually apply since the head echoes can 
only be detected from meteors producing very high line densities. 

The majority of the echoes observed on radio echo equipments operating 
in the wavelength range 4-10 m are of the short duration decay type and 
are of the predicted form. The persistent echoes are also observed, with 
durations up to several hundred seconds on 8m, and often the echo 
structure is complex, with discrete echoes at several ranges in addition to 
the one originating from the specular reflection point. Browne and 
Kaiser (1953) have pointed out that these may be produced in a similar 
way to the head echoes, deriving from discontinuities along the meteor 
trail produced by bursts of ionization associated with the visual flares 
observed with very bright meteors. A more important mechanism has 
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been established by Ellyett (1950) and Greenhow (1950) who attribute the 
effect to distortion of the meteor trails by non-uniform winds in the upper 
atmosphere in which case there may be several regions of a trail which are 
effective in reflection. If the range difference between two such regions 
is less than the range discrimination of the apparatus, the interference 
between waves reflected from various regions in relative line of sight motion 


gives rise to the deep fading which is commonly exhibited by enduring 
echoes. This will be dealt with in § 6. 


§5. EXPERIMENTS ON THE SCATTERING OF RADIO WAVES FROM METEOR 
TRAILS 


Closs, Clegg and Kaiser (1953) using a 72 Mc/s equipment with twin 
aerials designed to illuminate a meteor trail alternately with transverse 
and parallel polarization, have demonstrated that the majority of the 
decay type echoes, exhibit enhanced scattering with transverse polariza- 
tion. The mean resonant polarization ratio, 7=(9,/9u)res, for these 
echoes was equal to 1-7. Previous measurements by McKinley did not 
support the resonance theory, however, the smallest polarization ratio 
which McKinley could have detected was of the order of 4:0. 
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Resonant polarization ratio, 7, as a function of electron line density. 


Theoretical. 
© Experimental. 


The electron line density of the individual decay type echoes may be 
estimated independently of the echo amplitude from the decay constant 
and the instant t,,, at which g,/g, reaches its maximum (resonant) value. 
A similar method can be applied to echoes exhibiting a rapid rise to a 
saturation value followed later by an exponential decay. It is assumed 
that these correspond to values of « above 10! cm™? (see fig. 8) and the 
line densities may be estimated from the ratio of the total echo duration to 
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the decay constant. Closs et al. determined in this way the variation of 
with electron line density leading to the result depicted in fig. 9, which 
contains, for comparison, the theoretical curve derived from $4. The 
agreement between experiment and theory is remarkably good ; for the 
smallest values of «(~2x10™) the experimental value of 7 tends to 
about 2-1 which is within experimental error of the theoretical value of 
2-0 for a trail with a Gaussian electron distribution. 

This result has a twofold significance. (a) Kaiser and Closs assumed the - 
Sellmeyer formula for the dielectric constant in a plasma in arriving at the 
theoretical value 7=2-0. If the Lorentz formula containing the polariza- 
tion term had been used, « would no longer be a Gaussian function of r 


Fig. 10 
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Mean echo amplitude (corrected to R=100 km) as a function 
of electron line density. 
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Theoretical. 


(for a Gaussian electron distribution) and the gradient of « in the vicinity 
of «=0 would be increased, leading to 74. It would thus appear to be 
ruled out by the experimental result. (b) If the heating of the trail were 
significant, the temperature would be a decreasing function of radius and 
time, leading to a trail with sharper boundaries, hence again requiring 
7>2 in contradiction with experiment. 

We may compare this apparent absence of a temperature effect with 
prediction. Using (3.7) together with the theory of §4 we find that the 
increase in the temperature of the trail at the instant of resonance would 


be of the order of 
dT 6:5 x 10-*8 pv?/(Bn,A?). 
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Putting 4=50a.m.u. and taking 0-1 as a lower limit to B (see $11), 
v=40 km/sec, we obtain 87'=3 x 10%3 n,-! which is small compared with 

the temperature of the surrounding atmosphere (~200°) provided n, >101 
which is well satisfied throughout the meteor ionization region. Thus the 
temperature effect should be negligible for all decay type echoes. In fact 
this applies with even greater force to the enduring echoes, since although 
these are produced from trails with high densities, hence producing heat 
energy at a greater rate, we find that the trail temperature can be 
significantly greater than that of the surrounding atmosphere for only a 
small proportion of the total echo duration. 

Closs et al. have determined the echo amplitude vs. electron line density 
for transverse and parallel polarization. The result is given in fig. 10, 
together with the theoretical values. The agreement is very satisfactory, 
and the saturation in echo amplitude due to the change in the echo 
characteristics at «~101 is particularly well marked. The smallest line 
density group in fig. 10 is .=1-6510'1cm-!. The mean range of these 
echoes was 130 km and the echo amplitude with transverse polarization 
was about 5xnoise. This may be compared with the value «=1:5x 104, 
for 5 x noise echo, calculated from the equipment parameters. 

From the decay of the Arietid echoes, the most probable value for the 
diffusion coefficient was found to be 2:3 104+. The most recent measure- 
ments of meteor heights (see § 10) indicate that this corresponds to a height 
of 92:5 km where the scale height and particle density are 6-3 km and 
7x 1013 cm- respectively. Hence, from (3.5) we obtain A,=0-45 cm at 
92-5 km, and the mean collision diameter, o,=8:5x10-§cm. Thus, in 


general D=2x 10% n,-1 HY? 


5.1 
A,3-7X 108 n,-1 om 


and in the meteor ionization region, 
EST a, OG) 


The variation of echo duration with wavelength has been studied by 
Greenhow (1952 a) using equipments on 4 and 8 m. - The results show that 
the duration is proportional to X” where n=2 for most echoes, but 
decreases with increasing duration to about unity for 7~100 sec on 8 m. 
Greenhow suggests that the decrease is due to the influence of atmospheric 
turbulence, however, it is possible that loss of electrons by attachment to 
neutral molecules may have contributed (see §§ 7 and 15). 


§ 6. Fapine or Mzetnor Ecuors anp Upper ATMOSPHERIC WINDS 


Echoes of duration in excess of about a second rarely exhibit the ideal 
form given by the theory of § 4 but in most cases are subject to deep fading, 
sometimes simply periodic, more often complex. These amplitude 
fluctuations have been studied, with pulsed equipments, by Ellyett (1950) 
and Greenhow (1950, 1952 a, 1952b). The period of the fluctuation is 
proportional to the radio wavelength and for A=8 m, varies between 0-04 
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and 0:4 sec. The periodic fluctuations may be explained as due to reflec- 
tion from two distinct points on the trail, which are in relative motion and 
for the more complex echo patterns the echo amplitude follows a Rayleigh 
distribution consistent with the existence of a multiplicity of moving 
reflecting centres. Greenhow (1952 a) has suggested that these centres are 
formed through distortion of the meteor trail by non-uniform winds in the 
upper atmosphere. Strong evidence for this mechanism is that the 
complexity of the echo pattern tends to increase during the echo lifetime, 
consistent with increasing disorder along the trail. The distortion may be 
sufficient to enable an echo to be received from a trail which initially did 
not produce ionization in the vicinity of the intersection between the 
normal from the observing station and the path of the meteor. In this 
case there is a delay between the passage of the meteor and the appearance 
of the echo. If a sufficiently short transmitter pulse is used the complex 
echo structure manifests itself in a broadening of the echo pulses which are 
sometimes even split into several components. 

From the fluctuations and range drifts of long enduring echoes, 
Greenhow (1952 b) has obtained information concerning the properties 
of the upper atmospheric winds. The observations were made at 
Jodrell Bank, Cheshire, England (latitude 53°) during the Geminid 
meteor shower, December 12th—14th, 1950. Between 2300 and 2400 hr 
G.M.T. winds of velocity 30-60 m sec~! appeared to be blowing from a 
northerly direction in the height range 80-90 km while after 0320 hr 
the wind direction had reversed. The velocity gradient ranged between 
1 and 10 m sec~! km~1, thus the velocity at heights separated by 5-10 km 
may differ by as much as 50 msec}. | 

Manning, Villard and Peterson (1950, 1953) at Stanford, U.S.A., have 
used a continuous wave equipment and observed echo fluctuations due to 
interference between the ground wave and the wave reflected from the 
moving trail, as well as those due to the multiplicity of reflecting centres. 
They obtained (1950) a mean wind velocity of 30-40 msec! with a 
predominant SSW-NNW direction and since the measurements were 
carried out after 0300 hr (local time) they are in general agreement with 
those of Greenhow. 

Elford and Robertson in South Australia (private communication) 
have studied the effects of winds on meteor echoes with a combined pulse— 
C.W. apparatus which enabled them to measure simultaneously the slant 
range and direction of arrival of the reflected wave. They have obtained 
extensive results for the months October-December 1952 over the height 
range 80-105km. The mean wind velocity is found to be of the order 
of 60m sec with a mean positive wind gradient of 3-6 m sec-! km-1. 
The results indicate a prevailing component towards E or NE and a 
periodic component with a 12 hr and a 24 hr variation, both representing 
an anticlockwise rotation of the velocity vector. The diurnal and semi- 


diurnal components were directed to the north at 1-8-+-2 and 0+1 hr 
(local time) respectively. 
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§ 7. RECOMBINATION AND ATTACHMENT IN METEOR TRAILS 


The effect of recombination in meteor trails has been briefly discussed. 
by Kaiser and Greenhow (1953) who conclude that the recombination 
coefficient is not greater than 5x 10-1 cm’ sec-! i.e. considerably less 
than the value of 10-* observed in the E region of the ionosphere. 

An upper limit to the recombination coefficient may be obtained from 
a consideration of the observed echo durations. We have seen (§ 4) that 
provided the initial electron line density is large compared with 1012 cm-}, . 
the echo will persist until the axial volume electron density falls to the 
critical value. Now from (3.10) the axial electron density at time ¢ 
after the formation of the trail cannot be larger than 


2[% 6 In (t/t) 1-4, 


where (4Dt,)1/? is the initial radius of the trail and will be of the order of 
the molecular mean free path, A,. Equating this to the critical density, 
we find that the maximum duration, r,,, which can be observed is given by 


ue nes 
rein (ADra[Xa—=— (Ha). s (TH) 
Since A, will be a few times greater than X,, we find, from (5.1), D ~ 104A, 
and In (4D7,,/A,7)=In (4x 1047,,) for reasonable values of A,. Observa- 
tions on 4:-2m (Greenhow 1952) show that 7, must be greater than 
60 sec, hence, from (7.1) we find «,<3-5 x 10-4 cm3 sec71. 

McKinley and Millman (1953) have reported echo durations of up to 
30 min on A=9-1m which will reduce the upper limit set for the re- 
combination coefficient to «,<4:5 x 10-12 em? sec=}. 

Theoretical calculations of the radiative recombination coefficient for 
various ions lead to values for «, of the order of 10-¥ and Massey (1952) 
suggests that for thermal electrons it is never likely to exceed 10-1. 

The discrepancy between the value of «, obtained in the E region of 
the ionosphere and the radiative recombination value has been explained 
on the assumption that the E region ionization is mainly molecular 
(O,+) and that dissociative recombination with a coefficient of the order 
of 10-8 is predominant (Bates 1950). 

The fact that the ionization potentials of the meteor atoms are con- 
siderably lower than those of the atmospheric constituents, and the 
absence of spectral evidence for atmospheric ions suggests that the 
ionization is mainly atomic. Hence, provided charge exchange between 
the ions and neutral molecules is negligible, dissociative recombination 
will be absent and we would expect the recombination coefficient to 
approach the radiative value, which is in agreement with the above 
estimates. 

We can estimate an upper limit to the effect of electron attachment to 
neutral molecules (which we will assume to be molecular oxygen) from 
(3.13) which shows that attachment must be negligible if Bm_!< 1 where 
Nm is the molecular oxygen particle density. From (5.1), with 
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H=6-5x 105 em, the condition that decay type echoes should be un- 
affected by attachment becomes 


Ny< LO24(B Ae aes 


On a wavelength of 10m, with B,=10-1 cm sec? (Mitra 1952) this 
gives n,<10 cm-* corresponding to a height of about 75 km which is 
below the normally observed height range of decay type echoes. For 
A\=100 m, the above condition becomes n,<1014 em? corresponding to 
a height of about 90 km, in which case the decay type echoes may be 
seriously affected. 

The condition that attachment be negligible for enduring echoes is 
likewise B,n,,7< 1 where 7 is the echo duration. The smallest meteor 
which will be affected is thus that for which 8,n,,7 becomes of the order 
of unity at the point on the trail at which, for 7 given by (4.7), «/D is 
a maximum. The relation between the maximum observable echo 
duration 7,, and &,,x for a given meteor is derived in §11 (eqn. (11.5)), 
hence we find that the smallest meteor which will be significantly affected 
will be that which in the absence of attachment would be capable of 
giving a maximum echo duration 


7221:6%10-™f Ae (B,* cos +) 17/o 5 en ee 


and this will be of the order of the shortest echo duration for which 
attachment is likely to be significant. It must be realized that (7.2) is 
only an order of magnitude relation and in particular gives no information 
as to the magnitude of the effect of attachment on the echo duration. 

As an example we will take A=4x10? cm, v=4xX108 cm sec}, 
cos?/> y~1 giving 

(a) 8 .=10-*;, cl sec (6) B,=10-1*, 790 sec. 

Thus we might expect only the longest duration echoes to be influenced. 

The main observable effects of attachment will be to reduce below 2-0 
the exponent in the power law relating echo duration and wavelength 
(this may have been a contributory factor for the echoes observed by 
Greenhow (see §5)) and to decrease the number of echoes with durations 
greater than that given by (7.2) (see § 15). 

It is clear from (7.2) that the larger the meteor velocity, the smaller 
the effect of attachment, also that meteors travelling in near horizontal 


paths are less affected. This latter is because they ionize at relatively 
greater heights. 


§ 8. CORRELATION witH VisuAL METEORS 


Since both the excitation, which gives rise to the luminosity of a 
meteor, and the ionization are produced by collisions between evaporated 
meteor atoms and air molecules, we expect them to be proportional. 
The visual magnitude scale is defined so that a change of one magnitude 
represents a factor of 2-5 in luminosity, and is referred to the point of 
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maximum luminosity, which will coincide with that of maximum 
ionization. Hence 


M=—2-5 logy %maxteonstant, : . . . . (8.1) 


where M is the visual (zenithal) magnitude. 

Simultaneous visual—-radio echo observations have been made by 
various groups of workers (Prentice, Lovell and Banwell 1947, Millman, 
McKinley and Burland 1948 (see also Millman 1950)). The most 
extensive results (700 coincidences) are those of Millman et al. during the 
Perseid meteor shower of 1947, who obtained a linear relation between the 
logarithm of the radio echo duration and the apparent meteor magnitudes. 
The results are reproduced in fig. 11. The dotted line gives log 7 as a 
function of zenithal magnitude, M, and is obtained by assuming a constant 
correction of one magnitude. The wavelength used in these observations 
was 9-1m. The durations in fig. 11 represent a mean for each magnitude 
group, however it is not stated whether this is the arithmetic mean or 
the mean logarithm of the duration. 


Fig. 11 
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Duration of 9:1 m echoes from Perseid meteor trails as a 
function of visual magnitude. 


Using (8.1) we find that the mean duration, 7, in fig. 11 is related 


to max by - 
GEC Gp eee oe ees teeta oie) 
and it is of interest to compare this with theoretical prediction. 
It is shown in §11 (eqn. (11.5)) that the maximum persistent echo 
duration, 7,,, which it is possible to observe from a meteor of given magni- 


tude is related to &nax by 
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Greenhow and Hawkins (1952) assumed the observed mean duration 
to be of the same order as 7,, and that the diffusion coefficient would have 
the same value (equal to that found for decay type echoes) for all echoes, 
ie. that 7,, Coax. They explained the discrepancy with the observed 
relation (8.2) by suggesting, on the basis of early calculations by Opik 
(1937), that the luminous efficiency increases with the brightness of the 
meteor, while the ionizing efficiency remains constant, thereby deducing 
that a 6th (zenithal) magnitude meteor produces a maximum of 10% 
electrons per cm, and, by assuming Opik’s values for the luminous 
efficiency that the probability of ionization, 8, for a 60 km/sec Perseid 
meteor is of the order of 0-2. 

When we take account of the variation of the diffusion coefficient the 
extent to which the relation between luminosity and ionization would 
have to deviate from linearity would seem impossibly large and so it is 
worthwhile to look for other effects which may reduce the mean observed 
duration relative to 7,, for the bright meteors. Amongst these are 
attachment of electrons to neutral molecules and the disruption of the 
trails by non-uniform winds. From a calculation similar to that made 
in §7, we find that attachment would commence to be effective at a dura- 
tion lying between 35 and 200 sec for an attachment coefficient between 
10-15 and 10-16. The fact, however, that the observed relation between 
log + and M remains linear up to z=100 sec could mean that attachment 
is not important over this range. The variation of the exponent in the 
duration—wavelength law observed by Greenhow (see § 5), suggests that 
one or other of these processes is occurring, but an estimate of the 
magnitude of the effect from these results indicates that it is unlikely 
to be sufficient to account completely for the discrepancy between (8.2) 
and (8.3). 

Another explanation for the discrepancy is that the relation between 
the observed mean duration and 7,, will not be linear. The form of this 
relation will depend on the selection of echoes involved in obtaining the 
correlations and since Millman et al. have presumably concerned themselves 
only with the persistent echoes, the selection may in effect be that only 
echoes with durations greater than some value 7,,;, were accepted. The 
value of 7,,;, Will lie between the decay type echo duration (~ 0-1 sec 
for the equipment used) and the smallest duration given in fig. 11 
(=I sec). Assuming such a selection, we would expect the increase of 
T with %,, to be less rapid than that of 7,,, since for a given meteor the 
observable range of echo durations will extend from 7,,,, to Tm: Lhe 
magnitude of the effect will depend somewhat on how the durations were 
averaged, If 7 is taken as the mean logarithmic value, and if it is suffi- 
ciently in excess of 7,,;,, we find approximately, 7 ocz,,!/2, which alone 
is sufficient to reconcile (8.2) and (8.3). 

In view of the above, we will take the luminosity to be directly propor- 
tional to «, and estimate the value of the constant in (8.1) from the values 
of 7 given in fig. 11 for the faintest meteors observed. In this way we 
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find that «,,.=10! corresponds to a 5th—6th (zenithal) magnitude Perseid 
meteor. A similar result is obtained (Greenhow and Hawkins 1952) by 
extrapolating the straight line of fig. 9 to the observed decay type echo 
duration. A comparison of visual and radio echo rates (§16) favours 
the former value, which we will adopt, hence for Perseid meteors (8.1) 
becomes 


ee ele ogi HON, G Nica sae. eos (8.4) 


It is interesting to note that the faintest visible meteor produces a 
maximum line density corresponding to the transition between the decay 
and persistent echo types. 


The theoretical relation between M and «,,,, may be shown to be 
M=25-35—2-5 logy, (KuvFanax/B), - . » - (8.5) 


where K isthe fraction of the kinetic energy converted into visible radiation. 
Hence, for Perseid meteors (v=60 km/sec), taking 150 a.m.u., 


eo 20x LOS, 
If we had taken «,,,,=10-! for M=6 we would get K/B=3 x 10-3. 


§ 9. THEORY OF THE MeTEOoR Hetcut DisTRIBUTION 


We will present here the results of a theory of the meteor height 
distribution (Kaiser, unpublished) based on the meteor ionization 
equations of §2 and the reflection properties of meteor trails. The 
discussion is restricted, in that the radio echo equipment sensitivity is 
assumed to be such that the majority of the echoes are of the decay type. 
The mass-frequency distribution of the incident meteors is assumed to be 


Vein ONt =. CU SOC) sh aoe hte > (9-1) 


where b is a constant which defines the absolute meteor flux. Two cases 
are considered : (a) shower meteors constituting a homogeneous velocity 
group with a well defined radiant point at zenith angle y; (6) a selected 
velocity group of sporadic meteors with reflection points in a narrow range 
of elevations about elevation $,. 

Meteors from a given radiant point travel in parallel paths, hence, due 
to the specular reflection properties of the meteor trails, the reflection 
points will all lie in a plane normal to the paths of the meteors and 
intersecting the observing station. If7/2—0 is the angle in the echo plane, 
measured from the horizontal, the normalized height distribution observed 
in the element 6—(6+d6) for meteors of the selected velocity, v, is 


Ber, 3(s—1) = 
Da (ty) = Es . ° ° ° ° ° (9.2) 


where «,=(h—h,)/H, h, is the height. at which the pressure, p,, and the 
minimum detectable line density, «,, at the mean observed range are 
related by (2.5) (with p,, « replacing Pmax) max) The absolute echo rate 
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Fig. 12 


1-O eS 2:0 rads 320 


(a) f(u)=(1+e7u?)-. (c) f(u)=cos? cu. 
(b) f(u)=exp (—c?u?). (d) S8(6@)=cos [(z/2) sin 6]. 


Mean height, 4 ; shower meteors. 


Fig. 13 


(a) S(9)=(1+c?S*)-. (b) S(#)=exp (—c29*), (c) S(d)=cos? cd. 
© A/2 dipole, 4/4 above ground, ¢,,=40°. 
Mean height, 4 ; sporadic meteors. 


ae 
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in the element 6—0+-d@ is approximately proportional to 
[cos®/? 6 S?(A)]s-1 cos 6, 


where S(6) is the aerial polar diagram (amplitude) in the echo plane. 

The complete height distribution for shower meteors is obtained by 
integrating the elementary height distribution (9.2) over the echo plane, 
taking aecount of the variation of h, (though «,) and of the absolute rate. 
That for sporadic meteors is obtained by further averaging the height 
distribution over all echo planes intersecting the aerial beam, taking 
account of the sporadic radiant distribution. The height distribution is, 
however, found to be quite insensitive to the exact form of the geocentric 
radiant distribution which will therefore be assumed uniform. 

The properties of the height distribution are as follows. For shower 
meteors, the mean height, h, is close to the height h,, at which the pressure, 


Dm, 18 given by 9 9 1/3 
¥ Pn= g@ (j Hoy cos? x) Be U0, Ay et tio EES) 


y2 


Fig. 14 


=(1+c?u?)—. (c) f(u)=cos? cu. 
8 ee (— chu?) : (d) S(@)=cos [(7/2) sin 6]. 


R.m.s. height deviation, 5h ; shower meteors. 


where «,, is the minimum detectable line sensitivity at the height hin the 
direction of the aerial beam maximum in the echo plane. A similar relation 
holds for a selected velocity group of sporadic meteors with reflection 
point elevations occupying a narrow range about elevation ¢,, except 
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that y is now replaced by ¢,,, and «,, is taken as the minimum detectable 
line density at range h/sin ¢,,, in the direction of the aerial beam maximum 
for elevation ¢,,. 

The small difference, (h—h,,), depends slightly on the form of the aerial 
polar diagram and on the exponent, s, in the mass distribution. It is 
given as a function of s in figs. 12 and 13 for shower and sporadic 
meteors respectively, with various assumed forms for the aerial polar 
diagram. For shower meteors, the polar diagram is assumed to be 
symmetrical about 6=0 and is defined in terms of 


f(u)=cos* 6 S(@), u=tan 6, -. 55s 3 (ee 


while for sporadics, S(#) represents the polar diagram in azimuth, # at the 
constant elevation ¢,,. 


Fig. 15 


30 


(a) S(B)=(1+c29)-1. (6) S(9) =exp (—c292).  (c) S(B)=cos? cd. 
© A/2 dipole, A/4 above ground, ¢,,=40°. 
R.m.s. height deviation, 8h ; sporadic meteors. 


The width of the height distribution, defined as the r.m.s. deviation, 5h, 
from the mean height, is a function only of the atmospheric scale height. 
and of the exponent in the meteor mass law. Figures 14 and 15 give dh/H 
as a function of s for showers and sporadics, respectively, and for the 
various assumed aerial polar diagrams. 

It is clear from figs. 12-15 that the properties of the radio echo height. 
distribution are relatively independent of the form of the aerial polar 
diagram, especially in the important range about s=2 (see § 13). 
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Because of the factor v-2 in (9:3), p,, will decrease, and hence / will 
increase, with increasing velocity ; p,, will also increase slightly with 
increasing «,, i.e. with decreasing equipment sensitivity. 


Fig. 18 
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Temperature vs. Altitude from the radio echo measurements, assuming 
a mean atmospheric molecular weight of 29 a.m.u. 


Fig. 19 
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§ 10. PROPERTIES OF THE UprER ATMOSPHERE FROM THE Murnor HEIGHT 
DISTRIBUTION 

Evans (unpublished), using the meteor height finding equipment 

described by Clegg and Davidson (1950), has measured height distributions 
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for selected velocity groups of sporadic meteors between v=20 and 
60 km sec~! and for several meteor showers. 

During the period January to July, 1952, the heights and velocities of 
381 sporadic meteors were measured. The observing station was located 
at Jodrell Bank, Cheshire, latitude 53°. The observations were divided 
into 7 velocity groups for each of which the mean height and r.m.s. 
deviation, 5h, was determined. The effects of the aerial polar diagram were 
estimated to be such that the appropriate curve for (h—h,,)/H and 8h/H 
would lie between (a) and (b) of figs. 13 and 15. Taking s=2-00 for 
sporadic meteors (see §13) the results of figs. 16 and 17, which give 


Fig. 20 
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Atmospheric density vs. Altitude. 


(a) Radio echo measurements (Jodrell Bank, England, latitude 53°). 
(b) Photographic meteor observations (Massachussetts, U.S.A., latitude 42°). 
(c) Rocket flights (New Mexico, U.S.A., latitude 32°). 


lim V8. v and H vs. h respectively, were obtained. The full curve in fig. 17 
is a straight line fitted by the method of least squares, and the dotted 
line gives, for comparison, the results obtained from rocket flights in 
New Mexico (U.S.A.) at latitude 32° (Rocket Panel, 1952). Figure 18 
gives the temperature-height variation based on the assumption of a 
constant mean atmospheric molecular weight of 29g molt. Due to 
the dissociation of oxygen in this region, the actual temperatures and 
temperature gradient will be somewhat greater than those of fig. 18. 

The absolute value of the atmospheric pressure at the height h,, can 
be calculated from (9.3). This requires knowledge of the function @ of 
the meteor constants and the probability of ionization (which may vary 
with v). Fortunately Q need be known only for a single velocity, giving 


2N2 
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the pressure at one height, from which the pressure—height variation may 
be obtained by integrating the equation 

din ee eee 

dh He 

using the known scale heights. Jacchia (1949) in deriving atmospheric 
pressures from visual meteor observations, has determined the value of 
parameters similar to Q, except that they contain the luminous efficiency, 
K, instead of the probability of ionization 8. The ratio K/B has been 
evaluated for Perseid meteors in §8, which, with Jacchia’s results, leads to 
log Q=4-24 for v=60 km/sec. The resulting pressure height curve is 
given in fig. 19. The individual points in fig. 19 are the pressures derived 
from the various velocity groups on the assumption that Q is independent 
of velocity. The close agreement indicates that the probability of 
ionization does not vary significantly with meteor velocity over the range 
20-60 km sec~1. Figure 20 gives the atmospheric density derived from 
these pressures with, for comparison, results obtained from the rocket 
flights and from photographic meteor observations made at Massachussetts 
(U.S.A., latitude 42°, Jacchia 1952). 


§ 11. DETERMINATION OF THE CONSTANTS IN THE METEOR IONIZATION 
EQUATIONS 
We will first estimate in different ways, the value of the probability of 
ionization, 8. Jacchia (1948) has assumed that the luminous efficiency, 


_ K, varies directly with velocity and gives log,, (K/v)=—9-1 although. 


values as small as —10-0 could be obtained from extreme assumptions. 
Now, in §8, we have estimated that K/B lies between 3 and 7-5 x 10-3, 
hence with log,, (K/v)=—9-1, 8 must lie between 0-6 and 1-6 while the 
smaller value of K/v gives 0:1<B<0-2. 

From (2.5) we have 

B=p{l/(AAQ)P. 
For spherical iron meteors, p=7:8, hence A=0:48. The latent heat, 
estimated from a modification of Trouton’s rule due to v. Wartenberg 
(see Roberts 1948), is J=7-5x10!° erg g™, hence, putting A=0-7, 
logy) V=4-24, u=56 a.m.u., we find B=0-2. For stone meteors both A 
_ and / will probably be larger than these estimates, but we may expect f to 
be of the same order. 

McKinley (1951 b) using the radio echo technique, has measured the 
deceleration of three meteors. One of these was a Perseid with the 
following characteristics. Height of reflection point and enduring 
echo=91km. Echo duration 285 sec. Velocity at reflection point=59 
km sec™', Deceleration a=—1-5 km sec~?. From the expression for the 
rate of loss of mass given by Whipple (1943) we obtain the following 
expression for the electron line density in terms of the deceleration, a, and 
the atmospheric density p, : 


a= 5 XK: LOT Bip.) GAs ie ena eee eke ee 


PORE 
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The value of « may be estimated directly from (4.7) using the observed 
echo duration. From the preceding section we find p,=4:7x10~° for 
h=91 km, hence from (5-1) D=1-7x104, giving «#>2-5x 10 cm-l, 
Substituting this in (11.1) with the above values for a and v we find Bale 
McKinley estimated the electron line density from the echo amplitude 
using the Lovell—Clege scattering formula and obtained a much smaller 
value for 6, however, this was obviously an enduring echo, hence his 
estimate of « was several orders of magnitude too small. 

The somewhat surprising result of the above is that although the value 
of 8 has not been determined precisely, it is clearly of the order of unity 
and almost certainly is not less than 0-1. It enters generally as 61/8 which 
for most purposes can clearly be taken as unity. 

The variation of 8 with velocity is obtained by substituting in (9-3) the 
atmospheric pressures found by integrating (10-1). Assuming that B 
varies as v", we find n=0 + 0-6 for 20<v<60 km sec“!, which corresponds 
to particle energies in the range 100-1000 ev. 

If we take 7 ev as the ionization potential of a meteor atom, the fraction 
of the total meteor energy converted into ionization is 


J=2-24 X 10712 B(uv?)-1, 


For v=60 km/sec, 0:1<8<1, we get 10°<J<10-%. From Jacchia’s 
result, quoted above, the luminous efficiency K, will likewise be within 
this range, however, according to Jacchia, K oc v, whereas from the above, 
ede 
Using the above values for f in (2-4), a meteor producing a maximum of 
1012 electrons per cm (5th magnitude) will have a mass between 10-% and 
10 g, ie., a radius of the order of 0:02 cm. Thus the initial radius is 
related to pax by ° 
ee LOS! ae o emerna sy why te bs (L122) 


The radius at the point of maximum ionization is r,,,,=37" (Herlofson 1948). 
Substituting numerical values in (2.5) (taking H—6-5 km) we find that 
the maximum ionization is produced at a pressure Prax given by 


Dee Ol Oi Con COS: Va 9/0" aiuto nomen (dL Sy 


Taking the molecular diameter for air to be 3x 10~* cm, the mean free 
path in the meteor ionization region is \,=5-65/p cm. ‘Thus at the point of 
maximum evaporation A,>?max for &max<4'10~° v? which represents the 
limit to the validity of the meteor ionization theory. For velocities of 60 
and 40 km sec-!, this corresponds t0 &max<10!, M>—5, o&nax<2°5 X 10%, 
M>—3 respectively. The limit may alternatively be expressed as 
Pmax< 1:3 X 108 v-}, ie., for v=40 kmsec™?, Pinax<30 dynecm~?. From the 
pressure data, we thus find that the equations of § 2 may be unreliable for 
meteors which penetrate below a height of about 75 km. 

If the reflection point lies in a region of the trail where «>2-4 x 10" the 
echo duration is given by (4:7). Taking D as proportional to p~', we find 
from (2.4) that 7 will be a maximum, 7,,, for a given trail when the 
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reflection point is at p= (3/2)Pmaxsi-€-, = (27/32) &max- Thus the duration of 
the echo from a point on the trail at pressure p 1s 


16 p \* 1 p if 
=e TS Se) ee eee ae 
eee Ge) ( ey Ae 
Using the previously determined values for D, etc., and (11.3) we obtain 
Tmn= 1X 10-10 A2 «8/3. cos?/? y/v. woe ee ee Lee 


It is assumed in the above that the only process reducing the electron 
density in the meteor trail is radial diffusion. 


§ 12. Turory oF THE Rapio Ecuo RATE 


The following summarizes the results of a theory (Kaiser, unpublished) 
which gives the observed echo rate in terms of the properties of the 
incident meteor flux and the radio equipment parameters. 

(a) Shower Meteors 

For the incident meteor flux defined by (9.1) the total echo rate in the 
sector 6—(6+dé) of the echo plane is V,d6, which is given in the notation 
of §9 by : 

Nie bhHI, ( cos 


(s—1)m,,° sin? x cos? 0,, \cos @,, 


where [,= v1*(x) dx,2-3(s—1-13)-12 ; h=mean height of radio 


)"se@, . (12.1) 


echo reflection points; m,=%(«,,uH)/(B cos x); bm,,1*/(s—1)=total 
incident flux of meteors with mass >m,, ; 0,, is the direction of the maxi- 
mum of the aerial polar diagram, S(@), in the echo plane, S(6,,)=1 3; %m is 
the minimum detectable electron line density at the mean height h in the 
direction 6,,,in the echo plane ; ¢;=(3s—7)/2,q,=2(s—1) fora«,,< 2-4 10; 
qi=6s—8, g2=8(s—1) for «,, >2°4X 1012, 


(b) Sporadic Meteors 


On the assumption of a uniform geocentric radiant distribution and an 
incident flux of meteors 


v,dm=bm~dm steradian—! em-2 sec-1, 
the following expression is obtained for the observed rate, N 6,9 dd dh 
of echoes with reflection points in the element —(6+-d¢d), #—(9+d9) : 


bhHI,I, Ccosdm 
$, 9= Gym a int gE 9), ST Lee} 


m 
subject to ¢ being sufficiently large that the range to height h may be 
written as h/sin ¢, which for h~100 km requires ¢>10°, 


nose (S08) (S22) SEN ag 


m 


71/2 


Le é cos dé; my,=2(a,,4H)/(B cos bm) ; 6, 9 are elevation and azi- 
muth respectively ; S(¢, ) represents the aerial polar diagram (amplitude) ; 


a, 48 the minimum detectable electron line density at range h/sin ¢,, in the 
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direction (¢,,,0). This latter is arbitrary, but may conveniently be taken as 
that for which | V4, 9d isa maximum, i.e.,¢,, is the optimum elevation. 


Since the spread of echo heights is fairly small compared with the 
mean height, the range distribution may be obtained from (12-3) by 
assuming a constant height, h, for all meteors. Thus the number of 
echoes in the range interval R—(R-++dR) is Nz dR where R ~h/sin ¢. 


woe OHI Wy 7008 ¢,, sin? db. * 
Ma E—Tpgt = wink cong | _, 2H 9) 40 
The finite spread of heights will cause (12.4) to be inaccurate when 
¢ approaches 90°, i.e., when Rh. 


(12.4) 


§ 13. Tue Mass DistRIBUTION OF THE INCIDENT METEORS 


From (12.1) and (12.2) we see that, provided the aerial polar diagram is 
maintained constant, the echo rate is proportional to hz,,"-8). Now «» 
is fixed by the sensitivity of the radio echo equipment, hence s may be 
determined by measuring the echo rate as a function of some convenient 
parameter of the equipment (e.g., transmitter power, receiver sensitivity, 
ete.). 

Considering only the case «,,< 2-4 101", we can, from (4.1) define an 
equipment sensitivity factor, F, by 


Hee PN \eg 2 ge eg kw tke os be 3(13.1) 


where €, is the minimum detectable echo power. Thus if we have two 
equipments, A and B, with aerials of identical polar diagram (and gain), 
but with sensitivities, /, and /’,, the echo rates N,, Nz, are related by 


N,|Ny»= (FE p> (h,[h,)e-?. 


The slight variation of h with equipment sensitivity may be estimated 
from the results given in § 9, leading to 


N,N p= (Fal Fs) 1+ 4(A/h) In(F/ Fy)... (13.2) 


The second factor in the right hand side of (13.2) is close to unity, hence 
we may use an approximate value for H/h which we will take as 0-07 
(h~100 km, H~7 km). 

McKinley (1951) has measured the sporadic meteor rate as a function of 
transmitter power (in one case with a slight change in wavelength). With 
values of F,/F, equal to 2-5, 4-0, and 43 he obtained NV,/N,=2-53, 3-79, 
and 45, respectively. The corresponding values obtained for s are 2-02, 
1:97, 2:02. Thus for sporadic meteors we will adopt s=2-00 + 0-02. 

To a good approximation we could have ignored the variation of h 
with equipment sensitivity and used N,/N,=(F,/F;,)** leading to 
s=1-99 + 0-02 which does not differ significantly from the above. 

McKinley’s observations were made with a horizontal half wave dipole 
aerial situated one-quarter wavelength above ground. The appropriate 
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value of ¢,, is 40° (see § 14) and the values of «,, for the various equipments 
are found to be 101, 2-5 102°, 4-0 101, 4-3x1011em~1. All of these 
values with the possible exception of the fourth (see §14) adequately 
satisfy the criterion «,< 2-4 x 10¥. 

Assuming that o,,,.—=101? for M=-++5, it follows from the above range 
of equipment sensitivities that the mass law is valid over the range of 
zenithal magnitudes +5 to +10. 

The inverse power law of meteor masses corresponds to a meteor 
magnitude distribution such that frequency increases by a constant ratio, 
a, per magnitude where 


8=1+2-5 logy, a. MEP, PRN ee ary MI 


Thus for sporadic meteors we obtain a=2-51+0-04. This is in good 
agreement with the most recent value given for telescopic meteors by 
Watson (1940), namely a=2-5, but differs considerably from the value 
a—=2-08 derived by Williams (1940). In contrast to the radio echo 
measurements the visual and telescopic observations require large 
corrections, including subjective factors, in converting observed to absolute 
rates, and values of a ranging from 2-0 to 4:0 have been quoted. 

McKinley has also measured the relative rates during the Quadrantid 
meteor shower, using the two equipments for which F,,/F,=—43. Owing 
to the high sporadic rate on the high power equipment, the ratio of the 
echo rates was observed to vary during the period of observation. De- 
pending upon how the correction is made for the sporadic background th 
mass law exponent, s, lies between 1-67 and 1-76. 

The value s=2-00 has been used ($10) to derive the atmospheric scale 
height from the sporadic meteor height distributions. Since the width 
of the height distribution for meteor showers is also a function only of H 
and s, the latter may be determined for the various showers. Preliminary 
measurements of the height distribution for several showers lead to the 
values of s and a given in the table. 


Shower 8 a 
Daytime Arietids, 1952 2-25+0:1 3-2 +0:3 
Perseids, 1952 1-56 +0-05 1-67 +0-07 
Geminids, 1951 and 1952 1-66 +0-04 1:84+0-07 
Quadrantids, 1953 1-68 +-0-06 1:87 -+0-10 


These values were obtained with an equipment for which «,,=2-3 x 104, 
hence they are valid for magnitudes in the vicinity of +6 to +7. The 
value for the Quadrantids is consistent with that estimated from McKinley’s 
results. 

With the exception of the Arietids the exponent s is considerably less 
than the sporadic value, i.e., the showers contain a relatively greater 
proportion of bright meteors, in agreement with telescopic and visual 
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observations. This can be explained if it is assumed that the Poynting- 
Robertson effect selectively removes the smaller particles from the streams. 
The daytime Arietids constitute an interesting exception. It is possible 
that this stream is of recent origin and that insufficient time has elapsed 
for the Poynting—Robertson effect to influence the mass distribution. 
However, this would appear to be contradicted by the diffuse character 
of the radiant (Aspinall and Hawkins 1951). Since the Arietids are a 
daytime stream, no visual data are available, however, it would appear 
that, although the shower is extremely active when observed by radio 
echo equipments («,,<10'), there would be relatively few bright visual 
meteors. 


§ 14. Taz Rapio Ecuo Rare as a FuncTION OF THE EQUIPMENT 
PARAMETERS 


In this section we will consider further the implications of the theoretical 
results given in §12. We will assume throughout that the mean height is 
effectively independent of «.,,. 


14.1. Aerial Polar Diagram and Gain 
(i) Sporadic Meteors 


The distribution of echoes in elevation and azimuth is obtained by 
putting s=2-0 in (12.2) and is a function of the aerial polar diagram. We 
will first consider the effect on the echo rate of the parameters of the 
aerial system. 

(a) Pencil Beam Aerial.—Provided the elevation of the narrow pencil 
beam is not too near 0° or 90°, the number of echoes per unit solid angle 
observed in the direction (¢, #) is simply proportional to S® (¢, #) and is a 
maximum in the direction of the beam axis. The total echo rate, N, is 
proportional to «,,~1 [f S% (¢, 3) dé d#, from which we may derive N as a 
function of aerial gain. 

m< 2:4 107?.—In this case the echo rate in the direction (4, #) is 
simply proportional to the aerial power polar diagram. Both «,, and 
ff S2(¢, 3) dé dd are inversely proportional to G, thus the echo rate is 
independent of the aerial gain. Qualitatively, this results from the fact 
that although a higher gain aerial sees fainter meteors, it covers a smaller 
solid angle. 

For constant aerial gain, the echo rate varies with ¢,, as cos ¢,, sin~1/"6,, 
provided ¢,, is greater than about 10°. The increase in rate as ¢,, 
decreases is due to the larger volume of the meteor ionization region 
intersected by the beam which more than counterbalances the increase in 
%m due to the greater range. As ¢,, approaches zero, the effect of the 
earth’s curvature must be taken into account and the rate obviously tends 
to a finite limit. 

om, >24X 10. Since q.=8, the variation of echo rate with gain will 
depend to some extent on the exact form of the polar diagram, but we will 
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obtain a reasonable estimate by taking S(¢, #) as a two-dimensional 
Gaussian function hence f[f S? (¢, 0) dé dd is proportional to G~*, a, is 
proportional to G4 and the total echo rate increases as Ge 

(b) Broad Beam Aerial.—We must now take account of the other factors 
in E (d, #) and the maximum echo rate does not necessarily occur in the 
direction of the aerial beam maximum. As an example the form of the 


elevation distribution, V,= if N 4, 9 3 has been evaluated for an aerial 
consisting of a horizontal half ‘wave dipole situated A/4 above a reflecting 
plane, and is given in fig. 21. The optimum elevation is 40° or 64 


depending on whether «,, is small or large compared with 2-4 10”. 


Fig. 21 


1°) 30 60 $0 
Differential elevation distribution, Ny, of sporadic meteor echoes observed 
with a A/2 dipole situated A/4 above ground (assuming a uniform radiant 
distribution). 
(2) On <2'4 x 1072 cm-"*, 
(b) om S> 2:4 x 102 em-, 


The range distribution for this simple aerial system has been evaluated 
from (12.4) assuming h=100 km and is given in fig. 22 together with the 
experimental results of McKinley (1951a). The theoretical curve has 
been normalized to represent the same number of echoes as observed. 
The discrepancy at short range is due to the assumption made in deducing 
the theoretical curve that all echoes are observed at constant height, 
otherwise the agreement is good. McKinley also gives results for a verti- 
cally directed Yagi aerial with a gain of twice the dipole system, but which 
gave an echo rate of only one half that observed on the dipole system. 
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The Yagi polar diagram is not given, however, a reasonable assumption 


leads to good agreement with the theory for the range distribution and the 
echo rate. 


(ii) Shower Meteors 


For an aerial which has a narrow beam in the echo plane, the echo rate 
varies with @ as S*8-)(@) for «,<2:4x102% and as S8°-0(6) for 
Gm o>2'4xX 101, Since for most showers s is less than 2-0, the echo rate 
in the former case falls off less rapidly than the power polar diagram 


of the aerial. 


Fig. 22 


——> Range, km. 
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Range distribution of Sporadic meteor echoes observed with 
a dipole aerial (see text). 
————— Theoretical, assuming constant height and a uniform geocentric 
radiant distribution. 


O Experimental, after McKinley (1951 a). 


The variation with gain cannot be uniquely specified, since it will 
depend on the manner in which the change in the gain is obtained. For 
instance, if the gain is increased by decreasing the beam width transverse 
to the echo plane while maintaining S(@) constant, the total echo rate will 
increase as GS-1 for «,,<2:4X 10" and as G46-) for a, >2-4x10¥. If 
the gain is increased by changing the beam width by equal amounts in, 
and transverse to, the echo plane, the rate, for a narrow Gaussian beam, 
will vary as Gs-*/ and G4s-*/2 for the respective cases, while if the beam 
width is changed in the echo plane only, the respective factors are G*? 


and. Gs. 
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14.2. Echo Rate as a Function of the Equipment Sensitivity Factor, F 

If the aerial polar diagram is maintained constant we find (§ 4), 

(a) a, ocF-1, hence N oc F8-1 for %_,<2°4X 10745 

(b) x, 0cF-4, hence N ocF**-) for a_>2:4X 10”. 
The echo rate will vary with F as shown in fig. 23 (with, of course, a 
smooth transition between the two regions), although we may in practice 
expect a somewhat greater number of enduring echoes («>2-4x 1053 
due to the distortion of the meteor trails by winds (see §6). It is clear 
that we should observe a rapid falling off in the number of echoes when «,, 
increases above about 1012 cm7}?. 


Fig. 23 
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——> -LOgi, (QF) 
=e. =[ ° l 
Echo rate, V, as a function of the equipment sensitivity factor, F (c, and c, 
are constants, s is the exponent in the meteor mass distribution). The 
experimental points are after McKinley (1951 a). (a), (b), (¢), A=9-22 m ; 
(d) A=9-55 m ; (e) A=5-35m; (f) A=2-83 m. 


The value of F=F, at the intersection of the two straight lines of 
fig. 23 will depend slightly on the form of the aerial polar diagram and on 
the value of s. Considering the case of meteor trails observed with a 
narrow Gaussian aerial beam, and neglecting the polarization effect we find 


(a) Shower meteors, s=1:5, a,=5°8 xX 1011 cm-}, 
s=2-0, %=4'8 x 1014 om-1. 
(b) Sporadic meteors, s—2-0, a,=4:8 x 1011 em-!, 


where «, is the value of «,, obtained from the Lovell-Clegg scattering 
formula (4.1) with F=F,. For sporadic meteors observed with a half 
wave dipole, one-quarter wavelength above ground we get «,=3-8 x 1011, 
where we have here taken «, as the minimum detectable line density, at 
range h (~100 km) in the zenith and related to F, by the Lovell-Clegg 
formula ; the corresponding minimum detectable line density at elevation 
Om=40° is 1014 cm-, 
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The above values have been used to fit the experimental results for 
sporadic meteors given by McKinley (1951 a) to the theoretical curve in 
fig. 23. The points a, b,c, d were obtained with equipments operating on a 
wavelength of 9m from which the results already discussed in §13 were 
derived ; e and f were obtained with equipments operating on wavelengths 
of 5-35 and 2-83 m respectively. The agreement for the 9m results ig 
excellent. While the agreement appears good for the 5-35 m equipment, 
these observations included a substantial percentage of echoes’ from 
Perseid meteors and if we correct for these (taking s = 1-56 for Perseids) the 
point would lie somewhat below the theoretical curve, as does the 2:83 m 
value. McKinley suggests that a proportion of the echoes may have 
originated from meteors which failed to reach the echo plane (corresponding 
to the individual meteor radiant) and since the ratio between head echo 
and body echo amplitudes is proportional to A1/2 this would act in favour 
of the longer wavelengths. On the other hand, such meteors would tend 
to be considerably larger on the average than those seen by specular 
reflection, and thus should be correspondingly fewer in number and it 
seems unlikely that this mechanism alone could account for the factor of 2 
discrepancy on 2:83m. An additional factor which operates against the 
shorter wavelength equipments is the effect of the finite meteor velocity 
discussed in §4. Although the diffusion coefficient may be somewhat less 
for the short wavelengths (since the meteor will be seen at slightly lower 
heights due to the reduced sensitivity) this will be more than counter- 
balanced by the A~*/? factor in 4, hence the echo amplitude for parallel - 
scattering, and the plasma resonance effect may be reduced. It is possible 
that the true echo rate—sensitivity factor relation passes somewhere between 
d and f, the former being somewhat enhanced by the plasma resonance 
and the latter reduced due to the larger value of 4; in any case both of 
these points are close to the transition region where a relatively small 
charge in F may produce a large change in echo rate. If the minimum 
detectable electron line densities for the various equipments have been 
estimated too low, the points d and f may be even closer to the bend in the 
theoretical curve. 

The echo amplitude distribution follows from the above. The rate, N,, 
of the echoes with amplitude greater than A is obtained by putting A 
proportional to «,//? (for a receiver with a linear response) ; thus A ocd! 
and fig. 23 also gives the form of log NV, vs. log A. Figure 24 gives log 
N, vs. log A for the echoes observed with parallel polarization by Closs et 
al. (1953) during the 1951 Arietids. The variation is of the predicted 
form, the change in gradient being of the order of 4 times. IPftwo straight 
lines are fitted to the points, this intersection should correspond to 
x =5>1011 thus the smallest detectable line density (corresponding to 
the smallest amplitude in fig. 24) is found to be about tas 1014 em—+, in 
good agreement with the values 1-65 x 10" and 1-5 101 calculated by 
Closs e¢ al. from the plasma resonance with transverse polarization and 
from the equipment parameters respectively. From the initial gradient 


534 T. R. Kaiser on 


in fig. 24 we find s=2-0 which is rather less than the value obtained in 
1952 from the Arietid height distribution (§13). The value of these 
particular results is somewhat limited due to the rather narrow range of 
amplitudes involved. Extension of the measurements toa wider range 
of amplitudes (say 100 to 1) requires a special radio receiving system but 
would provide an excellent method of studying the meteor mass 


distributions. 
Fig. 24 
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Rate, V4, of echoes with amplitude greater than A (Arietid 
meteor shower, 1951, A=4 m). 


§ 15. Tur DistRIBUTION OF RapIo EcHo DuRATIONS 


The theoretical differential duration distribution for echoes of the long 
enduring type, observed with an equipment sufficiently sensitive that 
Om < 2°4X 1022, ig 

Ve TT RETIN Sete CE Ade ee ee 
subject to the conditions that the effects of recombination, attachment 
and trail distortion by non-uniform winds are negligible. 

Figure 25 gives log v, vs. log + obtained with a 4m radio echo equipment 
during the Perseids (1950-1), Arietids (1951-2) and Geminids (1950-1). 
The straight lines give the theoretical gradients for the Perseids and Arietids 
obtained from (15.1) with the values of the exponent, s, given in § 11. 


(b) Perseids 
The fit between theory and experiment is reasonable up to r=30 sec. 


Beyond this the numbers are too few to determine the distribution, how- 
ever, an examination of the experimental results indicates that there were 
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far fewer than would be predicted by theory. In fact, there should be 
130 echoes with +>30 sec. whereas only 12 were observed. 


(6) Arietids 
The initial gradient defined by the experimental results agrees fairly 
well with the theory, but above r=9 sec the observed numbers fall off 
more rapidly than predicted. There should be 100 echoes with durations 
in excess of 10 seconds whereas only 4 were observed. 
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Differential distribution of echo durations measured with A=4 m (v, number of 
echoes per one second duration increment). 
e@—Arietids (1951-52). 
O—Perseids (1950-51). 
x—Geminids (1950-51). 
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(c) Geminids 

Tt would appear that the initial gradient defined by the experimental 
points is closer to that of the Arietids than would be expected from (15.1) 
with s=1-66 (see §11). It is possible that the Geminid mass distribution 
in the magnitude range responsible for echoes of several seconds duration 
differs from that for the smaller meteors from which the values s=1-66 was 
obtained ; on the other hand the duration, 2 sec, of the first group of 
Geminid meteors in fig. 25 was close to the limit of resolution of the 
apparatus used for these observations and the value of v, may be in error 
by.a greater amount than indicated by its standard deviation. In this 
case the rate falls off rapidly for durations above 5-6 sec. 
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The decrease in the Perseid echo rate for 7>30sec means either 
that the mass distribution cuts off above the mass corresponding, from 
(11.5), tO Omax=2 X 101 (i.e. magnitude —1) or else some mechanism other 
than radial diffusion is affecting the long echo durations, and similarly 
for the Arietids for r>9 sec. (corresponding to «,,,,>4x 101%). Ifwe rule 
out the possibility of such a drastic change in the mass law we are led to 
conclude that whatever process is operating depends on the meteor 
velocity (60 and 36 km sec~! for the Perseids and Arietids respectively), 
and thus must be pressure dependent since the higher the velocity the 
smaller is the pressure at which, on the average, an echo of given duration 
is obtained. The most probable cause would seem to be loss of electrons 
through attachment to neutral molecules, although it is possible that the 
disruption of the trails through atmospheric turbulence is less at greater 
heights. If attachment is responsible we can estimate the coefficient, 

'B,, from (7.2) ; putting 7=30, v=6 x 10° and 7=9, v=3-6 x 10°, we obtain 
B,=8-5x10-1* and 1:9x10-1}° respectively, which are of the expected 
order of magnitude. The maximum line densities, «,,.—=2x 10! and 
41018 for the Perseids and Arietids respectively, both occur (from 11.3) 
at a pressure of about 6 dynes, i.e., at a height of about 80 km, thus the 
estimates for 8, will apply to molecular rather than atomic oxygen. 


§ 16. Torat INcIDENT FLUX AND Space Density oF METEORS 


The relation between the radio echo rate and the total incident flux 
is obtained by integrating (12-1) or (12.2) for shower or sporadic meteors 
respectively. We will assume in this section that the radio echo equipment 
sensitivity is such that «,,< 2-4 1022. 


(a) Sporadic Meteors 
For a beamed aerial, elevated at an angle ¢,,, the total observed echo rate, 
N, is given, subject to the assumption of a uniform radiant distribution, 


by 
10chH b 


Vega (=| < one ON wees eee (16.1) 
where b/m,, steradian~! cm~? sec! is the incident flux of meteors of mass 
greater than m,, (defined in §12). The rate N is of course independent of 
G since m,,cG-} . 

The corresponding expression for the half wave dipole aerial system 
considered previously is (putting ¢,,—40°) 

N 2280 GH (Olina foc ea oe ed ay 


We would, in fact, be in error by less than a factor of two if we used 
(16-1) for the dipole aerial, taking G as the gain at elevation ¢,,. 
(b) Shower Meteors 


The total observed rate is obtained by integrating (12.1) with respect 
to@. It will be dependent on the zenith angle directly, as well as indirectly 
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through the variation of S(0) and 6,, as the echo plane passes through the 
aerial beam. 

We will evaluate the rate with the dipole aerial system for the simple 
case for which S(@) is symmetrical about =0 (e.g., at transit of the shower 
radiant for a beam directed north or south). To a sufficient approxima- 
tion we can write S(?)=cos [(7z/2) sin 0] and for s between 1-5 and 2 we 
obtain ae 9-541 b 
and bm,,‘1~®/(s—1) is the incident flux of meteors with masses greater than 
Mm 
The sporadic echo rate observed with the height finding equipment 
at Jodrell Bank is about 10 hr~! (averaged over 24 hours). The aerial 
gain is of the order of 10, ¢,,30 and «,,=2:3101!. Taking h=100 km, 
H=6-5 km, we find that the incident flux of meteors capable of producing 
line densities in excess of 2-3 10!! is 1-2 10-* steradian-1! km-? hr-!. 
Due to the actual concentration of geocentric radiants towards the apex 
of the earth’s way, these figures represent a mean over the whole sky. 
McKinley (1951 a) observes a midnight sporadic rate of 600-1000 hr with 
a half wave dipole aerial and a sensitivity such that «,,=10!° cm-!. 
Since the sporadic rate is a maximum at about 0600 and a minimum at 
about 1800 hr, we may take this as mean rate, hence from (16.2) the mean 
incident flux for «,,,, in excess of 101° cm~1is 3 to 5 x 10-? steradian~! km-? 
hr-!. The corresponding value, extrapolated from the Jodrell Bank 
result, is 3x 10-? ster.-! km~? hr“, in excellent agreement. Using this 
value, we find that a total of 1-4 10° such meteors are incident on the 
earth per day. 

Fletcher Watson (1941) from a study of visual and telescopic meteor 
rates, estimates that the number of meteors incident per day, brighter 
than visual magnitude +10, is 4-5 10°. This is in agreement with the 
value obtained from the radio echo rate if we assume that «,,,.—= 101% em~! 
corresponds to a meteor of magnitude about +4. The analysis of radio 
echo-visual correlations (§ 8) has already shown, however, that «,,,,=10" 
corresponds to a visual magnitude between +5 and +6; we will conse- 
quently adopt the former value and assume that the remaining discrepancy 
derives from the methods used in evaluating the incident flux from the 
visual and radio echo data. Thus, from the radio echo results, the mean 
incident sporadic flux of meteors with mass greater than m is 
3x 10-4 (m;/m) steradian-! km~? hr}, and the flux per day over the whole 
earth is 1-4 107 (m;/m), where m; is the mass of a 5th magnitude meteor. 

It has been shown in §11 that a 5th magnitude meteor weighs between 
10-3-10-4 g, thus the mass of meteoric dust incident on the earth per 
day lies between one and ten kilograms per visual magnitude, and the 
range over which this applies may extend from negative visual magnitudes 
to M—-+20 to +30 (beyond which the particles are probably removed 
from the vicinity of the earth’s orbit by solar radiation pressure). 
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The mean geocentric meteor velocity will be of the order of 40 km sec“? 
hence, using the values for the incident flux obtained above we obtain the 
space densities, P,,, P,, of meteoric particles with masses and radii 
greater than m and r respectively : 

P,,=3'5 X 10-3 (m,/m) cm? ; P,=3-5 x 10-8 (r;/r)? em? , (16.4) 
where 1; is the radius of a 5th magnitude meteor. 

The space density of dust particles has been estimated by various 
authors from observations of the zodiacal light. Van de Hulst (1947) 
gives r~0-035 cm for the particles responsible for the zodiacal light and 
deduces an integrated mass distribution P,—2%x10-*r-1® thus for 
r=0:035, P,=3x10-18. This radius is about the same as deduced for 
5th magnitude meteors, thus the density calculated by van de Hulst is 
about 10° times greater than that estimated from the meteor observations. 
The discrepancy between densities is less difficult to explain than that 
between the particle size distributions, since this would require the 
simultaneous coexistence of a large group of particles with P,art® 
whose members never strike the earth and a small one with P, ocr~3 (over 
the same range of radii) which is alone responsible for the meteor activity. 
Siedentopf, Behr, and Elsisser (1953) suggest that the zodiacal light is 
due to scattering particles of about 10-* cm radius with a space density 
of the order of 10-1° cm-, with P,,ocr-*°. This radius distribution is in 
fair agreement with that obtained for meteors. If we take r;=0-02 cm 
(§11) then the space density, extrapolated from the meteor data, of 
particles with radii in excess of 10-?cm is 3x 10-!%cm~-%, which still 
represents only a small fraction of the particles assumed to be responsible 
for the zodiacal light. A possible explanation is that there are a sub- 
stantial number of particles moving in approximately circular orbits 
which contribute to the zodiacal light but not to the meteor flux (see, 
however, Opik 1951). 

Radio echo rates have been measured for the Geminid and Quadrantid 
meteor streams on the height finding equipment at Jodrell Bank using 
an aerial system and a geometry which reasonably satisfies the assumptions 
made in (16.3). For the Geminid shower, the rate was about 200 hr-} 
for 730°, «,—2x10!1cm-!. The incident flux of Geminid meteors 
capable of producing electron densities in excess of this value is found from 
(16.3) to be 1-5 10-* km? hr~1, consequently the flux of masses greater 
than m is 5 x 10-8 (m;/m)"°* km~ hr-4, i.e., 1-4 107 (ms5/m)"®5 over the 
whole earth per day. The geocentric Geminid velocity is 36 km sec}, 
thus the space density for mass greater than m is 


Pn=4X 10-28 (m,;/m)*** cm-8. 
A similar calculation for the Quadrantid shower leads to results of the same 
order of magnitude. ; 
We see from the above that the space density and total incident rate on 
the earth for Geminids and sporadics are about equal for meteors brighter 
than 5th magnitude. Below zero magnitude there are 5 times as many 
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Geminids as Sporadics, while the converse holds for magnitudes less than 


+10. 


§ 17, Ionization PRopuUCED By METEORS 


Since the radio echo equipment measures the amount of ionization in 
the meteor trails, the totalerate of ionization may be estimated without 
introducing the electron line density—visual magnitude—mass relations. 

For sporadic meteors, we find that the mean rate of ion production is 
independent of height (over the range of heights covered by meteors in 
the range of sizes obeying the mass law with exponent s=2). It is given 
in terms of the echo rate and equipment parameters by 


27 
= (—-). eee ce: (17:1) 


4 cos ¢,, \m,, 


where b/m,, is related to the echo rate by ( 16.1) or (16.2) (for the aerial 
systems considered in §16). From §16 we have b/m,,=3x10-? 
steradian~! km~? hr-1 with «,,=—10!°cm—!, ¢,,=40°, thus g=2-3x 10-4 
em~> sec—!, 

If this rate of ionization were uniformly distributed, the resulting 
equilibrium volume density of ionization is n,=(q/«,)!/2 where «, is the 
effective recombination coefficient of the electrons. The volume density 
of meteor ionization will thus contribute significantly to the night time # 
region (n,~10*) only if «,~10-!2. Although this latter is of the order 
estimated during the initial decay of a discrete meteor echo (see §7), the 
effective value for the equilibrium distribution of ionization may be con- 
siderably greater (due to attachment of electrons, etc.). Itis not, however, 
likely to be greater than 10-8 (the observed ionospheric # region value), so 
that the contribution of the sporadic meteors to H region ionization will be 
between 10? and 104 electrons per cm’, i.e., between 1% and 100% of th 
observed night time density. : 

It is interesting to speculate on the character of the ionospheric layer 
that will be set up by the sporadic meteors. The condition that the 
production of electrons may be considered continuous over a large volume 
is that, on the average, an incident meteor will pass through the distri- 
bution of ionization remaining after the passage of a previous meteor 
before the total number of electrons produced by this latter has been 
significantly reduced below its initial value. If m, be the mass of that 
meteor which produces its maximum ionization, ,,,,=%, at pressure p, it 
can be shown that about 30° of the mean volume density of electrons at p 
is produced by meteors for which m<my,. In order to obtain a rough 
estimate of the condition for continuous electron production we will 
assume a mean meteor velocity of 40 km sec“? and neglect the variation 
- of «, with zenith distance. The rate of meteors with %px>%1 is propor- 
tional to «,~4, and for «,—101% is, (§ 16), 3 x 10~? steradian—! km~? hr-}, i.€., 
about 3X 10-25cm-2sec-! normal to the earth’s surface. Now «, is 
related to p by (11-3) thus the rate of such meteors at pressure p is of the 
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order of 4 10-16 p-3 and in time ¢ one meteor enters each area of 2:5 x 10% 
p> t-1. The crosssectional area of a meteor trail at time ¢ after the passage 
of the meteor is 47Dt=6 x 10° p-¥t (from (5.2)) and equating this to the 
above we find that the ionization at p may be regarded as continuous 
provided that the total number of electrons in a single trail is not reduced 
significantly in time 10° p?. This is likely tg occur through attachment, 
which, from (3.13) requires a time of at least t~(8,n,,)~1 to be effective. 
Assuming that attachment takes place to oxygen (atoms or molecules) and 
taking 7’=300°, this becomes t=2 x 10~18 (8,p)~! and equating it to 10°p? 
we find that the lower limit to a continuous band of ionization will be 
at p~10-* B,-18 or, putting B,~10-' to 10-1%, p~10-1. The upper 
height limit will depend upon the size of the smallest meteor capable of 
producing ionization (assuming the mass law to hold down to this limit). 
Whipple (1952) suggests that this is of the order of a few microns which, 
using the estimates of meteor size given in § 11, corresponds to «,..~10°® 
ie., from (11.3) p22x10-dyne cm~*. Using the data given by the 
Rocket Panel (1952) these pressures define the height range 115-130 km. 

Summarizing, we would expect the sporadic meteor flux to produce a 
continuous band of ionization with a volume electron density in excess of 
10? em-*, extending from about 130 down to 115 km, below which the 
individual meteors will not join up to form a continuum and the ionization 
will be distributed in discrete patches. Apart from the uncertainty in the 
volume electron density, this description corresponds remarkably to the 
night time # region of the ionosphere, the lower parts of which give rise to 
the ‘ Sporadic #’’ reflections observed at long wavelengths. 

A similar calculation to the above for a shower such as the Geminids, 

assuming s=2 and the incident flux determined in §16 would give 
q=7x 10-4 cm? sec“1, i.e., of the same order as for sporadic meteors. 
However, all of the meteor showers with the possible exception of the 
Arietids have s<2, i.e., contain relatively fewer small meteors than in 
the sporadic flux and hence would contribute to a much smaller extent to 
the continuous band of ionization. In general, therefore, we might 
expect during a meteor shower to observe an increase in the incidence of 
sporadic # reflection, but negligible contribution to the normal E region 
ionization. 
Experimental evidence indicating a correlation between the ionization 
in the # region and meteor activity has been summarized by Lovell 
(1948). While little evidence exists in the case of the normal F regions, 
the seasonal variations of sporadic ionization shows a correlation with that 
of the transient echoes, while during certain intense meteor shower 
activity, notably the Leonids, 1933, 1936, and the Giacobinids, 1946, a 
considerable increase in the occurrence of sporadic EZ ionization was 
observed. In the latter case, Lovell et al. (1947) observed a maximum echo 
rate of 10* hr~? on an equipment with a normal sporadic rate of 2 hr-1 and 
Pierce (1947) indicates that these meteors produced a reflection layer which 
persisted for several hours on a frequency of 3-5 Mc/s. 
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§ 18. SumMARY 


The development of a satisfactory theory of the scattering of radio 
waves from meteor trails, with the previously established theory of 
evaporation from meteors, makes possible a detailed interpretation of the 
radio echoes in terms of the properties of the incident meteors, of the 
atmosphere in the meteor ionization region, and of the processes connected 
with the production and decay of ionization in a meteor trail. 

The small value obtained for the electronic recombination coefficient 
in meteor trails suggests that dissociative recombination, which is probably 
responsible for the high recombination rate in the E region of the 
ionosphere, is not important, from which it is concluded that the ions are 
predominantly atomic due to ionization of the evaporated meteor atoms. 
The high value (0-1— 1) obtained for the probability, 8, that an evaporated 
meteor atom will produce a free electron in subsequent collisions with air 
molecules is at first sight surprising, as is also its relative independence of 
velocity over the range 20-60km/sec (100-1000ev). A plausible 
explanation is that the cross section for single ionization of several or all 
of the meteor constituents is large compared with the momentum transfer 
cross section over the whole velocity range and that all other ionization 
cross sections are small (since the first ionization potential of a meteor 
is much smaller than its second and than that of the atmospheric 
molecules). In this case the ionization probability for the constituents 
concerned will be unity and velocity independent (even though the cross 
section may vary) and the measured probability will just represent their 
relative abundance. 

An experimental study of the decay type echoes supports the theoretical 
prediction that under certain circumstances the scattering of a wave with 
its incident electric vector transverse to the meteor trail is enhanced by a 
plasma resonance phenomenon. From the decay of these echoes, the 
collision diameter for the ion-neutral molecule collisions is estimated as 
8-5 x 10-8 cm, i.e., several times greater than the classical kinetic theory 
value (~3 x 10-8) for collisions between neutral molecules. A study of the 
durations of enduring echoes indicates that some process other than radial 
diffusion is effective in decreasing the volume density of ionization in the 
trails formed by the brighter meteors. It is suggested that this may be 
due to the attachment of electrons to oxygen molecules and the coefficient 
of attachment is estimated to be of the order of 10-4 cm? sec™?. 

The maximum electron line density, «,,,x, produced along a meteor trail 
plays the same role in the radio echo studies as does the equivalent stellar 
magnitude in visual and telescopic observations. The relation between 
these parameters may depend on the meteor velocity and is given by 
(8.4) for v=60kmsec. It is desirable that simultaneous visual-radio 
echo observations be carried out during the activity of major showers 
other than the Perseids with a view to investigating the velocity depend- 
ence of the ratio (luminous efficiency)/(probability of ionization). 
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The radio echo technique is now established as a powerful tool for 
investigation of the atmosphere in the region 80-110 km. Measurements 
of the Doppler effect due to line of sight motion of the meteor trails have 
been interpreted in terms of the upper atmospheric winds, while a recently 
developed theory of the radio echo height distribution makes possible the 
independent measurement of scale height and pressure. Preliminary 
values of these quantities obtained at Jodrell Bank (latitude 53°) are close 
to those from rocket observations at latitude 32°, whereas the pressure 
values obtained from visual meteor observations at Massachussetts 
(latitude 42°) are several times higher. Some 2-300 echoes are necessary 
to obtain reasonable accuracy in the pressure and scale height. Assuming 
that the velocity can be measured with sufficient precision for only 10% 
of the echoes, this would require about 10 days continuous operation with 
an equipment giving a mean sporadic rate of 10 hr-t.. In spite of the 
restricted height range of the radio echo measurements, the relative 
simplicity of the technique recommends it as a means for studying the 
circulation in the upper atmosphere through the establishment of wind 
measuring and height finding equipments at a number of suitable locations. 

The distribution of sporadic meteor masses and the total incident 
flux has been determined from the observed echo rates. The latter leads 
to an estimate of the space density of meteoric particles in the vicinity 
of the earth which is insufficient to account for the zodiacal light. One 
explanation is that a large proportion of the interplanetary particles 
may move in approximately circular orbits, alternatively particles of 
about 10 microns radius may be about 1000 times more frequent than would 
be predicted by extrapolation of the radio echo measurements to this 
meteor size. The mass distributions have been obtained from the height 
measurements for a number of meteor showers and in all cases except one _ 
these have values of the exponent s which are less than 2-0, indicating that 
there are relatively fewer small meteors than in the sporadic distribution. 
A consequence is that in comparing the echo rates of the various showers 
and of the sporadic meteors, account must be taken of the difference in the 
mass distributions, the equipment sensitivity, and the properties of the 
aerial system. The mean heights obtained for the various showers, when 
suitably corrected for the geometry of the particular aerial system used, 
agree fairly well with those obtained for sporadic meteors of the same 
velocity, thus there is as yet no evidence on these grounds for a difference 
in composition between shower and sporadic meteors. On the other hand 
the mean height may be fairly insensitive to the composition, since B 
enters (2.5) only as the cube root, and the product lp?/> may not vary 
significantly between stone and iron meteors. 

The estimated incident flux of sporadic meteors is sufficient to maintain 
a continuous ionospheric layer extending from about 130 down to 115 km, 
below which the ionization will be distributed in patches becoming fewer 
(and larger) with decreasing height. The electron density in the continuum 
will be in excess of 10? cm~!, however, if the second explanation given above 
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for the discrepancy between the meteor and zodiacal light estimates of the 
space density of interplanetary dust is correct, the rate of incidence of 
particles too small to be detected by radio echo or telescopic observation 
would be sufficient to maintain the night time H region even if the effective 
recombination coefficient is as high as 10-’ em’ sec~!. 

It will be appreciated that many of the previously unpublished results 
(especially in § § 9-17) are of apreliminary nature and indicate the direction 
which future work will take. It is intended that the theory of the meteor 
height distribution and echo rate, and the experimental determination of 


atmospheric pressure, scale height, etc., should be published in greater 
detail in the near future. 
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